Hydrodynamical simulations of galaxy formation and evolution attempt to fully model the physics that shapes galaxies. The agreement between the morphology of simulated and real galaxies, and the way the morphological types are distributed across galaxy scaling relations are important probes of our knowledge of galaxy formation physics. Here we propose an unsupervised deep learning approach to perform a stringent test of the fine morphological structure of galaxies coming from the Illustris and IllustrisTNG (TNG100 and TNG50) simulations against observations from a subsample of the Sloan Digital Sky Survey. Our framework is based on PixelCNN, an autoregressive model for image generation with an explicit likelihood. We adopt a strategy that combines the output of two PixelCNN networks in a metric that isolates the fine morphological details of galaxies from the sky background. We are able to \emph{quantitatively} identify the improvements of IllustrisTNG, particularly in the high-resolution TNG50 run, over the original Illustris. However, we find that the fine details of galaxy structure are still different between observed and simulated galaxies. This difference is driven by small, more spheroidal, and quenched galaxies which are globally less accurate regardless of resolution and which have experienced little improvement between the three simulations explored. We speculate that this disagreement, that is less severe for quenched disky galaxies, may stem from a still too coarse numerical resolution, which struggles to properly capture the inner, dense regions of quenched spheroidal galaxies.

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.

A significant fraction of high redshift star-forming disc galaxies are known to host giant clumps, whose formation, nature and role in galaxy evolution are yet to be understood. In this work we first present a new automated method based on deep neural networks to detect clumps in galaxy images and show that it is more sensitive and faster than previous proposed methods. We then use this method to systematically detect clumps in the rest-frame optical and UV images of a complete sample of $\sim1500$ star forming galaxies at $1<z<3$ in the CANDELS survey as well as in images from the VELA forward modeled zoom-in cosmological simulations. We show that observational effects have a dramatic impact on the derived clump properties leading to an overestimation of the clump mass up to a factor of 10, which highlights the importance of fair and direct comparisons between observations and simulations. After correcting for these effects with a simple mixture density network, we estimate that the clump stellar mass function follows a power-law down to the completeness limit ($10^{7}$ solar masses) with the majority of the clumps being less massive than $10^9$ solar masses. We find that between $20-40\%$ of $\sim10^{10}$ solar mass star-forming galaxies contain at least one massive ($M_c>10^{7}M_\odot$) off-centered clump but only $\sim2-5\%$ of the total galaxy stellar mass is in those clumps. We also show indications that the contribution of clumps to the stellar mass is more important in extended and low mass galaxies. The simulations explored in this work overall reproduce the shape of the observed clump stellar mass function when confronted under the same conditions although tend to lie in the lower limit of the confidence intervals of the observations. This agreement suggests that most of the observed clumps are formed in-situ through violent disk instabilities.

Using analytic modeling and simulations, we address the origin of an abundance of star-forming, clumpy, extended gas rings about massive central bodies in massive galaxies at $z \!<\! 4$. Rings form by high-angular-momentum streams and survive in galaxies of $M_{\rm star} \!>\! 10^{9.5-10} M_\odot$ where merger-driven spin flips and supernova feedback are ineffective. The rings survive after events of compaction to central nuggets. Ring longevity was unexpected based on inward mass transport driven by torques from violent disc instability. However, evaluating the torques from a tightly wound spiral structure, we find that the timescale for transport per orbital time is long and $\propto\! \delta_{\rm d}^{-3}$, with $\delta_{\rm d}$ the cold-to-total mass ratio interior to the ring. A long-lived ring forms when the ring transport is slower than its replenishment by accretion and the interior depletion by SFR, both valid for $\delta_{\rm d} \!<\! 0.3$. The central mass that lowers $\delta_{\rm d}$ is a compaction-driven bulge and/or dark matter, aided by the lower gas fraction at $z \!<\! 4$, provided that it is not too low. The ring is Toomre unstable for clump and star formation. The high-$z$ dynamic rings are not likely to arise form secular resonances or collisions. AGN feedback is not expected to affect the rings. Mock images of simulated rings through dust indicate qualitative consistency with observed rings about bulges in massive $z\!\sim\!0.5\!-\!3$ galaxies, in $H_{\alpha}$ and deep HST imaging. ALMA mock images indicate that $z\!\sim\!0.5\!-\!1$ rings should be detectable. We quote expected observable properties of rings and their central nuggets.

Existing models of galaxy formation have not yet explained striking correlations between structure and star-formation activity in galaxies, notably the sloped and moving boundaries that divide star-forming from quenched galaxies in key structural diagrams. This paper uses these and other relations to ``reverse-engineer'' the quenching process for central galaxies. The basic idea is that star-forming galaxies with larger radii (at a given stellar mass) have lower black-hole masses due to lower central densities. Galaxies cross into the green valley when the cumulative effective energy radiated by their black hole equals $\sim4\times$ their halo-gas binding energy. Since larger-radii galaxies have smaller black holes, one finds they must evolve to higher stellar masses in order to meet this halo-energy criterion, which explains the sloping boundaries. A possible cause of radii differences among star-forming galaxies is halo concentration. The evolutionary tracks of star-forming galaxies are nearly parallel to the green-valley boundaries, and it is mainly the sideways motions of these boundaries with cosmic time that cause galaxies to quench. BH-scaling laws for star-forming, quenched, and green-valley galaxies are different, and most BH mass growth takes place in the green valley. Implications include: the radii of star-forming galaxies are an important second parameter in shaping their black holes; black holes are connected to their halos but in different ways for star-forming, quenched, and green-valley galaxies; and the same BH-halo quenching mechanism has been in place since $z \sim 3$. We conclude with a discussion of black hole-galaxy co-evolution, the origin and interpretation of BH scaling laws.

We predict, analytically and by simulations, that gas discs tend to survive only in haloes above a threshold mass $\sim 2 \times 10^{11} M_\odot$ (stellar mass $\sim 10^9 M_\odot$), with only a weak redshift dependence. At lower masses, the disc spins typically flip in less than an orbital time due to mergers associated with a change in the pattern of the feeding cosmic-web streams. This threshold arises from the halo merger rate when accounting for the mass dependence of the ratio of galactic baryons and halo mass. Above the threshold, wet compactions lead to massive central nuggets that allow the longevity of extended clumpy gas rings. Supernova feedback has a major role in disrupting discs below the critical mass, by driving the stellar-to-halo mass ratio that affects the merger rate, by stirring up turbulence and suppressing high-angular-momentum gas supply, and by confining major compactions to the critical mass. Our predictions seem consistent with current observed fractions of gas discs, to be explored by future observations that will resolve galaxies below $10^9 M_\odot$ at high redshifts, e.g. by JWST.

This paper studies pseudo-bulges (P-bulges) and classical bulges (C-bulges) in Sloan Digital Sky Survey central galaxies using the new bulge indicator $\Delta\Sigma_1$, which measures relative central stellar-mass surface density within 1 kpc. We compare $\Delta\Sigma_1$ to the established bulge-type indicator $\Delta\langle\mu_e\rangle$ from Gadotti (2009) and show that classifying by $\Delta\Sigma_1$ agrees well with $\Delta\langle\mu_e\rangle$. $\Delta\Sigma_1$ requires no bulge-disk decomposition and can be measured on SDSS images out to $z = 0.07$. Bulge types using it are mapped onto twenty different structural and stellar-population properties for 12,000 SDSS central galaxies with masses 10.0 < log $M_*$/$M_{\odot}$ < 10.4. New trends emerge from this large sample. Structural parameters show fairly linear log-log relations vs. $\Delta\Sigma_1$ and $\Delta\langle\mu_e\rangle$ with only moderate scatter, while stellar-population parameters show a highly non-linear "elbow" in which specific star-formation rate remains roughly flat with increasing central density and then falls rapidly at the elbow, where galaxies begin to quench. P-bulges occupy the low-density end of the horizontal arm of the elbow and are universally star-forming, while C-bulges occupy the elbow and the vertical branch and exhibit a wide range of star-formation rates at fixed density. The non-linear relation between central density and star-formation rate has been seen before, but this mapping onto bulge class is new. The wide range of star-formation rates in C-bulges helps to explain why bulge classifications using different parameters have sometimes disagreed in the past. The elbow-shaped relation between density and stellar indices suggests that central structure and stellar-populations evolve at different rates as galaxies begin to quench.

As part of the AGORA High-resolution Galaxy Simulations Comparison Project (Kim et al. 2014, 2016) we have generated a suite of isolated Milky Way-mass galaxy simulations using 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical galaxy formation community. In these simulations we adopted identical galactic disk initial conditions, and common physics models (e.g., radiative cooling and ultraviolet background by a standardized package). Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production were carefully constrained. Here we release the simulation data to be freely used by the community. In this release we include the disk snapshots at 0 and 500Myr of evolution per each code as used in Kim et al. (2016), from simulations with and without star formation and feedback. We encourage any member of the numerical galaxy formation community to make use of these resources for their research - for example, compare their own simulations with the AGORA galaxies, with the common analysis yt scripts used to obtain the plots shown in our papers, also available in this release.

We use bulge-type classifications of 809 representative SDSS galaxies by Gadotti (2009) to classify a large sample of galaxies into real bulges (classical or elliptical) and pseudobulges using Random Forest. We use structural and stellar population predictors that can easily be measured without image decomposition. Multiple parameters such as the central mass density with 1 kpc, concentration index, S\'{e}rsic index and velocity dispersion do result in accurate bulge classifications when combined together. We classify $\sim 44,500$ face-on galaxies above stellar mass of 10$^{10}$ M$_\odot$ and redshift $ 0.02 < z < 0.07$ into real bulges or pseudobulges with $93 \pm 2$\% accuracy. We show that $\sim 75 - 90\%$ of AGNs identified by the optical line ratio diagnostic are hosted by real bulges. The pseudobulge fraction significantly decreases with AGN signature as the line ratios change from indicating pure star formation ($\sim 54 \pm 4$ \%), to composite of star formation and AGN ($\sim 18 \pm 3$\%), and to AGN-dominated galaxies ($\sim 5 \pm 3$\%). Using the dust-corrected [\ion{O}{3}] luminosity as an AGN accretion indicator, and the stellar mass and radius as proxies for a black hole mass, we find that AGNs in real bulges have lower Eddington ratios than AGNs in pseudobulges. Real bulges have a wide range of AGN and star formation activities, although most of them are weak AGNs. For both bulge-types, their Eddington ratios are correlated with specific star formation rates (SSFR). Real bulges have lower specific accretion rate but higher AGN fraction than pseudobulges do at similar SSFRs.

As computational resolution of modern cosmological simulations reach ever so close to resolving individual star-forming clumps in a galaxy, a need for "resolution-appropriate" physics for a galaxy-scale simulation has never been greater. To this end, we introduce a self-consistent numerical framework that includes explicit treatments of feedback from star-forming molecular clouds (SFMCs) and massive black holes (MBHs). In addition to the thermal supernovae feedback from SFMC particles, photoionizing radiation from both SFMCs and MBHs is tracked through full 3-dimensional ray tracing. A mechanical feedback channel from MBHs is also considered. Using our framework, we perform a state-of-the-art cosmological simulation of a quasar-host galaxy at z~7.5 for ~25 Myrs with all relevant galactic components such as dark matter, gas, SFMCs, and an embedded MBH seed of ~> 1e6 Ms. We find that feedback from SFMCs and an accreting MBH suppresses runaway star formation locally in the galactic core region. Newly included radiation feedback from SFMCs, combined with feedback from the MBH, helps the MBH grow faster by retaining gas that eventually accretes on to the MBH. Our experiment demonstrates that previously undiscussed types of interplay between gas, SFMCs, and a MBH may hold important clues about the growth and feedback of quasars and their host galaxies in the high-redshift Universe.

This paper is one in a series that explores the importance of radius as a second parameter in galaxy evolution. The topic investigated here is the relationship between star formation rate (SFR) and galaxy radius ($R_{\rm e}$) for main-sequence star-forming galaxies. The key observational result is that, over a wide range of stellar mass and redshift in both CANDELS and SDSS, there is little trend between SFR and $R_{\rm e}$ at fixed stellar mass. The Kennicutt-Schmidt law, or any similar density-related star formation law, then implies that smaller galaxies must have lower gas fractions than larger galaxies (at fixed $M_{\ast}$), and this is supported by observations of local star-forming galaxies. We investigate the implication by adopting the equilibrium "bathtub" model: the ISM gas mass is assumed to be constant over time and the net star formation rate is the difference between the accretion rate of gas onto the galaxy from the halo and the outflow rate due to winds. To match the observed null correlation between SFR and radius, the bathtub model requires that smaller galaxies at fixed mass have weaker galactic winds. Our hypothesis is that galaxies are a 2-dimensional family whose properties are set mainly by halo mass and concentration. Galaxy radius and accretion rate plausibly both depend on halo concentration, which predicts how wind strength should vary with $R_{\rm e}$ and SFR.

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.

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.

Theory suggests that there are two primary modes of accretion through which dark matter halos acquire the gas to form and fuel galaxies, hot and cold mode accretion. In cold mode accretion, gas streams along cosmic web filaments to the center of the halo, allowing for the efficient delivery of star-forming fuel. Recently, two QSO-illuminated HI Lyman alpha (Ly{\alpha}) emitting objects were reported to have properties of cold, rotating structures (Martin et al. 2015, Martin et al. 2016). However, the spatial and spectral resolution available was insufficient to constrain radial flows associated with connecting filaments. With the Keck Cosmic Web Imager (KCWI) we now have eight times the spatial resolution, permitting the detection of these in-spiraling flows. In order to detect these inflows, we introduce a suite of models which incorporate zonal radial flows, demonstrate their performance on a numerical simulation that exhibits coldflow accretion, and show that they are an excellent match to KCWI velocity maps of two Ly{\alpha} emitters observed around high-redshift quasars. These Multi-Filament Inflow models kinematically isolate zones of radial inflow that correspond to extended filamentary emission. The derived gas flux and inflow path is sufficient to fuel the inferred central galaxy star formation rate and angular momentum. Thus, our kinematic emission maps provide strong evidence for the inflow of gas from the cosmic web building galaxies at the peak of star formation.

The role of major mergers in galaxy evolution remains a key open question. Existing empirical merger identification methods use non-parametric and subjective visual classifications which can pose systematic challenges to constraining merger histories. As a first step towards overcoming these challenges, we develop and share publicly a new Python-based software tool that identifies and extracts the flux-wise and area-wise significant contiguous regions from the model-subtracted "residual" images produced by popular parametric light-profile fitting tools (e.g., GALFIT). Using Hubble Space Telescope ($HST$) $H$-band single-S\'ersic residual images of $17$ CANDELS galaxies, we demonstrate the tool's ability to measure the surface brightness and improve the qualitative identification of a variety of common residual features (disk structures, spiral substructures, plausible tidal features, and strong gravitational arcs). We test our method on synthetic $HST$ observations of a $z\sim 1.5$ major merger from the VELA hydrodynamic simulations. We extract $H$-band residual features corresponding to the birth, growth, and fading of tidal features during different stages and viewing orientations at CANDELS depths and resolution. We find that the extracted features at shallow depths have noisy visual appearance and are susceptible to viewing angle effects. For a VELA $z\sim 3$ major merger, we find that James Webb Space Telescope NIRCam observations can probe high-redshift tidal features with considerable advantage over existing $HST$ capabilities. Further quantitative analysis of plausible tidal features extracted with our new software hold promise for the robust identification of hallmark merger signatures and corresponding improvements to merger rate constraints.

We address a simple model where the Kennicutt-Schmidt (KS) relation between the macroscopic densities of star-formation rate (SFR, $\rho_{\rm sfr}$) and gas ($n$) in galactic discs emerges from self-regulation of the SFR via supernova feedback. It arises from the physics of supernova bubbles, insensitive to the microscopic SFR recipe and not explicitly dependent on gravity. The key is that the filling factor of SFR-suppressed supernova bubbles self-regulates to a constant, $f\sim 0.5$. Expressing the bubble fading radius and time in terms of $n$, the filling factor is $f \propto S\,n^{-s}$ with $s\sim 1.5$, where $S$ is the supernova rate density. A constant $f$ thus refers to $\rho_{\rm sfr} \propto n^{1.5}$, with a density-independent SFR efficiency per free-fall time $\sim 0.01$. The self-regulation to $f \sim 0.5$ and the convergence to a KS relation independent of the local SFR recipe are demonstrated in cosmological and isolated-galaxy simulations using different codes and recipes. In parallel, the spherical analysis of bubble evolution is generalized to clustered supernovae, analytically and via simulations, yielding $s \simeq 1.5 \pm 0.5$. An analysis of photo-ionized bubbles about pre-supernova stars yields a range of KS slopes but the KS relation is dominated by the supernova bubbles. Superbubble blowouts may lead to an alternative self-regulation by outflows and recycling. While the model is over-simplified, its simplicity and validity in the simulations may argue that it captures the origin of the KS relation.

The majority of massive star-forming galaxies at $z\sim2$ have velocity gradients suggestive of rotation, in addition to large amounts of disordered motions. In this paper, we demonstrate that it is challenging to distinguish the regular rotation of a disk galaxy from the orbital motions of merging galaxies with seeing-limited data. However, the merger fractions at $z\sim2$ are likely too low for this to have a large effect on measurements of disk fractions. To determine how often mergers pass for disks, we look to galaxy formation simulations. We analyze $\sim$24000 synthetic images and kinematic maps of 31 high-resolution simulations of isolated galaxies and mergers at $z\sim2$. We determine if the synthetic observations pass criteria commonly used to identify disk galaxies, and whether the results are consistent with their intrinsic dynamical states. Galaxies that are intrinsically mergers pass the disk criteria for anywhere from 0 to 100$\%$ of sightlines. The exact percentage depends strongly on the specific disk criteria adopted, and weakly on the separation of the merging galaxies. Therefore, one cannot tell with certainty whether observations of an individual galaxy indicate a merger or a disk. To estimate the fraction of mergers passing as disks in current kinematics samples, we combine the probability that a merger will pass as a disk with theoretical merger fractions from a cosmological simulation. Taking the latter at face-value, the observed disk fractions are overestimated by small amounts: at most by $5\%$ at high stellar mass ($10^{10-11}$ M$_{\odot}$) and $15\%$ at low stellar mass ($10^{9-10}$ M$_{\odot}$).

We study the rest-frame optical mass-size relation of bulges and discs from z~2 to z~0 for a complete sample of massive galaxies in the CANDELS fields using 2 component S\'ersic models (Dimauro et a. 2018). Discs and star forming galaxies follow similar mass-size relations. The mass-size relation of bulges is less steep than the one of quiescent galaxies (best fit slope of ~ 0.7 for quiescent galaxies against ~ 0.4 for bulges). We find little dependence of the structural properties of massive bulges and discs with the global morphology of galaxies (disc vs. bulge dominated) and the star formation activity (star-forming vs. quiescent). This result suggests similar bulge formation mechanisms for most massive galaxies and also that the formation of the bulge component does not significantly affect the disc structure. Our findings pose a challenge to models envisioning multiple channels for massive bulge growth, such as disc instabilities and mergers.

We study the components of cool and warm/hot gas in the circumgalactic medium (CGM) of simulated galaxies and address the relative production of OVI by photoionization versus collisional ionization, as a function of halo mass, redshift, and distance from the galaxy halo center. This is done utilizing two different suites of zoom-in hydro-cosmological simulations, VELA (6 halos; $z>1$) and NIHAO (18 halos; to $z=0$), which provide a broad theoretical basis because they use different codes and physical recipes for star formation and feedback. In all halos studied in this work, we find that collisional ionization by thermal electrons dominates at high redshift, while photoionization of cool or warm gas by the metagalactic radiation takes over near $z\sim2$. In halos of $\sim 10^{12}M_{\odot}$ and above, collisions become important again at $z<0.5$, while photoionization remains significant down to $z=0$ for less massive halos. In halos with $M_{\textrm v}>3\times10^{11}~M_{\odot}$, at $z\sim 0$ most of the photoionized OVI is in a warm, not cool, gas phase ($T\lesssim 3\times 10^5$~K). We also find that collisions are dominant in the central regions of halos, while photoionization is more significant at the outskirts, around $R_{\textrm v}$, even in massive halos. This too may be explained by the presence of warm gas or, in lower mass halos, by cool gas inflows.

We model the projected b/a-log a distributions of CANDELS main sequence star-forming galaxies, where a (b) is the semi-major (semi-minor) axis of the galaxy images. We find that smaller-a galaxies are rounder at all stellar masses M and redshifts, so we include a when analyzing b/a distributions. Approximating intrinsic shapes of the galaxies as triaxial ellipsoids and assuming a multivariate normal distribution of galaxy size and two shape parameters, we construct their intrinsic shape and size distributions to obtain the fractions of prolate, oblate and spheroidal galaxies in each redshift and mass bin. We find that galaxies tend to be prolate at low m and high redshifts, and oblate at high M and low redshifts, qualitatively consistent with van der Wel et al. (2014), implying that galaxies tend to evolve from prolate to oblate. These results are consistent with the predictions from simulations (Ceverino et al. 2015, Tomassetti et al. 2016) that the transition from prolate to oblate is caused by a compaction event at a characteristic mass range, making the galaxy center baryon dominated. We give probabilities of a galaxy's being prolate, oblate or spheroidal as a function of its M, redshift, projected b/a and a, which can facilitate target selections of galaxies with specific shapes at hight redshifts. We also give predicted optical depths of galaxies, which are qualitatively consistent with the expected correlation that AV should be higher for edge-on disk galaxies in each log a slice at low redshift and high mass bins.

The similarity between the distributions of spins for galaxies ($\lambda_{\rm g}$) and for dark-matter haloes ($\lambda_{\rm h}$), indicated both by simulations and observations, is naively interpreted as a one-to-one correlation between the spins of a galaxy and its host halo. This is used to predict galaxy sizes in semi-analytic models via $R_{\rm e}\simeq\lambda_{\rm h} R_{\rm v}$, with $R_{\rm e}$ the half-mass radius of the galaxy and $R_{\rm v}$ the halo radius. Utilizing two different suites of zoom-in cosmological simulations, we find that $\lambda_{\rm g}$ and $\lambda_{\rm h}$ are in fact only barely correlated, especially at $z\geq 1$. A general smearing of this correlation is expected based on the different spin histories, where the more recently accreted baryons through streams gain and then lose significant angular momentum compared to the gradually accumulated dark matter. Expecting the spins of baryons and dark matter to be correlated at accretion into $R_{\rm v}$, the null correlation at the end reflects an anti-correlation between $\lambda_{\rm g}/\lambda_{\rm h}$ and $\lambda_{\rm h}$, which can partly arise from mergers and a compact star-forming phase that many galaxies undergo. On the other hand, the halo and galaxy spin vectors tend to be aligned, with a median $\cos\theta=0.6$-0.7 between galaxy and halo, consistent with instreaming within a preferred plane. The galaxy spin is better correlated with the spin of the inner halo, but this largely reflects the effect of the baryons on the halo. Following the null spin correlation, $\lambda_{\rm h}$ is not a useful proxy for $R_{\rm e}$. While our simulations reproduce a general relation of the sort $R_{\rm e}=AR_{\rm vir}$, in agreement with observational estimates, the relation becomes tighter with $A=0.02(c/10)^{-0.7}$, where $c$ is the halo concentration, which in turn introduces a dependence on mass and redshift.

We use machine learning to identify in color images of high-redshift galaxies an astrophysical phenomenon predicted by cosmological simulations. This phenomenon, called the blue nugget (BN) phase, is the compact star-forming phase in the central regions of many growing galaxies that follows an earlier phase of gas compaction and is followed by a central quenching phase. We train a Convolutional Neural Network (CNN) with mock "observed" images of simulated galaxies at three phases of evolution: pre-BN, BN and post-BN, and demonstrate that the CNN successfully retrieves the three phases in other simulated galaxies. We show that BNs are identified by the CNN within a time window of $\sim0.15$ Hubble times. When the trained CNN is applied to observed galaxies from the CANDELS survey at $z=1-3$, it successfully identifies galaxies at the three phases. We find that the observed BNs are preferentially found in galaxies at a characteristic stellar mass range, $10^{9.2-10.3} M_\odot$ at all redshifts. This is consistent with the characteristic galaxy mass for BNs as detected in the simulations, and is meaningful because it is revealed in the observations when the direct information concerning the total galaxy luminosity has been eliminated from the training set. This technique can be applied to the classification of other astrophysical phenomena for improved comparison of theory and observations in the era of large imaging surveys and cosmological simulations.

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.

Understanding how bulges grow in galaxies is critical step towards unveiling the link between galaxy morphology and star-formation. To do so, it is necessary to decompose large sample of galaxies at different epochs into their main components (bulges and disks). This is particularly challenging, especially at high redshifts, where galaxies are poorly resolved. This work presents a catalog of bulge-disk decompositions of the surface brightness profiles of ~17.600 H-band selected galaxies in the CANDELS fields (F160W<23, 0<z<2) in 4 to 7 filters covering a spectral range of 430-1600nm. This is the largest available catalog of this kind up to z = 2. By using a novel approach based on deep-learning to select the best model to fit, we manage to control systematics arising from wrong model selection and obtain less contaminated samples than previous works. We show that the derived structural properties are within $\sim10-20\%$ of random uncertainties. We then fit stellar population models to the decomposed SEDs (Spectral Energy Distribution) of bulges and disks and derive stellar masses (and stellar mass bulge-to-total ratios) as well as rest-frame colors (U,V,J) for bulges and disks separately. All data products are publicly released with this paper and through the web page https://lerma.obspm.fr/huertas/form_CANDELS and will be used for scientific analysis in forthcoming works.

Using an isolated Milky Way-mass galaxy simulation, we compare results from 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the standardized package Grackle) and common analysis toolkit yt, all of which are publicly available. Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production are carefully constrained across code platforms. With numerical accuracy that resolves the disk scale height, we find that the codes overall agree well with one another in many dimensions including: gas and stellar surface densities, rotation curves, velocity dispersions, density and temperature distribution functions, disk vertical heights, stellar clumps, star formation rates, and Kennicutt-Schmidt relations. Quantities such as velocity dispersions are very robust (agreement within a few tens of percent at all radii) while measures like newly-formed stellar clump mass functions show more significant variation (difference by up to a factor of ~3). Systematic differences exist, for example, between mesh-based and particle-based codes in the low density region, and between more diffusive and less diffusive schemes in the high density tail of the density distribution. Yet intrinsic code differences are generally small compared to the variations in numerical implementations of the common subgrid physics such as supernova feedback. Our experiment reassures that, if adequately designed in accordance with our proposed common parameters, results of a modern high-resolution galaxy formation simulation are more sensitive to input physics than to intrinsic differences in numerical schemes.

SubHalo Abundance Matching (SHAM) assumes that one (sub)halo property, such as mass Mvir or peak circular velocity Vpeak, determines properties of the galaxy hosted in each (sub)halo such as its luminosity or stellar mass. This assumption implies that the dependence of Galaxy Luminosity Functions (GLFs) and the Galaxy Stellar Mass Function (GSMF) on environmental density is determined by the corresponding halo density dependence. In this paper, we test this by determining from an SDSS sample the observed dependence with environmental density of the ugriz GLFs and GSMF for all galaxies, and for central and satellite galaxies separately. We then show that the SHAM predictions are in remarkable agreement with these observations, even when the galaxy population is divided between central and satellite galaxies. However, we show that SHAM fails to reproduce the correct dependence between environmental density and g-r color for all galaxies and central galaxies, although it better reproduces the color dependence on environmental density of satellite galaxies.

Studying giant star-forming clumps in distant galaxies is important to understand galaxy formation and evolution. At present, however, observers and theorists have not reached a consensus on whether the observed "clumps" in distant galaxies are the same phenomenon that is seen in simulations. In this paper, as a step to establish a benchmark of direct comparisons between observations and theories, we publish a sample of clumps constructed to represent the commonly observed "clumps" in the literature. This sample contains 3193 clumps detected from 1270 galaxies at $0.5 \leq z < 3.0$. The clumps are detected from rest-frame UV images, as described in our previous paper. Their physical properties, e.g., rest-frame color, stellar mass (M*), star formation rate (SFR), age, and dust extinction, are measured by fitting the spectral energy distribution (SED) to synthetic stellar population models. We carefully test the procedures of measuring clump properties, especially the method of subtracting background fluxes from the diffuse component of galaxies. With our fiducial background subtraction, we find a radial clump U-V color variation, where clumps close to galactic centers are redder than those in outskirts. The slope of the color gradient (clump color as a function of their galactocentric distance scaled by the semi-major axis of galaxies) changes with redshift and M* of the host galaxies: at a fixed M*, the slope becomes steeper toward low redshift, and at a fixed redshift, it becomes slightly steeper with M*. Based on our SED-fitting, this observed color gradient can be explained by a combination of a negative age gradient, a negative E(B-V) gradient, and a positive specific star formation rate gradient of the clumps. We also find that the color gradients of clumps are steeper than those of intra-clump regions. [Abridged]

We present high spatial resolution (FWHM$\sim$0.14'') observations of the CO($8-7$) line in GDS-14876, a compact star-forming galaxy at $z=2.3$ with total stellar mass of $\log(M_{\star}/M_{\odot})=10.9$. The spatially resolved velocity map of the inner $r\lesssim1$~kpc reveals a continous velocity gradient consistent with the kinematics of a rotating disk with $v_{\rm rot}(r=1\rm kpc)=163\pm5$ km s$^{-1}$ and $v_{\rm rot}/\sigma\sim2.5$. The gas-to-stellar ratios estimated from CO($8-7$) and the dust continuum emission span a broad range, $f^{\rm CO}_{\rm gas}=M_{\rm gas}/M_{\star}=13-45\%$ and $f^{\rm cont}_{\rm gas}=50-67\%$, but are nonetheless consistent given the uncertainties in the conversion factors. The dynamical modeling yields a dynamical mass of$\log(M_{\rm dyn}/M_{\odot})=10.58^{+0.5}_{-0.2}$ which is lower, but still consistent with the baryonic mass, $\log$(M$_{\rm bar}$= M$_{\star}$ + M$^{\rm CO}_{\rm gas}$/M$_{\odot}$)$=11.0$, if the smallest CO-based gas fraction is assumed. Despite a low, overall gas fraction, the small physical extent of the dense, star-forming gas probed by CO($8-7$), $\sim3\times$ smaller than the stellar size, implies a strong concentration that increases the gas fraction up to $f^{\rm CO, 1\rm kpc}_{\rm gas}\sim 85\%$ in the central 1 kpc. Such a gas-rich center, coupled with a high star-formation rate, SFR$\sim$ 500 M$_{\odot}$ yr$^{-1}$, suggests that GDS-14876 is quickly assembling a dense stellar component (bulge) in a strong nuclear starburst. Assuming its gas reservoir is depleted without replenishment, GDS-14876 will quickly ($t_{\rm depl}\sim27$ Myr) become a compact quiescent galaxy that could retain some fraction of the observed rotational support.

We investigate the nature and evolution of large-scale structure within the SSA22 protocluster region at $z=3.09$ using cosmological simulations. A redshift histogram constructed from current spectroscopic observations of the SSA22 protocluster reveals two separate peaks at $z = 3.065$ (blue) and $z = 3.095$ (red). Based on these data, we report updated overdensity and mass calculations for the SSA22 protocluster. We find $\delta_{b,gal}=4.8 \pm 1.8$, $\delta_{r,gal}=9.5 \pm 2.0$ for the blue and red peaks, respectively, and $\delta_{t,gal}=7.6\pm 1.4$ for the entire region. These overdensities correspond to masses of $M_b = (0.76 \pm 0.17) \times 10^{15} h^{-1} M_{\odot}$, $M_r = (2.15 \pm 0.32) \times 10^{15} h^{-1} M_{\odot}$, and $M_t=(3.19 \pm 0.40) \times 10^{15} h^{-1} M_{\odot}$ for the red, blue, and total peaks, respectively. We use the Small MultiDark Planck (SMDPL) simulation to identify comparably massive $z\sim 3$ protoclusters, and uncover the underlying structure and ultimate fate of the SSA22 protocluster. For this analysis, we construct mock redshift histograms for each simulated $z\sim 3$ protocluster, quantitatively comparing them with the observed SSA22 data. We find that the observed double-peaked structure in the SSA22 redshift histogram corresponds not to a single coalescing cluster, but rather the proximity of a $\sim 10^{15}h^{-1} M_{\odot}$ protocluster and at least one $>10^{14} h^{-1} M_{\odot}$ cluster progenitor. Such associations in the SMDPL simulation are easily understood within the framework of hierarchical clustering of dark matter halos. We finally find that the opportunity to observe such a phenomenon is incredibly rare, with an occurrence rate of $7.4h^3 \mbox{ Gpc}^{-3}$.

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 examine the fraction of massive ($M_{*}>10^{10} M_{\odot}$), compact star-forming galaxies (cSFGs) that host an active galactic nucleus (AGN) at $z\sim2$. These cSFGs are likely the direct progenitors of the compact quiescent galaxies observed at this epoch, which are the first population of passive galaxies to appear in large numbers in the early Universe. We identify cSFGs that host an AGN using a combination of Hubble WFC3 imaging and Chandra X-ray observations in four fields: the Chandra Deep Fields, the Extended Groth Strip, and the UKIDSS Ultra Deep Survey field. We find that $39.2^{+3.9}_{-3.6}$\% (65/166) of cSFGs at $1.4<z<3.0$ host an X-ray detected AGN. This fraction is 3.2 times higher than the incidence of AGN in extended star-forming galaxies with similar masses at these redshifts. This difference is significant at the $6.2\sigma$ level. Our results are consistent with models in which cSFGs are formed through a dissipative contraction that triggers a compact starburst and concurrent growth of the central black hole. We also discuss our findings in the context of cosmological galaxy evolution simulations that require feedback energy to rapidly quench cSFGs. We show that the AGN fraction peaks precisely where energy injection is needed to reproduce the decline in the number density of cSFGs with redshift. Our results suggest that the first abundant population of massive, quenched galaxies emerged directly following a phase of elevated supermassive black hole growth and further hints at a possible connection between AGN and the rapid quenching of star formation in these galaxies.

We present a simple, efficient and robust approach to improve cosmological redshift measurements. The method is based on the presence of a reference sample for which a precise redshift number distribution (dN/dz) can be obtained for different pencil-beam-like sub-volumes within the original survey. For each sub-volume we then impose: (i) that the redshift number distribution of the uncertain redshift measurements matches the reference dN/dz corrected by their selection functions; and (ii) the rank order in redshift of the original ensemble of uncertain measurements is preserved. The latter step is motivated by the fact that random variables drawn from Gaussian probability density functions (PDFs) of different means and arbitrarily large standard deviations satisfy stochastic ordering. We then repeat this simple algorithm for multiple arbitrary pencil-beam-like overlapping sub-volumes; in this manner, each uncertain measurement has multiple (non-independent) "recovered" redshifts which can be used to estimate a new redshift PDF. We refer to this method as the Stochastic Order Redshift Technique (SORT). We have used a state-of-the-art N-body simulation to test the performance of SORT under simple assumptions and found that it can improve the quality of cosmological redshifts in an efficient and robust manner. Particularly, SORT redshifts are able to recover the distinctive features of the 'cosmic web' and can provide unbiased measurement of the two-point correlation function on scales > 4 Mpc/h. Given its simplicity, we envision that a method like SORT can be incorporated into more sophisticated algorithms aimed to exploit the full potential of large extragalactic photometric surveys.

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 present new determinations of the stellar-to-halo mass relation (SHMR) at $z=0-10$ that match the evolution of the galaxy stellar mass function, the SFR$-M_*$ relation,and the cosmic star formation rate. We utilize a compilation of 40 observational studies from the literature and correct them for potential biases. Using our robust determinations of halo mass assembly and the SHMR, we infer star formation histories, merger rates, and structural properties for average galaxies, combining star-forming and quenched galaxies. Our main findings: (1) The halo mass $M_{50}$ above which 50\% of galaxies are quenched coincides with sSFR/sMAR$\sim1$, where sMAR is the specific halo mass accretion rate. (2) $M_{50}$ increases with redshift, presumably due to cold streams being more efficient at high redshift while virial shocks and AGN feedback become more relevant at lower redshifts. (3) The ratio sSFR/sMAR has a peak value, which occurs around $M_{\rm vir}\sim2\times10^{11}M_{\odot}$. (4) The stellar mass density within 1 kpc, $\Sigma_1$, is a good indicator of the galactic global sSFR. (5) Galaxies are statistically quenched after they reach a maximum in $\Sigma_1$, consistent with theoretical expectations of the gas compaction model; this maximum depends on redshift. (6) In-situ star formation is responsible for most galactic stellar mass growth, especially for lower-mass galaxies. (7) Galaxies grow inside out. The marked change in the slope of the size--mass relation when galaxies became quenched, from $d\log R_{\rm eff}/d\log M_*\sim0.35$ to $\sim2.5$, could be the result of dry minor mergers.

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 explore the evolution of the internal gas kinematics of star-forming galaxies from the peak of cosmic star-formation at $z\sim2$ to today. Measurements of galaxy rotation velocity $V_{rot}$, which quantify ordered motions, and gas velocity dispersion $\sigma_g$, which quantify disordered motions, are adopted from the DEEP2 and SIGMA surveys. This sample covers a continuous baseline in redshift from $z=2.5$ to $z=0.1$, spanning 10 Gyrs. At low redshift, nearly all sufficiently massive star-forming galaxies are rotationally supported ($V_{rot}>\sigma_g$). By $z=2$, the percentage of galaxies with rotational support has declined to 50$\%$ at low stellar mass ($10^{9}-10^{10}\,M_{\odot}$) and 70$\%$ at high stellar mass ($10^{10}-10^{11}M_{\odot}$). For $V_{rot}\,>\,3\,\sigma_g$, the percentage drops below 35$\%$ for all masses. From $z\,=\,2$ to now, galaxies exhibit remarkably smooth kinematic evolution on average. All galaxies tend towards rotational support with time, and it is reached earlier in higher mass systems. This is mostly due to an average decline in $\sigma_g$ by a factor of 3 since a redshift of 2, which is independent of mass. Over the same time period, $V_{rot}$ increases by a factor of 1.5 for low mass systems, but does not evolve for high mass systems. These trends in $V_{rot}$ and $\sigma_g$ with time are at a fixed stellar mass and should not be interpreted as evolutionary tracks for galaxy populations. When galaxy populations are linked in time with abundance matching, not only does $\sigma_g$ decline with time as before, but $V_{rot}$ strongly increases with time for all galaxy masses. This enhances the evolution in $V_{rot}/\sigma_g$. These results indicate that $z\,=\,2$ is a period of disk assembly, during which the strong rotational support present in today's massive disk galaxies is only just beginning to emerge.

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.

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 results of a high-resolution zoom cosmological simulation of the evolution of a low-mass galaxy with a maximum velocity of V=100 km/s at z=0, using the initial conditions from the AGORA project (Kim et al. 2014). The final disc-dominated galaxy is consistent with local disc scaling relations, such as the stellar-to-halo mass relation and the baryonic Tully-Fisher. The galaxy evolves from a compact, dispersion-dominated galaxy into a rotation-dominated but dynamically hot disc in about 0.5 Gyr (from z=1.4 to z=1.2). The disc dynamically cools down for the following 7 Gyr, as the gas velocity dispersion decreases over time, in agreement with observations. The primary cause of this slow evolution of velocity dispersion in this low-mass galaxy is stellar feedback. It is related to the decline in gas fraction, and to the associated gravitational disk instability, as the disc slowly settles from a global Toomre Q>1 turbulent disc to a marginally unstable disc (Q=1).

We study the correlation of galaxy structural properties with their location relative to the SFR-M* correlation, also known as the star formation "main sequence" (SFMS), in the CANDELS and GAMA surveys and in a semi-analytic model (SAM) of galaxy formation. We first study the distribution of median Sersic index, effective radius, star formation rate (SFR) density and stellar mass density in the SFR-M* plane. We then define a redshift dependent main sequence and examine the medians of these quantities as a function of distance from this main sequence, both above (higher SFRs) and below (lower SFRs). Finally, we examine the distributions of distance from the main sequence in bins of these quantities. We find strong correlations between all of these galaxy structural properties and the distance from the SFMS, such that as we move from galaxies above the SFMS to those below it, we see a nearly monotonic trend towards higher median Sersic index, smaller radius, lower SFR density, and higher stellar density. In the semi-analytic model, bulge growth is driven by mergers and disk instabilities, and is accompanied by the growth of a supermassive black hole which can regulate or quench star formation via Active Galactic Nucleus (AGN) feedback. We find that our model qualitatively reproduces the trends described above, supporting a picture in which black holes and bulges co-evolve, and AGN feedback plays a critical role in moving galaxies off of the SFMS.

We study the evolution of giant clumps in high-z disc galaxies using AMR cosmological simulations at redshifts z=6-1. Our sample consists of 34 galaxies, of halo masses 10^{11}-10^{12}M_s at z=2, run with and without radiation pressure (RP) feedback from young stars. While RP has little effect on the sizes and global stability of discs, it reduces the amount of star-forming gas by a factor of ~2, leading to a decrease in stellar mass by a similar factor by z~2. Both samples undergo violent disc instability (VDI) and form giant clumps of masses 10^7-10^9M_s at a similar rate, though RP significantly reduces the number of long-lived clumps. When RP is (not) included, clumps with circular velocity <40(20)km/s, baryonic surface density <200(100)M_s/pc^2 and baryonic mass <10^{8.2}(10^{7.3})M_s are short-lived, disrupted in a few free-fall times. The more massive and dense clumps survive and migrate toward the disc centre over a few disc orbital times. In the RP simulations, the distribution of clump masses and star-formation rates (SFRs) normalized to their host disc is very similar at all redshifts. They exhibit a truncated power-law with a slope slightly shallower than -2. Short-lived clumps preferentially have young stellar ages, low masses, high gas fractions and specific SFRs (sSFR), and they tend to populate the outer disc. The sSFR of massive, long-lived clumps declines with age as they migrate towards the disc centre, producing gradients in mass, stellar age, gas fraction, sSFR and metallicity that distinguish them from short-lived clumps. Ex situ mergers make up ~37% of the mass in clumps and ~29% of the SFR. They are more massive and with older stellar ages than the in situ clumps, especially near the disc edge. Roughly half the galaxies at redshifts z=4-1 are clumpy over a wide range of stellar mass, with clumps accounting for ~3-30% of the SFR but ~0.1-3% of the stellar mass.

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.

We present spatially-resolved Atacama Large Millimeter/sub-millimeter Array (ALMA) 870 $\mu$m dust continuum maps of six massive, compact, dusty star-forming galaxies (SFGs) at $z\sim2.5$. These galaxies are selected for their small rest-frame optical sizes ($r_{\rm e, F160W}\sim1.6$ kpc) and high stellar-mass densities that suggest that they are direct progenitors of compact quiescent galaxies at $z\sim2$. The deep observations yield high far-infrared (FIR) luminosities of L$_{\rm IR}=10^{12.3-12.8}$ L$_{\odot}$ and star formation rates (SFRs) of SFR$=200-700$ M$_{\odot}$yr$^{-1}$, consistent with those of typical star-forming "main sequence" galaxies. The high-spatial resolution (FWHM$\sim$0.12"-0.18") ALMA and HST photometry are combined to construct deconvolved, mean radial profiles of their stellar mass and (UV+IR) SFR. We find that the dusty, nuclear IR-SFR overwhelmingly dominates the bolometric SFR up to $r\sim5$ kpc, by a factor of over 100$\times$ from the unobscured UV-SFR. Furthermore, the effective radius of the mean SFR profile ($r_{\rm e, SFR}\sim1$ kpc) is $\sim$30% smaller than that of the stellar mass profile. The implied structural evolution, if such nuclear starburst last for the estimated gas depletion time of $\Delta t=\pm100$ Myr, is a 4$\times$ increase of the stellar mass density within the central 1 kpc and a 1.6$\times$ decrease of the half-mass radius. This structural evolution fully supports dissipation-driven, formation scenarios in which strong nuclear starbursts transform larger, star-forming progenitors into compact quiescent galaxies.

We examine a sample of 1495 galaxies in the CANDELS fields to determine the evolution of two component galaxies, including bulges and discs, within massive galaxies at the epoch 1 < z < 3 when the Hubble sequence forms. We fit all of our galaxies' light profiles with a single S\'ersic fit, as well as with a combination of exponential and S\'ersic profiles. The latter is done in order to describe a galaxy with an inner and an outer component, or bulge and disc component. We develop and use three classification methods (visual, F-test and the RFF) to separate our sample into 1-component galaxies (disc/spheroids-like galaxies) and 2-component galaxies (galaxies formed by an 'inner part' or bulge and an 'outer part' or disc). We then compare the results from using these three different ways to classify our galaxies. We find that the fraction of galaxies selected as 2-component galaxies increases on average 50 per cent from the lowest mass bin to the most massive galaxies, and decreases with redshift by a factor of four from z = 1 to z = 3. We find that single S\'ersic 'disc-like' galaxies have the highest relative number densities at all redshifts, and that 2-component galaxies have the greatest increase and become at par with S\'ersic discs by z = 1. We also find that the systems we classify as 2-component galaxies have an increase in the sizes of their outer components, or 'discs', by about a factor of three from z = 3 to z = 1.5, while the inner components or 'bulges' stay roughly the same size. This suggests that these systems are growing from the inside out, whilst the bulges or protobulges are in place early in the history of these galaxies. This is also seen to a lesser degree in the growth of single 'disc-like' galaxies vs. 'spheroid-like' galaxies over the same epoch.

The velocity structure imprinted in the H{\alpha} emission line profiles contains valuable information about galactic outflows. Using a set of high-resolution zoom-in cosmological simulations of galaxies at z=2, we generate H{\alpha} emission line profiles, taking into account the temperature-dependent H{\alpha} emissivity, as well as dust extinction. The H{\alpha} line can be described as a sum of two gaussians, as typically done with observations. In general, its properties are in good agreement with those observed in local isolated galaxies with similar masses and star formation rates, assuming a spatially constant clumping factor of c=24. Blueshifted outflows are very common in the sample. They extend several kpc above the galaxy discs. They are also spread over the full extent of the discs. However, at small radii, the material with high velocities tends to remain confined within a thick disc, as part of galactic fountains or a turbulent medium, most probably due to the deeper gravitational potential at the galaxy centre.

We present the stellar mass ($M_{*}$)--gas-phase metallicity relation (MZR) and its scatter at intermediate redshifts ($0.5\leq z\leq0.7$) for 1381 field galaxies collected from deep spectroscopic surveys. The star formation rate (SFR) and color at a given $M_{*}$ of this magnitude-limited ($R\lesssim24$ AB) sample are representative of normal star-forming galaxies. For masses below $10^9 M_\odot$, our sample of 237 galaxies is $\sim$10 times larger than those in previous studies beyond the local universe. This huge gain in sample size enables superior constraints on the MZR and its scatter in the low-mass regime. We find a power-law MZR at $10^{8} M_\odot < M_{*} < 10^{11} M_\odot$: ${12+log(O/H) = (5.83\pm0.19) + (0.30\pm0.02)log(M_{*}/M_\odot)}$. Our MZR shows good agreement with others measured at similar redshifts in the literature in the intermediate and massive regimes, but is shallower than the extrapolation of the MZRs of others to masses below $10^{9} M_\odot$. The SFR dependence of the MZR in our sample is weaker than that found for local galaxies (known as the Fundamental Metallicity Relation). Compared to a variety of theoretical models, the slope of our MZR for low-mass galaxies agrees well with predictions incorporating supernova energy-driven winds. Being robust against currently uncertain metallicity calibrations, the scatter of the MZR serves as a powerful diagnostic of the stochastic history of gas accretion, gas recycling, and star formation of low-mass galaxies. Our major result is that the scatter of our MZR increases as $M_{*}$ decreases. Our result implies that either the scatter of the baryonic accretion rate or the scatter of the $M_{*}$--$M_{halo}$ relation increases as $M_{*}$ decreases. Moreover, our measures of scatter at $z=0.7$ appears consistent with that found for local galaxies.

We study the evolution of global shapes of galaxies using cosmological simulations. The shapes refer to the components of dark matter (DM), stars and gas at the stellar half-mass radius. Most galaxies undergo a characteristic compaction event into a blue nugget at $z\sim2-4$, which marks the transition from a DM-dominated central body to a self-gravitating baryonic core. We find that in the high-$z$, DM-dominated phase, the stellar and DM systems tend to be triaxial, preferentially prolate and mutually aligned. The elongation is supported by an anisotropic velocity dispersion that originates from the assembly of the galaxy along a dominant large-scale filament. We estimate that torques by the dominant halo are capable of inducing the elongation of the stellar system and its alignment with the halo. Then, in association with the transition to self-gravity, small-pericenter orbits puff up and the DM and stellar systems evolve into a more spherical and oblate configuration, aligned with the gas disc and associated with rotation. This transition typically occurs when the stellar mass is $\sim 10^9$ M$_\odot$ and the escape velocity in the core is $\sim 100$ km s$^{-1}$, indicating that supernova feedback may be effective in keeping the core DM-dominated and the system prolate. The early elongated phase itself may be responsible for the compaction event, and the transition to the oblate phase may be associated with the subsequent quenching in the core.

Using cosmological simulations, we address the interplay between structure and star formation in high-redshift galaxies via the evolution of surface density profiles. Our sample consists of 26 galaxies evolving in the redshift range $z=7-1$, spanning the stellar mass range $(0.2-6.4)\times 10^{10}M_\odot$ at $z=2$. We recover the main trends by stacking the profiles in accordance to their evolution phases. Following a wet compaction event that typically occurs when the stellar mass is $\sim10^{9.5}~M_{\odot}$ at $z\sim2-4$, the gas develops a cusp inside the effective radius, associated with a peak in star-formation rate (SFR). The SFR peak and the associated feedback, in the absence of further gas inflow to the centre, marks the onset of gas depletion from the central 1 kpc, leading to quenching of the central SFR. An extended, star-forming ring that forms by fresh gas during the central quenching process shows as a rising specific SFR (sSFR) profile, which is interpreted as inside-out quenching. Before quenching, the stellar density profile grows self-similarly, maintaining its log-log shape because the sSFR is similar at all radii. During the quenching process, the stellar density saturates to a constant value, especially in the inner 1 kpc. The stellar mass and SFR profiles deduced from observations show very similar shapes, consistent with the scenario of wet compaction leading to inside-out quenching and the subsequent saturation of a dense stellar core. We predict a cuspy gas profile during the blue nugget phase, and a gas-depleted core, sometimes surrounded by a ring, in the post-blue nugget phase.

We present estimates of intrinsic scatter in the Star Formation Rate (SFR) - Stellar Mass (M*) correlation in the redshift range 0.5 < z < 3.0 and in the mass range 10^7 < M* < 10^11 Msun. We utilize photometry in the Hubble Ultradeep Field (HUDF12), Ultraviolet Ultra Deep Field (UVUDF) campaigns and CANDELS/GOODS-S. We estimate SFR, M* from broadband Spectral Energy Distributions (SEDs) and the best available redshifts. The maximum depth of the HUDF photometry (F160W 29.9 AB, 5 sigma depth) probes the SFR-M* correlation down to M* ~ 10 ^7 Msun, a factor of 10-100X lower in M* than previous studies, and comparable to dwarf galaxies in the local universe. We find the slope of the SFR-M* relationship to be near unity at all redshifts and the normalization to decrease with cosmic time. We find a moderate increase in intrinsic scatter with cosmic time from 0.2 to 0.4 dex across the epoch of peak cosmic star formation. None of our redshift bins show a statistically significant increase in intrinsic scatter at low mass. However, it remains possible that intrinsic scatter increases at low mass on timescales shorter than ~ 100 Myr. Our results are consistent with a picture of gradual and self-similar assembly of galaxies across more than three orders of magnitude in stellar mass from as low as 10^7 Msun.

Gas inflow feeds galaxies with low metallicity gas from the cosmic web, sustaining star formation across the Hubble time. We make a connection between these inflows and metallicity inhomogeneities in star-forming galaxies, by using synthetic narrow-band images of the Halpha emission line from zoom-in AMR cosmological simulations of galaxies with stellar masses of $M \simeq 10^9 $Msun at redshifts z=2-7. In $\sim$50\% of the cases at redshifts lower than 4, the gas inflow gives rise to star-forming, Halpha-bright, off-centre clumps. Most of these clumps have gas metallicities, weighted by Halpha luminosity, lower than the metallicity in the surrounding interstellar medium by $\sim$0.3 dex, consistent with observations of chemical inhomogeneities at high and low redshifts. Due to metal mixing by shear and turbulence, these metallicity drops are dissolved in a few disc dynamical times. Therefore, they can be considered as evidence for rapid gas accretion coming from cosmological inflow of pristine gas.

Important but rare and subtle processes driving galaxy morphology and star-formation may be missed by traditional spiral, elliptical, irregular or S\'ersic bulge/disk classifications. To overcome this limitation, we use a principal component analysis of non-parametric morphological indicators (concentration, asymmetry, Gini coefficient, $M_{20}$, multi-mode, intensity and deviation) measured at rest-frame $B$-band (corresponding to HST/WFC3 F125W at 1.4 $< z <$ 2) to trace the natural distribution of massive ($>10^{10} M_{\odot}$) galaxy morphologies. Principal component analysis (PCA) quantifies the correlations between these morphological indicators and determines the relative importance of each. The first three principal components (PCs) capture $\sim$75 per cent of the variance inherent to our sample. We interpret the first principal component (PC) as bulge strength, the second PC as dominated by concentration and the third PC as dominated by asymmetry. Both PC1 and PC2 correlate with the visual appearance of a central bulge and predict galaxy quiescence. PC1 is a better predictor of quenching than stellar mass, as as good as other structural indicators (S\'ersic-n or compactness). We divide the PCA results into groups using an agglomerative hierarchical clustering method. Unlike S\'ersic, this classification scheme separates compact galaxies from larger, smooth proto-elliptical systems, and star-forming disk-dominated clumpy galaxies from star-forming bulge-dominated asymmetric galaxies. Distinguishing between these galaxy structural types in a quantitative manner is an important step towards understanding the connections between morphology, galaxy assembly and star-formation.

Using cosmological simulations, we address the properties of high-redshift star-forming galaxies (SFGs) across their main sequence (MS) in the plane of star-formation rate (SFR) versus stellar mass. We relate them to the evolution of galaxies through phases of gas compaction, depletion, possible replenishment, and eventual quenching. We find that the high-SFR galaxies in the upper envelope of the MS are compact, with high gas fractions and short depletion times ("blue nuggets"), while the lower-SFR galaxies in the lower envelope have lower central gas densities, lower gas fractions and longer depletion times, consistent with observed gradients across the MS. Stellar-structure gradients are negligible. The SFGs oscillate about the MS ridge on timescales $\sim0.4~t_{\mathrm{Hubble}}$ ($\sim1$ Gyr at $z\sim3$). The propagation upwards is due to gas compaction, triggered, e.g., by mergers, counter-rotating streams, and/or violent disc instabilities. The downturn at the upper envelope is due to central gas depletion by peak star formation and outflows while inflow from the shrunken gas disc is suppressed. An upturn at the lower envelope can occur once the extended disc has been replenished by fresh gas and a new compaction can be triggered, namely as long as the replenishment time is shorter than the depletion time. The mechanisms of gas compaction, depletion and replenishment confine the SFGs to the narrow ($\pm0.3$ dex) MS. Full quenching occurs in massive haloes ($M_{\mathrm{vir}}>10^{11.5}~M_\odot$) and/or at low redshifts ($z<3$), where the replenishment time is long compared to the depletion time, explaining the observed bending down of the MS at the massive end.

We utilize zoom-in cosmological simulations to study the nature of violent disc instability (VDI) in clumpy galaxies at high redshift, $z=1$--$5$. Our simulated galaxies are not in the ideal state assumed in Toomre instability, of linear fluctuations in an isolated, uniform, rotating disk. There, instability is characterised by a $Q$ parameter below unity, and lower when the disk is thick. Instead, the high-redshift discs are highly perturbed. Over long periods they consist of non-linear perturbations, compact massive clumps and extended structures, with new clumps forming in inter-clump regions. This is while the galaxy is subject to frequent external perturbances. We compute the local, two-component $Q$ parameter for gas and stars, smoothed on a $\sim1~{\rm kpc}$ scale to capture clumps of $10^{8-9}~{\rm M}_\odot$. The $Q<1$ regions are confined to collapsed clumps due to the high surface density there, while the inter-clump regions show $Q$ significantly higher than unity. Tracing the clumps back to their relatively smooth Lagrangian patches, we find that $Q$ prior to clump formation typically ranges from unity to a few. This is unlike the expectations from standard Toomre instability. We discuss possible mechanisms for high-$Q$ clump formation, e.g. rapid turbulence decay leading to small clumps that grow by mergers, non-axisymmetric instability, or clump formation induced by non-linear perturbations in the disk. Alternatively, the high-$Q$ non-linear VDI may be stimulated by the external perturbations such as mergers and counter-rotating streams. The high $Q$ may represent excessive compressive modes of turbulence, possibly induced by tidal interactions.

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.

The spatial fluctuations of the extragalactic background light trace the total emission from all stars and galaxies in the Universe. A multi-wavelength study can be used to measure the integrated emission from first galaxies during reionization when the Universe was about 500 million years old. Here we report arcminute-scale spatial fluctuations in one of the deepest sky surveys with the Hubble Space Telescope in five wavebands between 0.6 and 1.6 $\mu$m. We model-fit the angular power spectra of intensity fluctuation measurements to find the ultraviolet luminosity density of galaxies at $z$ > 8 to be $\log \rho_{\rm UV} = 27.4^{+0.2}_{-1.2}$ erg s$^{-1}$ Hz$^{-1}$ Mpc$^{-3}$ $(1\sigma)$. This level of integrated light emission allows for a significant surface density of fainter primeval galaxies that are below the point source detection level in current surveys.

We report the identification of elongated (triaxial or prolate) galaxies in cosmological simulations at $z\simeq2$. These are preferentially low-mass galaxies ($M_s \le 10^{9.5} \ M_\odot$), residing in dark-matter (DM) haloes with strongly elongated inner parts, a common feature of high-redshift DM haloes in the $\Lambda$CDM cosmology. Feedback slows formation of stars at the centres of these halos, so that a dominant and prolate DM distribution gives rise to galaxies elongated along the DM major axis. As galaxies grow in stellar mass, stars dominate the total mass within the galaxy half-mass radius, making stars and DM rounder and more oblate. A large population of elongated galaxies produces a very asymmetric distribution of projected axis ratios, as observed in high-z galaxy surveys. This indicates that the majority of the galaxies at high redshifts are not discs or spheroids but rather galaxies with elongated morphologies.

We study the evolution of the core (r<1 kpc) and effective (r<r_e) stellar-mass surface densities, in star-forming and quiescent galaxies. Since z=3, both populations occupy distinct, linear relations in log(Sigma_e) and log(Sigma_1) vs. log(M). These structural relations exhibit slopes and scatter that remain almost constant with time while their normalizations decline. For SFGs, the normalization declines by less than a factor of 2 from z=3, in both Sigma_e and Sigma_1. Such mild declines suggest that SFGs build dense cores by growing along these relations. We define this evolution as the structural main sequence (Sigma-MS). Quiescent galaxies follow different relations (Sigma^Q_e, Sigma^Q_1) off the Sigma-MS by having higher densities than SFGs of the same mass and redshift. The normalization of Sigma^Q_e declines by a factor of 10 since z=3, but only a factor of 2 in Sigma^Q_1. Thus, the common denominator for quiescent galaxies at all redshifts is the presence of a dense stellar core, and the formation of such cores in SFGs is the main requirement for quenching. Expressed in 2D as deviations off the SFR-MS and off Sigma^Q_1 at each redshift, the distribution of massive galaxies forms a universal, L-shaped sequence that relates two fundamental physical processes: compaction and quenching. Compaction is a process of substantial core-growth in SFGs relative to that in the Sigma-MS. This process increases the core-to-total mass and Sersic index, thereby, making compact SFGs. Quenching occurs once compact SFGs reach a maximum central density above Sigma^Q_1 > 9.5 M_sun/kpc^2. This threshold provides the most effective selection criterion to identify the star-forming progenitors of quiescent galaxies at all redshifts.

From mock Hubble Space Telescope images, we quantify non-parametric statistics of galaxy morphology, thereby predicting the emergence of relationships among stellar mass, star formation, and observed rest-frame optical structure at 1 < z < 3. We measure automated diagnostics of galaxy morphology in cosmological simulations of the formation of 22 central galaxies with 9.3 < log10 M_*/M_sun < 10.7. These high-spatial-resolution zoom-in calculations enable accurate modeling of the rest-frame UV and optical morphology. Even with small numbers of galaxies, we find that structural evolution is neither universal nor monotonic: galaxy interactions can trigger either bulge or disc formation, and optically bulge-dominated galaxies at this mass may not remain so forever. Simulated galaxies with M_* > 10^10 M_sun contain relatively more disc-dominated light profiles than those with lower mass, reflecting significant disc brightening in some haloes at 1 < z < 2. By this epoch, simulated galaxies with specific star formation rates below 10^-9.7 yr^-1 are more likely than normal star-formers to have a broader mix of structural types, especially at M_* > 10^10 M_sun. We analyze a cosmological major merger at z ~ 1.5 and find that the newly proposed MID morphology diagnostics trace later merger stages while G-M20 trace earlier ones. MID is sensitive also to clumpy star-forming discs. The observability time of typical MID-enhanced events in our simulation sample is less than 100 Myr. A larger sample of cosmological assembly histories may be required to calibrate such diagnostics in the face of their sensitivity to viewing angle, segmentation algorithm, and various phenomena such as clumpy star formation and minor mergers.

The TMT Detailed Science Case describes the transformational science that the Thirty Meter Telescope will enable. Planned to begin science operations in 2024, TMT will open up opportunities for revolutionary discoveries in essentially every field of astronomy, astrophysics and cosmology, seeing much fainter objects much more clearly than existing telescopes. Per this capability, TMT's science agenda fills all of space and time, from nearby comets and asteroids, to exoplanets, to the most distant galaxies, and all the way back to the very first sources of light in the Universe. More than 150 astronomers from within the TMT partnership and beyond offered input in compiling the new 2015 Detailed Science Case. The contributing astronomers represent the entire TMT partnership, including the California Institute of Technology (Caltech), the Indian Institute of Astrophysics (IIA), the National Astronomical Observatories of the Chinese Academy of Sciences (NAOC), the National Astronomical Observatory of Japan (NAOJ), the University of California, the Association of Canadian Universities for Research in Astronomy (ACURA) and US associate partner, the Association of Universities for Research in Astronomy (AURA).

We examine the spheroid growth and star formation quenching experienced by galaxies from z~3 to the present by studying the evolution with redshift of the quiescent and spheroid-dominated fractions of galaxies from the CANDELS and GAMA surveys. We compare the observed fractions with predictions from a semi-analytic model which includes prescriptions for bulge growth and AGN feedback due to mergers and disk instabilities. We facilitate direct morphological comparison by converting our model bulge-to-total stellar mass ratios to Sersic indices. We then subdivide our population into the four quadrants of the sSFR-Sersic index plane and study the buildup of each of these subpopulations. We find that the fraction of star forming disks declines steadily, while the fraction of quiescent spheroids builds up over cosmic time. The fractions of star forming spheroids and quiescent disks are both non-negligible, and stay nearly constant over the period we have studied, at about 10% and 15-20% respectively. Our model is qualitatively successful at reproducing the evolution of the two "main" populations (star forming disk-dominated galaxies and quiescent spheroid-dominated galaxies), and approximately reproduces the relative fractions of all four types, but predicts a stronger decline in star forming spheroids, and increase in quiescent disks, than seen in the observations. A model with an additional channel for bulge growth via disk instabilities agrees better overall with the observations than a model in which bulges may grow only through mergers. We study evolutionary tracks of some individual galaxies as they experience morphological transformation and quenching, and examine the importance of different physical drivers of this transformation (major and minor mergers and disk instabilities). We find that complex histories with multiple transformative events are the norm.

LCDM is remarkably successful in predicting the cosmic microwave background and large-scale structure, and LCDM parameters have been determined with only mild tensions between different types of observations. Hydrodynamical simulations starting from cosmological initial conditions are increasingly able to capture the complex interactions between dark matter and baryonic matter in galaxy formation. Simulations with relatively low resolution now succeed in describing the overall galaxy population. For example, the EAGLE simulation in volumes up to 100 cubic Mpc reproduces the observed local galaxy mass function nearly as well as semi-analytic models. It once seemed that galaxies are pretty smooth, that they generally grow in size as they evolve, and that they are a combination of disks and spheroids. But recent HST observations combined with high-resolution hydrodynamic simulations are showing that most star-forming galaxies are very clumpy; that galaxies often undergo compaction which reduces their radius and increases their central density; and that most lower-mass star-forming galaxies are not spheroids or disks but are instead elongated when their centers are dominated by dark matter. We also review LCDM challenges on smaller scales: cusp-core, "too big to fail," and substructure issues. Although starbursts can rapidly drive gas out of galaxy centers and thereby reduce the dark matter density, it remains to be seen whether this or other baryonic physics can explain the observed rotation curves of the entire population of dwarf and low surface brightness galaxies. If not, perhaps more complicated physics such as self-interacting dark matter may be needed. But standard LCDM appears to be successful in predicting the dark matter halo substructure that is now observed via gravitational lensing and breaks in cold stellar streams, and any alternative theory must do at least as well.

We present Keck-I MOSFIRE spectroscopy in the Y and H bands of GDN-8231, a massive, compact, star-forming galaxy (SFG) at a redshift $z\sim1.7$. Its spectrum reveals both H$_{\alpha}$ and [NII] emission lines and strong Balmer absorption lines. The H$_{\alpha}$ and Spitzer MIPS 24 $\mu$m fluxes are both weak, thus indicating a low star formation rate of SFR $\lesssim5-10$ M$_{\odot}$ yr$^{-1}$. This, added to a relatively young age of $\sim700$ Myr measured from the absorption lines, provides the first direct evidence for a distant galaxy being caught in the act of rapidly shutting down its star formation. Such quenching allows GDN-8231 to become a compact, quiescent galaxy, similar to 3 other galaxies in our sample, by $z\sim1.5$. Moreover, the color profile of GDN-8231 shows a bluer center, consistent with the predictions of recent simulations for an early phase of inside-out quenching. Its line-of-sight velocity dispersion for the gas, $\sigma^{\rm{gas}}_{\!_{\rm LOS}}=127\pm32$ km s$^{-1}$, is nearly 40% smaller than that of its stars, $\sigma^{\star}_{\!_{\rm LOS}}=215\pm35$ km s$^{-1}$. High-resolution hydro-simulations of galaxies explain such apparently colder gas kinematics of up to a factor of $\sim1.5$ with rotating disks being viewed at different inclinations and/or centrally concentrated star-forming regions. A clear prediction is that their compact, quiescent descendants preserve some remnant rotation from their star-forming progenitors.

We use cosmological simulations to study a characteristic evolution pattern of high redshift galaxies. Early, stream-fed, highly perturbed, gas-rich discs undergo phases of dissipative contraction into compact, star-forming systems (blue nuggets) at z~4-2. The peak of gas compaction marks the onset of central gas depletion and inside-out quenching into compact ellipticals (red nuggets) by z~2. These are sometimes surrounded by gas rings or grow extended dry stellar envelopes. The compaction occurs at a roughly constant specific star-formation rate (SFR), and the quenching occurs at a constant stellar surface density within the inner kpc ($\Sigma_1$). Massive galaxies quench earlier, faster, and at a higher $\Sigma_1$ than lower-mass galaxies, which compactify and attempt to quench more than once. This evolution pattern is consistent with the way galaxies populate the SFR-radius-mass space, and with gradients and scatter across the main sequence. The compaction is triggered by an intense inflow episode, involving (mostly minor) mergers, counter-rotating streams or recycled gas, and is commonly associated with violent disc instability. The contraction is dissipative, with the inflow rate >SFR, and the maximum $\Sigma_1$ anti-correlated with the initial spin parameter, as predicted by Dekel & Burkert (2014). The central quenching is triggered by the high SFR and stellar/supernova feedback (possibly also AGN feedback) due to the high central gas density, while the central inflow weakens as the disc vanishes. Suppression of fresh gas supply by a hot halo allows the long-term maintenance of quenching once above a threshold halo mass, inducing the quenching downsizing.

Although giant clumps of stars are crucial to galaxy formation and evolution, the most basic demographics of clumps are still uncertain, mainly because the definition of clumps has not been thoroughly discussed. In this paper, we study the basic demographics of clumps in star-forming galaxies (SFGs) at 0.5<z<3, using our proposed physical definition that UV-bright clumps are discrete star-forming regions that individually contribute more than 8% of the rest-frame UV light of their galaxies. Clumps defined this way are significantly brighter than the HII regions of nearby large spiral galaxies, either individually or blended, when physical spatial resolution and cosmological dimming are considered. Under this definition, we measure the fraction of SFGs that contain at least one off-center clump (Fclumpy) and the contributions of clumps to the rest-frame UV light and star formation rate of SFGs in the CANDELS/GOODS-S and UDS fields, where our mass-complete sample consists of 3239 galaxies with axial ratio q>0.5. The redshift evolution of Fclumpy changes with the stellar mass (M*) of the galaxies. Low-mass (log(M*/Msun)<9.8) galaxies keep an almost constant Fclumpy of about 60% from z~3.0 to z~0.5. Intermediate-mass and massive galaxies drop their Fclumpy from 55% at z~3.0 to 40% and 15%, respectively, at z~0.5. We find that (1) the trend of disk stabilization predicted by violent disk instability matches the Fclumpy trend of massive galaxies; (2) minor mergers are a viable explanation of the Fclumpy trend of intermediate-mass galaxies at z<1.5, given a realistic observability timescale; and (3) major mergers are unlikely responsible for the Fclumpy trend in all masses at z<1.5. The clump contribution to the rest-frame UV light of SFGs shows a broad peak around galaxies with log(M*/Msun)~10.5 at all redshifts, possibly linked to the molecular gas fraction of the galaxies. (Abridged)

Despite recent success in forming realistic present-day galaxies, simulations still form the bulk of their stars earlier than observations indicate. We investigate the process of stellar mass assembly in low-mass field galaxies, a dwarf and a typical spiral, focusing on the effects of radiation from young stellar clusters on the star formation (SF) histories. We implement a novel model of SF with a deterministic low efficiency per free-fall time, as observed in molecular clouds. Stellar feedback is based on observations of star-forming regions, and includes radiation pressure from massive stars, photoheating in H II regions, supernovae and stellar winds. We find that stellar radiation has a strong effect on the formation of low-mass galaxies, especially at z > 1, where it efficiently suppresses SF by dispersing cold and dense gas, preventing runaway growth of the stellar component. This behaviour is evident in a variety of observations but had so far eluded analytical and numerical models without radiation feedback. Compared to supernovae alone, radiation feedback reduces the SF rate by a factor of ~100 at z < 2, yielding rising SF histories which reproduce recent observations of Local Group dwarfs. Stellar radiation also produces bulgeless spiral galaxies and may be responsible for excess thickening of the stellar disc. The galaxies also feature rotation curves and baryon fractions in excellent agreement with current data. Lastly, the dwarf galaxy shows a very slow reduction of the central dark matter density caused by radiation feedback over the last ~7 Gyr of cosmic evolution.

We study the angular-momentum (AM) buildup in high-$z$ massive galaxies using high-resolution cosmological simulations. The AM originates in co-planar streams of cold gas and merging galaxies tracing cosmic-web filaments, and it undergoes four phases of evolution. (I) Outside the halo virial radius ($R_{\rm v}\!\sim\!100\,{\rm kpc}$), the elongated streams gain AM by tidal torques with a specific AM (sAM) $\sim\!1.7$ times the dark-matter (DM) spin due to the gas' higher quadrupole moment. This AM is expressed as stream impact parameters, from $\sim\!0.3R_{\rm v}$ to counter rotation. (II) In the outer halo, while the incoming DM mixes with the existing halo of lower sAM to a spin $\lambda_{\rm dm}\!\sim\!0.04$, the cold streams transport the AM to the inner halo such that their spin in the halo is $\sim\!3\lambda_{\rm dm}$. (III) Near pericenter, the streams dissipate into an irregular rotating ring extending to $\sim\!0.3R_{\rm v}$ and tilted relative to the inner disc. Torques exerted partly by the disc make the ring gas lose AM, spiral in, and settle into the disc within one orbit. The ring is observable with 30\% probability as a damped Lyman-$\alpha$ absorber. (IV) Within the disc, $<\!0.1R_{\rm v}$, torques associated with violent disc instability drive AM out and baryons into a central bulge, while outflows remove low-spin gas, introducing certain sensitivity to feedback strength. Despite the different AM histories of gas and DM, the disc spin is comparable to the DM-halo spin. Counter rotation can strongly affect disc evolution.

We address the formation of massive stellar spheroids between redshifts $z=4$ and 1 using a suite of AMR hydro-cosmological simulations. The spheroids form as bulges, and the spheroid mass growth is partly driven by violent disc instability (VDI) and partly by mergers. A kinematic decomposition to disc and spheroid yields that the mass fraction in the spheroid is between 50\% and 90\% and is roughly constant in time, consistent with a cosmological steady state of VDI discs that are continuously fed from the cosmic web. The density profile of the spheroid is typically "classical", with a Sersic index $n = 4.5\pm 1$, independent of whether it grew by mergers or VDI and independent of the feedback strength. The disc is characterized by $n=1.5\pm 0.5$, and the whole galaxy by $n=3\pm 1$. The high-redshift spheroids are compact due to the dissipative inflow of gas and the high universal density. The stellar surface density within the effective radius of each galaxy as it evolves remains roughly constant in time after its first growth. For galaxies of a fixed stellar mass, the surface density is higher at higher redshifts.

We estimate the galaxy density along lines of sight to hard extragalactic gamma-ray sources by correlating source positions on the sky with a void catalog based on the Sloan Digital Sky Survey (SDSS). Extragalactic gamma-ray sources that are detected at very high energy (VHE; E>100 GeV) or have been highlighted as VHE-emitting candidates in the Fermi Large Area Telescope hard source catalog (together referred to as "VHE-like" sources) are distributed along underdense lines of sight at the 2.4 sigma level. There is also a less suggestive correlation for the Fermi hard source population (1.7 sigma). A correlation between 10-500 GeV flux and underdense fraction along the line of sight for VHE-like and Fermi hard sources is found at 2.4 sigma and 2.6 sigma, respectively. The preference for underdense sight lines is not displayed by gamma-ray emitting galaxies within the second Fermi catalog, containing sources detected above 100 MeV, or the SDSS DR7 quasar catalog. We investigate whether this marginal correlation might be a result of lower extragalactic background light (EBL) photon density within the underdense regions and find that, even in the most extreme case of a entirely underdense sight line, the EBL photon density is only 2% less than the nominal EBL density. Translating this into gamma-ray attenuation along the line of sight for a highly attenuated source with opacity tau(E,z) ~5, we estimate that the attentuation of gamma-rays decreases no more than 10%. This decrease, although non-neglible, is unable to account for the apparent hard source correlation with underdense lines of sight.

We determine the intrinsic, 3-dimensional shape distribution of star-forming galaxies at 0<z<2.5, as inferred from their observed projected axis ratios. In the present-day universe star-forming galaxies of all masses 1e9 - 1e11 Msol are predominantly thin, nearly oblate disks, in line with previous studies. We now extend this to higher redshifts, and find that among massive galaxies (M* > 1e10 Msol) disks are the most common geometric shape at all z < 2. Lower-mass galaxies at z>1 possess a broad range of geometric shapes: the fraction of elongated (prolate) galaxies increases toward higher redshifts and lower masses. Galaxies with stellar mass 1e9 Msol (1e10 Msol) are a mix of roughly equal numbers of elongated and disk galaxies at z~1 (z~2). This suggests that galaxies in this mass range do not yet have disks that are sustained over many orbital periods, implying that galaxies with present-day stellar mass comparable to that of the Milky Way typically first formed such sustained stellar disks at redshift z~1.5-2. Combined with constraints on the evolution of the star formation rate density and the distribution of star formation over galaxies with different masses, our findings imply that, averaged over cosmic time, the majority of stars formed in disks.

Two-dimensional integral field surveys such as ATLAS^3D are producing rich observational data sets yielding insights into galaxy formation. These new kinematic observations have highlighted the need to understand the evolutionary mechanisms leading to a spectrum of fast-rotators and slow-rotators in early-type galaxies. We address the formation of slow and fast rotators through a series of controlled, comprehensive hydrodynamical simulations sampling idealized galaxy merger scenarios constructed from model spiral galaxies. Idealized and controlled simulations of this sort complement the more 'realistic' cosmological simulations by isolating and analyzing the effects of specific parameters, as we do in this paper. We recreate minor and major binary mergers, binary merger trees with multiple progenitors, and multiple sequential mergers. Within each of these categories of formation history, we correlate progenitor gas fraction, mass ratio, orbital pericenter, orbital ellipticity, and spin with remnant kinematic properties. We create kinematic profiles of these 95 simulations comparable to ATLAS^3D data. By constructing remnant profiles of the projected specific angular momentum (lambda_R = <R|V|> / <sqrt(V^2+sigma^2)>, triaxiality, and measuring the incidences of kinematic twists and kinematically decoupled cores, we distinguish between varying formation scenarios. We find that binary mergers nearly always form fast rotators. Slow rotators can be formed from zero initial angular momentum configurations and gas-poor mergers, but are not as round as the ATLAS^3D galaxies. Remnants of binary merger trees are triaxial slow rotators. Sequential mergers form round slow rotators that most resemble the ATLAS^3D rotators.

Recent observations have probed the formation histories of nearby elliptical galaxies by tracking correlations between the stellar population parameters, age and metallicity, and the structural parameters that enter the Fundamental Plane, radius and velocity dispersion sigma. These studies have found intriguing correlations between these four parameters. In this work, we make use of a semi-analytic model, based on halo merger trees extracted from the Bolshoi cosmological simulation, that predicts the structural properties of spheroid-dominated galaxies based on an analytic model that has been tested and calibrated against an extensive suite of hydrodynamic+N-body binary merger simulations. We predict the radius, sigma, luminosity, age, and metallicity of spheroid-dominated galaxies, enabling us to compare directly to observations. Our model predicts a strong correlation between age and sigma for early-type galaxies, and no significant correlation between age and radius, in agreement with observations. In addition we predict a strong correlation between metallicity and sigma, and a weak correlation between metallicity and radius, in qualitative agreement with observations. We find that the correlations with sigma arise as a result of the strong link between sigma and the galaxy assembly time. Minor mergers produce a large change in radius while leaving sigma nearly the same, which explains the weaker trends with radius.

We study the properties of giant clumps and their radial gradients in high-$z$ disc galaxies using AMR cosmological simulations. Our sample consists of 770 snapshots in the redshift range $z=4-1$ from 29 galaxies that at $z=2$ span the stellar mass range $(0.2-3)\times 10^{11}M_{\odot}$. Extended gas discs exist in 83% of the snapshots. Clumps are identified by gas density in 3D and their stellar and dark matter components are considered thereafter. While most of the overdensities are diffuse and elongated, 91% of their mass and 83% of their star-fromation rate (SFR) are in compact round clumps. Nearly all galaxies have a central, massive bulge clump, while 70% of the discs show off-center clumps, 3-4 per galaxy. The fraction of clumpy discs peaks at intermediate disc masses. Clumps are divided based on dark-matter content into $\textit{in-situ}$ and $\textit{ex-situ}$, originating from violent disc instability (VDI) and minor mergers respectively. 60% of the discs are in a VDI phase showing off-center $\textit{in-situ}$ clumps, which contribute 1-7% of the disc mass and 5-45% of its SFR. The $\textit{in-situ}$ clumps constitute 75% of the off-center clumps in terms of number and SFR but only half the mass, each clump containing on average 1% of the disc mass and 6% of its SFR. They have young stellar ages, $100-400 {\rm Myr}$, and high specific SFR (sSFR), $1-10 {\rm Gyr}^{-1}$. They exhibit gradients resulting from inward clump migration, where the inner clumps are somewhat more massive and older, with lower gas fraction and sSFR and higher metallicity. Similar observed gradients indicate that clumps survive outflows. The \exsitu clumps have stellar ages $0.5-3 {\rm Gyr}$ and sSFR $\sim 0.1-2 {\rm Gyr}^{-1}$, and they exhibit weaker gradients. Massive clumps of old stars at large radii are likely \exsitu mergers, though half of them share the disc rotation.

We use a large suite of hydrodynamical simulations of binary galaxy mergers to construct and calibrate a physical prescription for computing the effective radii and velocity dispersions of spheroids. We implement this prescription within a semi-analytic model embedded in merger trees extracted from the Bolshoi Lambda-CDM N-body simulation, accounting for spheroid growth via major and minor mergers as well as disk instabilities. We find that without disk instabilities, our model does not predict sufficient numbers of intermediate mass early-type galaxies in the local universe. Spheroids also form earlier in models with spheroid growth via disk instabilities. Our model correctly predicts the normalization, slope, and scatter of the low-redshift size-mass and Fundamental Plane relations for early type galaxies. It predicts a degree of curvature in the Faber-Jackson relation that is not seen in local observations, but this could be alleviated if higher mass spheroids have more bottom-heavy initial mass functions. The model also correctly predicts the observed strong evolution of the size-mass relation for spheroids out to higher redshifts, as well as the slower evolution in the normalization of the Faber-Jackson relation. We emphasize that these are genuine predictions of the model since it was tuned to match hydrodynamical simulations and not these observations.

Any successful model of galaxy formation needs to explain the low rate of star formation in the small progenitors of today's galaxies. This inefficiency is necessary for reproducing the low stellar-to-virial mass fractions, suggested by current abundance matching models. A possible driver of this low efficiency is the radiation pressure exerted by ionizing photons from massive stars. The effect of radiation pressure in cosmological, zoom-in galaxy formation simulations is modeled as a non-thermal pressure that acts only in dense and optically thick star-forming regions. We also include photoionization and photoheating by massive stars. The full photoionization of hydrogen reduces the radiative cooling in the $10^{4-4.5}$ K regime. The main effect of radiation pressure is to regulate and limit the high values of gas density and the amount of gas available for star formation. This maintains a low star formation rate of $\sim 1 \ {\rm M_\odot} \ {\rm yr}^{-1}$ in halos with masses about $10^{11} \ {M_\odot}$ at $z\simeq3$. Infrared trapping and photoionization/photoheating processes are secondary effects in this mass range. The galaxies residing in these low-mass halos contain only $\sim0.6\%$ of the total virial mass in stars, roughly consistent with abundance matching. Radiative feedback maintains an extended galaxy with a rising circular velocity profile.

We present Keck-I MOSFIRE near-infrared spectroscopy for a sample of 13 compact star-forming galaxies (SFGs) at redshift $2\leq z \leq2.5$ with star formation rates of SFR$\sim$100M$_{\odot}$ y$^{-1}$ and masses of log(M/M$_{\odot}$)$\sim10.8$. Their high integrated gas velocity dispersions of $\sigma_{\rm{int}}$=230$^{+40}_{-30}$ km s$^{-1}$, as measured from emission lines of H$_{\alpha}$ and [OIII], and the resultant M$_{\star}-\sigma_{\rm{int}}$ relation and M$_{\star}$$-$M$_{\rm{dyn}}$ all match well to those of compact quiescent galaxies at $z\sim2$, as measured from stellar absorption lines. Since log(M$_{\star}$/M$_{\rm{dyn}}$)$=-0.06\pm0.2$ dex, these compact SFGs appear to be dynamically relaxed and more evolved, i.e., more depleted in gas and dark matter ($<$13$^{+17}_{-13}$\%) than their non-compact SFG counterparts at the same epoch. Without infusion of external gas, depletion timescales are short, less than $\sim$300 Myr. This discovery adds another link to our new dynamical chain of evidence that compact SFGs at $z\gtrsim2$ are already losing gas to become the immediate progenitors of compact quiescent galaxies by $z\sim2$.

Cosmological simulations of galaxies have typically produced too many stars at early times. We study the global and morphological effects of radiation pressure (RP) in eight pairs of high-resolution cosmological galaxy formation simulations. We find that the additional feedback suppresses star formation globally by a factor of ~2. Despite this reduction, the simulations still overproduce stars by a factor of ~2 with respect to the predictions provided by abundance matching methods for halos more massive than 5E11 Msun/h (Behroozi, Wechsler & Conroy 2013). We also study the morphological impact of radiation pressure on our simulations. In simulations with RP the average number of low mass clumps falls dramatically. Only clumps with stellar masses Mclump/Mdisk <= 5% are impacted by the inclusion of RP, and RP and no-RP clump counts above this range are comparable. The inclusion of RP depresses the contrast ratios of clumps by factors of a few for clump masses less than 5% of the disk masses. For more massive clumps, the differences between and RP and no-RP simulations diminish. We note however, that the simulations analyzed have disk stellar masses below about 2E10 Msun/h. By creating mock Hubble Space Telescope observations we find that the number of clumps is slightly reduced in simulations with RP. However, since massive clumps survive the inclusion of RP and are found in our mock observations, we do not find a disagreement between simulations of our clumpy galaxies and observations of clumpy galaxies. We demonstrate that clumps found in any single gas, stellar, or mock observation image are not necessarily clumps found in another map, and that there are few clumps common to multiple maps.

Exploiting the deep high-resolution imaging of all 5 CANDELS fields, and accurate redshift information provided by 3D-HST, we investigate the relation between structure and stellar populations for a mass-selected sample of 6764 galaxies above 10^10 Msun, spanning the redshift range 0.5 < z < 2.5. For the first time, we fit 2-dimensional models comprising a single Sersic fit and two-component (i.e., bulge + disk) decompositions not only to the H-band light distributions, but also to the stellar mass maps reconstructed from resolved stellar population modeling. We confirm that the increased bulge prominence among quiescent galaxies, as reported previously based on rest-optical observations, remains in place when considering the distributions of stellar mass. Moreover, we observe an increase of the typical Sersic index and bulge-to-total ratio (with median B/T reaching 40-50%) among star-forming galaxies above 10^11 Msun. Given that quenching for these most massive systems is likely to be imminent, our findings suggest that significant bulge growth precedes a departure from the star-forming main sequence. We demonstrate that the bulge mass (and ideally knowledge of the bulge and total mass) is a more reliable predictor of the star-forming versus quiescent state of a galaxy than the total stellar mass. The same trends are predicted by the state-of-the-art semi-analytic model by Somerville et al. In the latter, bulges and black holes grow hand in hand through merging and/or disk instabilities, and AGN-feedback shuts off star formation. Further observations will be required to pin down star formation quenching mechanisms, but our results imply they must be internal to the galaxies and closely associated with bulge growth.

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

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

We analyze the star-forming and structural properties of 45 massive (log(M/Msun)>10) compact star-forming galaxies (SFGs) at 2<z<3 to explore whether they are progenitors of compact quiescent galaxies at z~2. The optical/NIR and far-IR Spitzer/Herschel colors indicate that most compact SFGs are heavily obscured. Nearly half (47%) host an X-ray bright AGN. In contrast, only about 10% of other massive galaxies at that time host AGNs. Compact SFGs have centrally-concentrated light profiles and spheroidal morphologies similar to quiescent galaxies, and are thus strikingly different from other SFGs. Most compact SFGs lie either within the SFR-M main sequence (65%) or below (30%), on the expected evolutionary path towards quiescent galaxies. These results show conclusively that galaxies become more compact before they lose their gas and dust, quenching star formation. Using extensive HST photometry from CANDELS and grism spectroscopy from the 3D-HST survey, we model their stellar populations with either exponentially declining (tau) star formation histories (SFHs) or physically-motivated SFHs drawn from semi-analytic models (SAMs). SAMs predict longer formation timescales and older ages ~2 Gyr, which are nearly twice as old as the estimates of the tau models. While both models yield good SED fits, SAM SFHs better match the observed slope and zero point of the SFR-M main sequence. Some low-mass compact SFGs (log(M/Msun)=10-10.6) have younger ages but lower sSFRs than that of more massive galaxies, suggesting that the low-mass galaxies reach the red sequence faster. If the progenitors of compact SFGs are extended SFGs, state-of-the-art SAMs show that mergers and disk instabilities are both able to shrink galaxies, but disk instabilities are more frequent (60% versus 40%) and form more concentrated galaxies. We confirm this result via high-resolution hydrodynamic simulations.

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

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.

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

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

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

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 evidence from a sample of 544 galaxies from the DEEP2 Survey for evolution of the internal kinematics of blue galaxies with stellar masses ranging 8.0 < log M* (M_Sun) < 10.7 over 0.2<z<1.2. DEEP2 provides galaxy spectra and Hubble imaging from which we measure emission-line kinematics and galaxy inclinations, respectively. Our large sample allows us to overcome scatter intrinsic to galaxy properties in order to examine trends in kinematics. We find that at a fixed stellar mass galaxies systematically decrease in disordered motions and increase in rotation velocity and potential well depth with time. Massive galaxies are the most well-ordered at all times examined, with higher rotation velocities and less disordered motions than less massive galaxies. We quantify disordered motions with an integrated gas velocity dispersion corrected for beam smearing (sigma_g). It is unlike the typical pressure-supported velocity dispersion measured for early type galaxies and galaxy bulges. Because both seeing and the width of our spectral slits comprise a significant fraction of the galaxy sizes, sigma_g integrates over velocity gradients on large scales which can correspond to non-ordered gas kinematics. We compile measurements of galaxy kinematics from the literature over 1.2<z<3.8 and do not find any trends with redshift, likely for the most part because these datasets are biased toward the most highly star-forming systems. In summary, over the last ~8 billion years since z=1.2, blue galaxies evolve from disordered to ordered systems as they settle to become the rotation-dominated disk galaxies observed in the Universe today, with the most massive galaxies being the most evolved at any time.

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

The definition of homology for single-component galaxies is clear, but for multi-component (luminous and dark matter) galaxies there is some ambiguity. We attempt to clarify the situation by carefully separating the different concepts of homology that have been used to date. We argue that the most useful definition is that a set of galaxies is homologous if they are the same in all respects up to a set of three dimensional scaling constants which may differ from one galaxy to the next. Noting that we are free to choose the dimensional constants, we find that a set of hydrodynamic simulated galaxy merger remnants is significantly closer to homologous when the dimensional length constant is taken to be the radius containing equal amounts of dark and baryonic matter rather than the usual observationally motivated choice of the baryonic half-mass radius. Once the correct dimensional scaling constants are used, the stellar velocity dispersion anisotropy is essentially the sole source of the variation in the kinematic structure of these simulated merger remnants. In order to facilitate the use of these scaling constants to analyse observed galaxies, we calculated the relationship between our preferred dimensional scaling constants and the typical observationally accessible quantities.

In the standard picture of disc galaxy formation, baryons and dark matter receive the same tidal torques, and therefore approximately the same initial specific angular momentum. However, observations indicate that disc galaxies typically have only about half as much specific angular momentum as their dark matter haloes. We argue this does not necessarily imply that baryons lose this much specific angular momentum as they form galaxies. It may instead indicate that galaxies are most directly related to the inner regions of their host haloes, as may be expected in a scenario where baryons in the inner parts of haloes collapse first. A limiting case is examined under the idealised assumption of perfect angular momentum conservation. Namely, we determine the density contrast Delta, with respect to the critical density of the Universe, by which dark matter haloes need to be defined in order to have the same average specific angular momentum as the galaxies they host. Under the assumption that galaxies are related to haloes via their characteristic rotation velocities, the necessary Delta is ~600. This Delta corresponds to an average halo radius and mass which are ~60% and ~75%, respectively, of the virial values (i.e., for Delta = 200). We refer to this radius as the radius of baryonic collapse R_BC, since if specific angular momentum is conserved perfectly, baryons would come from within it. It is not likely a simple step function due to the complex gastrophysics involved, therefore we regard it as an effective radius. In summary, the difference between the predicted initial and the observed final specific angular momentum of galaxies, which is conventionally attributed solely to angular momentum loss, can more naturally be explained by a preference for collapse of baryons within R_BC, with possibly some later angular momentum transfer.

Attenuation of high-energy gamma rays by pair-production with UV, optical and IR extragalactic background light (EBL) photons provides a link between the history of galaxy formation and high-energy astrophysics. We present results from our latest semi-analytic models (SAMs), which employ the main ingredients thought to be important to galaxy formation and evolution, as well as an improved model for reprocessing of starlight by dust to mid- and far-IR wavelengths. These SAMs are based upon a Lambda-CDM hierarchical structural formation scenario, and are successful in reproducing a large variety of observational constraints such as number counts, luminosity and mass functions, and color bimodality. Our fiducial model is based upon a WMAP5 cosmology, and treats dust emission using empirical templates. This model predicts a background flux considerably lower than optical and near-IR measurements that rely on subtraction of zodiacal and galactic foregrounds, and near the lower bounds set by number counts of resolvable sources at a large number of wavelengths. We also show the results of varying cosmological parameters and dust attenuation model used in our SAM. For each EBL prediction, we show how the optical depth due to electron-positron pair-production is affected by redshift and gamma-ray energy, and the effect of gamma-ray absorption on the spectra of a variety of extragalactic sources. We conclude with a discussion of the implications of our work, comparisons to other models and key measurements of the EBL and a discussion of how the burgeoning science of gamma-ray astronomy will continue to help constrain cosmology. The low EBL flux predicted by our fiducial model suggests an optimistic future for further studies of distant gamma-ray sources.

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

We combine a semi-analytic model of galaxy formation with simple analytic recipes describing the absorption and re-emission of starlight by dust in the interstellar medium of galaxies. We use the resulting models to predict galaxy counts and luminosity functions from the far-ultraviolet to the sub-mm, from redshift five to the present, and compare with an extensive compilation of observations. We find that in order to reproduce the rest-UV and optical luminosity functions at high redshift, we must assume an evolving normalization in the dust-to-metal ratio, implying that galaxies of a given bolometric luminosity (or metal column density) must be less extinguished than their local counterparts. In our best-fit model, we find remarkably good agreement with observations from rest-frame 1500 Angstroms to 250 microns. At longer wavelengths, most dramatically in the sub-mm, our models underpredict the number of bright galaxies by a large factor. The models reproduce the observed total IR luminosity function fairly well. We show the results of varying several ingredients of the models, including various aspects of the dust attenuation recipe, the dust emission templates, and the cosmology. We use our models to predict the integrated Extragalactic Background Light (EBL), and compare with an observationally-motivated EBL model and with other available observational constraints.

Recent low-redshift observations have attempted to determine the star formation histories of elliptical galaxies by tracking correlations between the stellar population parameters (age and metallicity) and the structural parameters that enter the fundamental plane (size and velocity dispersion). These studies have found that velocity dispersion, rather than effective radius or dynamical mass, is the main predictor of a galaxy's stellar age and metallicity. In this work, we apply an analytic model that predicts the structural properties of remnants formed in major mergers to progenitor disk galaxies with properties taken from two different semi-analytic models. We predict the effective radius, velocity dispersion, luminosity, age, and metallicity of the merger remnants, enabling us to compare directly to observations of early-type galaxies. While we find a tight correlation between age and velocity dispersion, we find a stronger dependence of age and metallicity on effective radius than observations report. The correlations arise as a result of the dependence of gas fraction, age, and metallicity on the stellar mass in the progenitor disk galaxies. These dependences induce a rotation in the radius-velocity plane between the correlations with effective radius and circular velocity in the disk galaxy progenitors, and the correlations with effective radius and velocity dispersion in the elliptical galaxy remnants. The differences between our results and those from observations suggest that major mergers alone cannot produce the observed lack of correlation between effective radius and stellar population parameters. Simulations have suggested that subsequent minor mergers introduce scatter in the effective radius while leaving the velocity dispersion essentially unchanged. Incorporating such minor mergers into the model may, then, bring the simulations into closer agreement with observations.

We present all-sky simulated Fermi maps of gamma-rays from dark matter decay and annihilation in the Local Universe. The dark matter distribution is obtained from a constrained cosmological simulation of the neighboring large-scale structure provided by the CLUES project. The dark matter fields of density and density squared are then taken as an input for the Fermi observation simulation tool to predict the gamma-ray photon counts that Fermi would detect in 5 years of all-sky survey for given dark matter models. Signal-to-noise sky maps have also been obtained by adopting the current Galactic and isotropic diffuse background models released by the Fermi collaboration. We point out the possibility for Fermi to detect a dark matter gamma-ray signal in local extragalactic structures. In particular, we conclude here that Fermi observations of nearby clusters (e.g. Virgo and Coma) and filaments are expected to give stronger constraints on decaying dark matter compared to previous studies. As an example, we find a significant signal-to-noise ratio in dark matter models with a decay rate fitting the positron excess as measured by PAMELA. This is the first time that dark matter filaments are shown to be promising targets for indirect detection of dark matter. On the other hand, the prospects for detectability of annihilating dark matter in local extragalactic structures are less optimistic even with extreme cross-sections. We make the dark matter density and density squared maps available online at http://www.clues-project.org/articles/darkmattermaps.html

We address the internal support against total free-fall collapse of the giant clumps that form by violent gravitational instability in high-z disc galaxies. Guidance is provided by an analytic model, where the proto-clumps are cut from a rotating disc and collapse to equilibrium while preserving angular momentum. This model predicts prograde clump rotation. This is confirmed in hydro-AMR zoom-in simulations of galaxies in a cosmological context. In most high-z clumps, the centrifugal force dominates the support, R=Vrot^2/Vcirc^2 > 0.5, where Vrot is the rotation velocity and Vcirc is the circular velocity. The clump spin indeed tends to be in the sense of the global disc angular momentum, but substantial tilts are frequent. Most clumps are in Jeans equilibrium, with the rest of the support provided by turbulence. Simulations of isolated gas-rich discs that resolve the clump substructure reveal that the cosmological simulations may overestimate R by ~30%, but the dominance of rotational support at high-z is not a resolution artifact. In turn, isolated gas-poor disc simulations produce at z=0 smaller gaseous non-rotating transient clouds, indicating that the difference in rotational support is associated with the fraction of cold baryons in the disc. In our current cosmological simulations, the clump rotation velocity is typically Vrot~100 km/s, but when beam smearing of \geq 0.1 arcsec is imposed, the rotation signal is reduced to a small gradient of \leq 30 km/s/kpc across the clump. The velocity dispersion in the simulated clumps is comparable to the disc dispersion so it is expected to leave only a marginal signal. Retrograde minor-merging galaxies could lead to massive clumps that do not show rotation.Testable predictions of the scenario as simulated are that the mean stellar age of the clumps, and the stellar fraction, are declining linearly with distance from the disc center.

We present a new measurement of the mass of the Milky Way (MW) based on observed properties of its largest satellite galaxies, the Magellanic Clouds (MCs), and an assumed prior of a {\Lambda}CDM universe. The large, high-resolution Bolshoi cosmological simulation of this universe provides a means to statistically sample the dynamical properties of bright satellite galaxies in a large population of dark matter halos. The observed properties of the MCs, including their circular velocity, distance from the center of the MW, and velocity within the MW halo, are used to evaluate the likelihood that a given halo would have each or all of these properties; the posterior PDF for any property of the MW system can thus be constructed. This method provides a constraint on the MW virial mass, 1.2 +0.7 -0.3(stat.) +0.3 -0.4 (sys.) x 10^12 M\odot (68% confidence), which is consistent with recent determinations that involve very different assumptions. In addition, we calculate the posterior PDF for the density profile of the MW and its satellite accretion history. Although typical satellites of 10^12 M\odot halos are accreted over a wide range of epochs over the last 10 Gyr, we find a \sim72% probability that the Magellanic Clouds were accreted within the last Gyr, and a 50% probability that they were accreted together.

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.

This paper describes the Hubble Space Telescope imaging data products and data reduction procedures for the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). This survey is designed to document the evolution of galaxies and black holes at $z\sim1.5-8$, and to study Type Ia SNe beyond $z>1.5$. Five premier multi-wavelength sky regions are selected, each with extensive multiwavelength observations. The primary CANDELS data consist of imaging obtained in the Wide Field Camera 3 / infrared channel (WFC3/IR) and UVIS channel, along with the Advanced Camera for Surveys (ACS). The CANDELS/Deep survey covers \sim125 square arcminutes within GOODS-N and GOODS-S, while the remainder consists of the CANDELS/Wide survey, achieving a total of \sim800 square arcminutes across GOODS and three additional fields (EGS, COSMOS, and UDS). We summarize the observational aspects of the survey as motivated by the scientific goals and present a detailed description of the data reduction procedures and products from the survey. Our data reduction methods utilize the most up to date calibration files and image combination procedures. We have paid special attention to correcting a range of instrumental effects, including CTE degradation for ACS, removal of electronic bias-striping present in ACS data after SM4, and persistence effects and other artifacts in WFC3/IR. For each field, we release mosaics for individual epochs and eventual mosaics containing data from all epochs combined, to facilitate photometric variability studies and the deepest possible photometry. A more detailed overview of the science goals and observational design of the survey are presented in a companion paper.

The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) is designed to document the first third of galactic evolution, over the approximate redshift (z) range 8--1.5. It will image >250,000 distant galaxies using three separate cameras on the Hubble Space Telescope, from the mid-ultraviolet to the near-infrared, and will find and measure Type Ia supernovae at z>1.5 to test their accuracy as standardizable candles for cosmology. Five premier multi-wavelength sky regions are selected, each with extensive ancillary data. The use of five widely separated fields mitigates cosmic variance and yields statistically robust and complete samples of galaxies down to a stellar mass of 10^9 M_\odot to z \approx 2, reaching the knee of the ultraviolet luminosity function (UVLF) of galaxies to z \approx 8. The survey covers approximately 800 arcmin^2 and is divided into two parts. The CANDELS/Deep survey (5\sigma\ point-source limit H=27.7 mag) covers \sim 125 arcmin^2 within GOODS-N and GOODS-S. The CANDELS/Wide survey includes GOODS and three additional fields (EGS, COSMOS, and UDS) and covers the full area to a 5\sigma\ point-source limit of H \gtrsim 27.0 mag. Together with the Hubble Ultra Deep Fields, the strategy creates a three-tiered "wedding cake" approach that has proven efficient for extragalactic surveys. Data from the survey are nonproprietary and are useful for a wide variety of science investigations. In this paper, we describe the basic motivations for the survey, the CANDELS team science goals and the resulting observational requirements, the field selection and geometry, and the observing design. The Hubble data processing and products are described in a companion paper.

Very high energy (VHE, i.e. $\gtrsim 10$ GeV) photons from Gamma-Ray Bursts (GRBs), as high as 90 GeV in rest frame energy, have been detected by the $Fermi$ Large Area Telescope (LAT). This provides hope for a high statistics GRB detection by a ground-based instrument in the VHE domain. We here report on our estimates of the expected GRB detection rate and $\gamma$-ray rate (in case of detection) for the next generation ground-based Imaging Air Cherenkov Telescopes (IACTs): the Cherenkov Telescope Array (CTA). Moreover, we investigated the effect critical design parameters of CTA have on our predictions and we compare the performance of various observing strategies for the array pointing. Our estimations are based on a phenomenological model which uses temporal and spectral information from GRBs detected by $Fermi$-LAT as well as other instruments operating at lower energy. While detection of VHE emission from GRBs has eluded ground-based instruments thus far, our results suggest it should be within reach of CTA with a rate between 0.35 and 1.6 GRBs/year depending on the characteristics of the true GRB population and the performance that CTA will eventually achieve.

It has long been regarded as difficult for a cosmological model to account simultaneously for the galaxy luminosity, mass, and velocity distributions. We revisit this issue using a modern compilation of observational data along with the best available large-scale cosmological simulation of dark matter. We find that the standard cosmological model, used in conjunction with halo abundance matching (HAM) and simple dynamical corrections, fits all basic statistics of galaxies with circular velocities Vcirc > 80 km/s. Our observational constraint is the luminosity-velocity relation which allows all types of galaxies to be included. We have compiled data for a variety of galaxies ranging from dwarf irregulars to giant ellipticals. The data present a clear monotonic luminosity-velocity relation from 50 km/s to 500 km/s, with a bend below 80 km/s and a systematic offset between late- and early-type galaxies. For comparison to theory, we employ our LCDM "Bolshoi" simulation of dark matter, which has unprecedented mass and force resolution. We use halo abundance matching to assign rank-ordered galaxy luminosities to the dark matter halos. The resulting predictions for the luminosity-velocity relation are in excellent agreement with the available data on both early-type and late-type galaxies for the luminosity range from Mr = -14-22. We also compare our predictions for the "cold" baryon mass (i.e., stars and cold gas) of galaxies as a function of circular velocity with the available observations, again finding a very good agreement. The predicted circular velocity function is in agreement with the galaxy velocity function for 80-400 km/s. However, we find that the dark matter halos with Vcirc < 80 km/s are much more abundant than observed galaxies with the same Vcirc . We find that the two-point correlation function of galaxies in our model matches very well the results from the SDSS.

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

We present the first results from the new Bolshoi N-body cosmological LCDM simulation that uses cosmological parameters favored by current observations. The Bolshoi simulation was done in a volume 250Mpc on a side using 8billion particles with mass and force resolution adequate to follow subhalos down to a completeness limit of Vcirc=50km/ s circular velocity. Using excellent statistics of halos and subhalos (10M at every moment and 50M over the whole history) we present accurate approximations for statistics such as the halo mass function, the concentrations for distinct halos and subhalos, abundance of halos as function of their circular velocity, the abundance and the spatial distribution of subhalos. We find that at high redshifts the concentration falls to a minimum of about 3.8 and then rises slightly for higher values of halo mass. We find that while the Sheth-Tormen approximation for the mass function of halos found by spherical overdensity is accurate at low redshifts, it over-predicts the abundance of halos by nearly an order of magnitude by z=10. We find that the number of subhalos scales with the circular velocity of the host halo as Vhost**0.5, and that subhalos have nearly the same radial distribution as dark matter particles at radii 0.3-2 times the host halo virial radius. The subhalo velocity function n(>V) behaves as V**(-3). We give normalization of this relation for different masses and redshifts. Finally, we use an abundance-matching procedure to assign r-band luminosities to dark matter halos as a function of halo Vcirc, and find that the luminosity-velocity relation is in remarkably good agreement with the observed Tully-Fisher relation for galaxies in the range 50-200km/s.

Calculating the galaxy merger rate requires both a census of galaxies identified as merger candidates, and a cosmologically-averaged `observability' timescale T_obs(z) for identifying galaxy mergers. While many have counted galaxy mergers using a variety of techniques, T_obs(z) for these techniques have been poorly constrained. We address this problem by calibrating three merger rate estimators with a suite of hydrodynamic merger simulations and three galaxy formation models. We estimate T_obs(z) for (1) close galaxy pairs with a range of projected separations, (2) the morphology indicator G-M20, and (3) the morphology indicator asymmetry A. Then we apply these timescales to the observed merger fractions at z < 1.5 from the recent literature. When our physically-motivated timescales are adopted, the observed galaxy merger rates become largely consistent. The remaining differences between the galaxy merger rates are explained by the differences in the range of mass-ratio measured by different techniques and differing parent galaxy selection. The major merger rate per unit co-moving volume for samples selected with constant number density evolves much more strongly with redshift (~ (1+z)^(+3.0 \pm 1.1)) than samples selected with constant stellar mass or passively evolving luminosity (~ (1+z)^(+0.1 \pm 0.4)). We calculate the minor merger rate (1:4 < M_{sat}/M_{primary} <~ 1:10) by subtracting the major merger rate from close pairs from the `total' merger rate determined by G-M20. The implied minor merger rate is ~3 times the major merger rate at z ~ 0.7, and shows little evolution with redshift.

Hydro cosmological simulations reveal that massive galaxies at high redshift are fed by long narrow streams of merging galaxies and a smoother component of cold gas. We post-process seven high-resolution simulated galaxies with radiative transfer to study the absorption characteristics of the gas in galaxies and streams, in comparison with the statistics of observed absorption-line systems. We find that much of the stream gas is ionized by UV radiation from background and local stellar sources, but still optically thick (N_HI > 10^17 cm^-2) so that the streams appear as Lyman-limit systems (LLSs). At z>3, the fraction of neutral gas in streams becomes non-negligible, giving rise to damped Lyman-alpha (DLA) absorbers as well. The gas in the central and incoming galaxies remains mostly neutral, responsible for DLAs. Within one (two) virial radii, the covering factor of optically thick gas is <25% (10%) for LLSs and <5% (1%) for DLAs, slowly declining with time following the universal expansion. Nevertheless, galaxies and their cold streams in the studied mass range, M_vir = 10^10 - 10^12 Msun, account for >30% of the observed absorbers in the foreground of quasars, the rest possibly arising from smaller galaxies or the intergalactic medium. The mean metallicity in the streams is ~1% solar, much lower than in the galaxies. The simulated galaxies reproduce the Lyalpha-absorption equivalent widths observed around Lyman-break galaxies, but they severely underpredict the equivalent widths in metal lines, suggesting that the latter may arise from outflows. We conclude that the observed metal-poor LLSs are likely detections of the predicted cold streams. Revised analysis of the observed LLSs kinematics and simulations with more massive outflows in conjunction with the inflows may enable a clearer distinction between the signatures of the various gas modes.

Data from (non-) attenuation of gamma rays from active galactic nuclei (AGN) and gamma ray bursts (GRBs) give upper limits on the extragalactic background light (EBL) from the UV to the mid-IR that are only a little above the lower limits from observed galaxies. These upper limits now rule out some EBL models and purported observations, with improved data likely to provide even stronger constraints. We present EBL calculations both based on multiwavelength observations of thousands of galaxies and also based on semi-analytic models, and show that they are consistent with these lower limits from observed galaxies and with the gamma-ray upper limit constraints. Such comparisons "close the loop" on cosmological galaxy formation models, since they account for all the light, including that from galaxies too faint to see. We compare our results with those of other recent works, and discuss the implications of these new EBL calculations for gamma ray attenuation. Catching a few GRBs with groundbased atmospheric Cherenkov Telescope (ACT) arrays or water Cherenkov detectors could provide important new constraints on the high-redshift star formation history of the universe.

There are strong correlations between the three structural properties of elliptical galaxies -- stellar mass, velocity dispersion and size -- in the form of a tight "fundamental plane" and a "scaling relation" between each pair. Major mergers of disk galaxies are assumed to be a mechanism for producing ellipticals, but semi-analytic galaxy formation models (SAM) have encountered apparent difficulties in reproducing the observed slope and scatter of the size-mass relation. We study the scaling relations of merger remnants using progenitor properties from two SAMs. We apply a simple merger model that includes gas dissipation and star formation based on theoretical considerations and simulations. Combining the SAMs and the merger model allows calculation of the structural properties of the remnants of major mergers that enter the population of elliptical galaxies at a given redshift. Without tuning the merger model parameters for each SAM, the results roughly match the slope and scatter in the observed scaling relations and their evolution in the redshift range $z=0-3$. Within this model, the observed scaling relations, including the tilt of the fundamental plane relative to the virial plane, result primarily from the decrease of gas fraction with increasing progenitor mass. The scatter in the size-mass relation of the remnants is reduced from that of the progenitors because of a correlation between progenitor size and gas fraction at a given mass.

We study the concentration of dark matter halos and its evolution in N-body simulations of the standard LCDM cosmology. The results presented in this paper are based on 4 large N-body simulations with about 10 billion particles each: the Millennium-I and II, Bolshoi, and MultiDark simulations. The MultiDark (or BigBolshoi) simulation is introduced in this paper. This suite of simulations with high mass resolution over a large volume allows us to compute with unprecedented accuracy the concentration over a large range of scales (about six orders of magnitude in mass), which constitutes the state-of-the-art of our current knowledge on this basic property of dark matter halos in the LCDM cosmology. We find that there is consistency among the different simulation data sets. We confirm a novel feature for halo concentrations at high redshifts: a flattening and upturn with increasing mass. The concentration c(M,z) as a function of mass and the redshift and for different cosmological parameters shows a remarkably complex pattern. However, when expressed in terms of the linear rms fluctuation of the density field sigma(M,z), the halo concentration c(sigma) shows a nearly-universal simple U-shaped behaviour with a minimum at a well defined scale at sigma=0.71. Yet, some small dependences with redshift and cosmology still remain. At the high-mass end (sigma < 1) the median halo kinematic profiles show large signatures of infall and highly radial orbits. This c-sigma(M,z) relation can be accurately parametrized and provides an analytical model for the dependence of concentration on halo mass. When applied to galaxy clusters, our estimates of concentrations are substantially larger -- by a factor up to 1.5 -- than previous results from smaller simulations, and are in much better agreement with results of observations. (abridged)

The extragalactic background light (EBL) is of fundamental importance both for understanding the entire process of galaxy evolution and for gamma-ray astronomy. However, the overall spectrum of the EBL between 0.1 and 1000 microns has never been determined directly neither from observed luminosity functions (LFs), over a wide redshift range, nor from any multiwavelength observation of galaxy spectral energy distributions (SEDs). The evolving, overall spectrum of the EBL is derived here utilizing a novel method based on observations only. The changing fractions of quiescent galaxies, star-forming galaxies, starburst galaxies and active galactic nucleus (AGN) galaxies from redshift 0.2 to 1 are estimated, and two alternative extrapolations of SED types to higher redshifts are considered. This allows calculation of the evolving EBL. The EBL uncertainties in our modelling based directly on the data are quantified, and their consequences for attenuation of very-high-energy gamma-rays due to pair production on the EBL are discussed. It is concluded that the EBL seems well constrained from the UV to the mid-IR at an intensity level roughly matching galaxy count data. Independent efforts from IR and gamma-ray astronomy are needed in order to reduce the uncertainties in the far-IR.

Fermi has shown GRBs to be a source of >10 GeV photons. We present an estimate of the detection rate of GRBs with a next generation Cherenkov telescope. Our predictions are based on the observed properties of GRBs detected by Fermi, combined with the spectral properties and redshift determinations for the bursts population by instruments operating at lower energies. While detection of VHE emission from GRBs has eluded ground-based instruments thus far, our results suggest that ground-based detection may be within reach of the proposed Cherenkov Telescope Array (CTA), albeit with a low rate, 0.25 - 0.5/yr. Such a detection would help constrain the emission mechanism of gamma-ray emission from GRBs. Photons at these energies from distant GRBs are affected by the UV-optical background light, and a ground-based detection could also provide a valuable probe of the Extragalactic Background Light (EBL) in place at high redshift.

The extragalactic background light (EBL) is of fundamental importance both for understanding the entire process of galaxy evolution and for gamma-ray astronomy, but the overall spectrum of the EBL between 0.1-1000 microns has never been determined directly from galaxy spectral energy distribution (SED) observations over a wide redshift range. The evolving, overall spectrum of the EBL is derived here utilizing a novel method based on observations only. This is achieved from the observed evolution of the rest-frame K-band galaxy luminosity function up to redshift 4 (Cirasuolo et al. 2010), combined with a determination of galaxy SED-type fractions. These are based on fitting SWIRE templates to a multiwavelength sample of about 6000 galaxies in the redshift range from 0.2 to 1 from the All-wavelength Extended Groth Strip International Survey (AEGIS). The changing fractions of quiescent galaxies, star-forming galaxies, starburst galaxies and AGN galaxies in that redshift range are estimated, and two alternative extrapolations of SED-types to higher redshifts are considered. This allows calculation of the evolution of the luminosity densities from the UV to the IR, the evolving star formation rate density of the universe, the evolving contribution to the bolometric EBL from the different galaxy populations including AGN galaxies and the buildup of the EBL. Our EBL calculations are compared with those from a semi-analytic model, from another observationally-based model and observational data. The EBL uncertainties in our modeling based directly on the data are quantified, and their consequences for attenuation of very high energy gamma-rays due to pair production on the EBL are discussed. It is concluded that the EBL is well constrained from the UV to the mid-IR, but independent efforts from infrared and gamma-ray astronomy are needed in order to reduce the uncertainties in the far-IR.

We describe the effect of AGN light on host galaxy optical and UV-optical colours, as determined from X-ray-selected AGN host galaxies at z~1, and compare the AGN host galaxy colours to those of a control sample matched to the AGN sample in both redshift and stellar mass. We identify as X-ray-selected AGNs 8.7 +4/-3 per cent of the red-sequence control galaxies, 9.8 +/-3 per cent of the blue-cloud control galaxies, and 14.7 +4/-3 per cent of the green-valley control galaxies. The nuclear colours of AGN hosts are generally bluer than their outer colours, while the control galaxies exhibit redder nuclei. AGNs in blue-cloud host galaxies experience less X-ray obscuration, while AGNs in red-sequence hosts have more, which is the reverse of what is expected from general considerations of the interstellar medium. Outer and integrated colours of AGN hosts generally agree with the control galaxies, regardless of X-ray obscuration, but the nuclear colours of unobscured AGNs are typically much bluer, especially for X-ray luminous objects. Visible point sources are seen in many of these, indicating that the nuclear colours have been contaminated by AGN light and that obscuration of the X-ray radiation and visible light are therefore highly correlated. Red AGN hosts are typically slightly bluer than red-sequence control galaxies, which suggests that their stellar populations are slightly younger. We compare these colour data to current models of AGN formation. The unexpected trend of less X-ray obscuration in blue-cloud galaxies and more in red-sequence galaxies is problematic for all AGN feedback models, in which gas and dust is thought to be removed as star formation shuts down. [See paper for full abstract.]

We use high-resolution cosmological hydrodynamical AMR simulations to predict the characteristics of La emission from the cold gas streams that fed galaxies in massive haloes at high redshift. The La luminosity in our simulations is powered by the release of gravitational energy as gas flows from the intergalactic medium into the halo potential wells. The UV background contributes only <20% to the gas heating. The La emissivity is due primarily to electron-impact excitation cooling radiation in gas ~2x10^4K. We calculate the La emissivities assuming collisional ionisation equilibrium (CIE) at all gas temperatures. The simulated streams are self-shielded against the UV background, so photoionisation and recombination contribute negligibly to the La line formation. We produce theoretical maps of the La surface brightnesses, assuming that ~85% of the La photons are directly observable. We find that typical haloes of mass Mv~10^12-13 Msun at z~3 emit as La blobs (LABs) with luminosities 10^43-44 erg/s. Most of the La comes from the extended narrow, partly clumpy, inflowing, cold streams that feed the growing galaxies. The predicted LAB morphology is therefore irregular, with dense clumps and elongated extensions. The linewidth is expected to range from 10^2 to more than 10^3 km/s with a large variance. The typical La surface brightness profile is proportional to r^-1.2 where r is the distance from the halo centre. Our simulated LABs are similar in luminosity, morphology and extent to the observed LABs, with distinct kinematic features. The predicted La luminosity function is consistent with observations, and the predicted areas and linewidths roughly recover the observed scaling relations. This mechanism for producing LABs appears inevitable in many high-z galaxies. Some of the LABs may thus be regarded as direct detections of the cold streams that drove galaxy evolution at high z.

We assess the effects of simulated active galactic nuclei (AGNs) on the colour and morphology measurements of their host galaxies. To test the morphology measurements, we select a sample of galaxies not known to host AGNs and add a series of point sources scaled to represent specified fractions of the observed V band light detected from the resulting systems; we then compare morphology measurements of the simulated systems to measurements of the original galaxies. AGN contributions >20 per cent bias most of the morphology measurements tested, though the extent of the apparent bias depends on the morphological characteristics of the original galaxies. We test colour measurements by adding to non-AGN galaxy spectra a quasar spectrum scaled to contribute specified fractions of the rest-frame B band light detected from the resulting systems. A quasar fraction of 5 per cent can move the NUV-r colour of an elliptical galaxy from the UV-optical red sequence to the green valley, and 20 per cent can move it into the blue cloud. Combining the colour and morphology results, we find that a galaxy/AGN system with an AGN contribution >20 per cent may appear bluer and more bulge-dominated than the underlying galaxy. We conclude that (1) bulge-dominated, E/S0/Sa, and early-type morphology classifications are accurate for red AGN host galaxies and may be accurate for blue host galaxies, unless the AGN manifests itself as a well-defined point source; and (2) although highly unobscured AGNs, such as the quasar used for our experiments, can significantly bias the measured colours of AGN host galaxies, it is possible to identify such systems by examining optical images of the hosts for the presence of a point source and/or measuring the level of nuclear obscuration.

When calculating the infrared spectral energy distributions (SEDs) of galaxies in radiation-transfer models, the calculation of dust grain temperatures is generally the most time-consuming part of the calculation. Because of its highly parallel nature, this calculation is perfectly suited for massively parallel general-purpose Graphics Processing Units (GPUs). This paper presents an implementation of the calculation of dust grain equilibrium temperatures on GPUs in the Monte-Carlo radiation transfer code Sunrise, using the CUDA API. The GPU can perform this calculation 69 times faster than the 8 CPU cores, showing great potential for accelerating calculations of galaxy SEDs.

The majority of galaxy mergers are expected to be minor mergers. The observational signatures of minor mergers are not well understood, thus there exist few constraints on the minor merger rate. This paper seeks to address this gap in our understanding by determining if and when minor mergers exhibit disturbed morphologies and how they differ from the morphology of major mergers. We simulate a series of unequal-mass moderate gas-fraction disc galaxy mergers. With the resulting g-band images, we determine how the time-scale for identifying galaxy mergers via projected separation and quantitative morphology (the Gini coefficient G, asymmetry A, and the second-order moment of the brightest 20% of the light M20) depends on the merger mass ratio, relative orientations and orbital parameters. We find that G-M20 is as sensitive to 9:1 baryonic mass ratio mergers as 1:1 mergers, with observability time-scales ~ 0.2-0.4 Gyr. In contrast, asymmetry finds mergers with baryonic mass ratios between 4:1 and 1:1 (assuming local disc galaxy gas-fractions). Asymmetry time-scales for moderate gas-fraction major disc mergers are ~ 0.2-0.4 Gyr, and less than 0.06 Gyr for moderate gas-fraction minor mergers. The relative orientations and orbits have little effect on the time-scales for morphological disturbances. Observational studies of close pairs often select major mergers by choosing paired galaxies with similar luminosities and/or stellar masses. Therefore, the various ways of finding galaxy mergers (G-M20, A, close pairs) are sensitive to galaxy mergers of different mass ratios. By comparing the frequency of mergers selected by different techniques, one may place empirical constraints on the major and minor galaxy merger rates.

Gas-rich galaxy mergers are more easily identified by their disturbed morphologies than mergers with less gas. Because the typical gas fraction of galaxy mergers is expected to increase with redshift, the under-counting of low gas-fraction mergers may bias morphological estimates of the evolution of galaxy merger rate. To understand the magnitude of this bias, we explore the effect of gas fraction on the morphologies of a series of simulated disc galaxy mergers. With the resulting g-band images, we determine how the time-scale for identifying major and minor galaxy mergers via close projected pairs and quantitative morphology (the Gini coefficient G, the second-order moment of the brightest 20% of the light M20, and asymmetry A) depends on baryonic gas fraction f(gas). Strong asymmetries last significantly longer in high gas-fraction mergers of all mass ratios, with time-scales ranging from >= 300 Myr for f(gas) ~ 20% to >= 1 Gyr for f(gas) ~ 50%. Therefore the strong evolution with redshift observed in the fraction of asymmetric galaxies may reflect evolution in the gas properties of galaxies rather than the global galaxy merger rate. On the other hand, the time-scale for identifying a galaxy merger via G-M20 is weakly dependent on gas-fraction (~ 200-400 Myr), consistent with the weak evolution observed for G-M20 mergers.

There is a large observational scatter toward low velocities in the stellar mass Tully-Fisher relation if disturbed and compact objects are included. However, this scatter can be eliminated if one replaces rotation velocity with $\rm S_{\rm 0.5}$, a quantity that includes a velocity dispersion term added in quadrature with the rotation velocity. In this work we use a large suite of hydrodynamic N-body galaxy merger simulations to explore a possible mechanism for creating the observed relations. Using mock observations of the simulations, we test for the presence of observational effects and explore the relationship between $\rm S_{\rm 0.5}$ and intrinsic properties of the galaxies. We find that galaxy mergers can explain the scatter in the TF as well as the tight $\rm S_{\rm 0.5}$-stellar mass relation. Furthermore, $\rm S_{\rm 0.5}$ is correlated with the total central mass of a galaxy, including contributions due to dark matter.

The abundance of dark matter satellites and subhalos, the existence of density cusps at the centers of dark matter halos, and problems producing realistic disk galaxies in simulations are issues that have raised concerns about the viability of the standard cold dark matter (LambdaCDM) scenario for galaxy formation. This article reviews these issues, and considers the implications for cold vs. various varieties of warm dark matter (WDM). The current evidence appears to be consistent with standard LambdaCDM, although improving data may point toward a rather tepid version of LambdaWDM -- tepid since the dark matter cannot be very warm without violating observational constraints.

The four lectures that I gave in the XIII Ciclo de Cursos Especiais at the National Observatory of Brazil in Rio in October 2008 were (1) a brief history of dark matter and structure formation in a LambdaCDM universe; (2) challenges to LambdaCDM on small scales: satellites, cusps, and disks; (3) data on galaxy evolution and clustering compared with simulations; and (4) semi-analytic models. These lectures, themselves summaries of much work by many people, are summarized here briefly.

Gamma rays from extragalactic sources are attenuated by pair-production interactions with diffuse photons of the extragalactic background light (EBL). Gamma-ray bursts (GRBs) are a source of high-redshift photons above 10 GeV, and could be therefore useful as a probe of the evolving UV background radiation. In this paper, we develop a simple phenomenological model for the number and redshift distribution of gamma-ray bursts that can be seen at GeV energies with the Fermi satellite and MAGIC atmospheric Cherenkov telescope. We estimate the observed number of gamma rays per year, and show how this result is modified by considering interactions with different realizations of the evolving EBL. We also discuss the bright Fermi GRB 080916C in the context of this model. We find that the LAT on Fermi can be expected to see a small number of photons above 10 GeV each year from distant GRBs. Annual results for ground-based instruments like MAGIC are highly variable due to the low duty cycle and sky coverage of the telescope. However, successfully viewing a bright or intermediate GRB from the ground could provide hundreds of photons from high redshift, which would almost certainly be extremely useful in constraining both GRB physics and the high-redshift EBL.

We present new calculations of the evolving UV background out to the epoch of cosmological reionization and make predictions for the amount of GeV gamma-ray attenuation by electron-positron pair production. Our results are based on recent semi-analytic models of galaxy formation, which provide predictions of the dust-extinguished UV radiation field due to starlight, and empirical estimates of the contribution due to quasars. We account for the reprocessing of ionizing photons by the intergalactic medium. We test whether our models can reproduce estimates of the ionizing background at high redshift from flux decrement analysis and proximity effect measurements from quasar spectra, and identify a range of models that can satisfy these constraints. Pair-production against soft diffuse photons leads to a spectral cutoff feature for gamma rays observed between 10 and 100 GeV. This cutoff varies with redshift and the assumed star formation and quasar evolution models. We find only negligible amounts of absorption for gamma rays observed below 10 GeV for any emission redshift. With observations of high-redshift sources in sufficient numbers by the Fermi Gamma-ray Space Telescope and new ground-based instruments it should be possible to constrain the extragalactic background light in the UV and optical portion of the spectrum.

One of the main themes in extragalactic astronomy for the next decade will be the evolution of galaxies over cosmic time. Many future observatories, including JWST, ALMA, GMT, TMT and E-ELT will intensively observe starlight over a broad redshift range, out to the dawn of the modern Universe when the first galaxies formed. It has, however, become clear that the properties and evolution of galaxies are intimately linked to the growth of their central black holes. Understanding the formation of galaxies, and their subsequent evolution, will therefore be incomplete without similarly intensive observations of the accretion light from supermassive black holes (SMBH) in galactic nuclei. To make further progress, we need to chart the formation of typical SMBH at z>6, and their subsequent growth over cosmic time, which is most effectively achieved with X-ray observations. Recent technological developments in X-ray optics and instrumentation now bring this within our grasp, enabling capabilities fully matched to those expected from flagship observatories at longer wavelengths.

This White Paper to the National Academy of Sciences Astro2010 Decadal Review Committee outlines some of the outstanding questions regarding the assembly history of Massive Black Holes in the nuclei of galaxies and the revolutionary contributions anticipated in this field from low-frequency gravitational wave astronomy.

There is a vast menagerie of plausible candidates for the constituents of dark matter, both within and beyond extensions of the Standard Model of particle physics. Each of these candidates may have scattering (and other) cross section properties that are consistent with the dark matter abundance, BBN, and the most scales in the matter power spectrum; but which may have vastly different behavior at sub-galactic "cutoff" scales, below which dark matter density fluctuations are smoothed out. The only way to quantitatively measure the power spectrum behavior at sub-galactic scales at distances beyond the local universe, and indeed over cosmic time, is through probes available in multiply imaged strong gravitational lenses. Gravitational potential perturbations by dark matter substructure encode information in the observed relative magnifications, positions, and time delays in a strong lens. Each of these is sensitive to a different moment of the substructure mass function and to different effective mass ranges of the substructure. The time delay perturbations, in particular, are proving to be largely immune to the degeneracies and systematic uncertainties that have impacted exploitation of strong lenses for such studies. There is great potential for a coordinated theoretical and observational effort to enable a sophisticated exploitation of strong gravitational lenses as direct probes of dark matter properties. This opportunity motivates this white paper, and drives the need for: a) strong support of the theoretical work necessary to understand all astrophysical consequences for different dark matter candidates; and b) tailored observational campaigns, and even a fully dedicated mission, to obtain the requisite data.

The full spatio-chemo-dynamical structure of galaxies of all types and environments at low redshift provides a critical accompaniment to observations of galaxy formation at high redshift. The next decade brings the observational opportunity to strongly constrain nearby galaxies' histories of star formation and assembly, especially in the spheroids that comprise the large majority of the stellar mass in the Universe but have until now been difficult to study. In order to constrain the pathways to building up the spheroidal "red-sequence", various standard techniques in photometry and spectroscopy, particularly with resolved tracer populations like globular clusters and planetary nebulae, can be scaled up to comprehensive surveys as improved wide-field instrumentation is increasingly available. At the same time, progress in adaptive optics on giant telescopes could for the first time permit deep, resolved photometric and spectroscopic analysis of large samples of individual stars in these systems, thereby revolutionizing galaxy studies. Strong theoretical support is needed in order to understand the new observational constraints via detailed modeling and self-consistent simulations of star and galaxy formation throughout cosmic time.

The abundance of dark matter satellites and subhalos, the existence of density cusps at the centers of dark matter halos, and problems producing realistic disk galaxies in simulations are issues that have raised concerns about the viability of the standard cold dark matter (LambdaCDM) scenario for galaxy formation. This talk reviews these issues, and considers the implications for cold vs. various varieties of warm dark matter (WDM). The current evidence appears to be consistent with standard LambdaCDM, although improving data may point toward a rather tepid version of LambdaWDM - tepid since the dark matter cannot be very warm without violating observational constraints.

(Brief Summary) What is the total radiative content of the Universe since the epoch of recombination? The extragalactic background light (EBL) spectrum captures the redshifted energy released from the first stellar objects, protogalaxies, and galaxies throughout cosmic history. Yet, we have not determined the brightness of the extragalactic sky from UV/optical to far-infrared wavelengths with sufficient accuracy to establish the radiative content of the Universe to better than an order of magnitude. Among many science topics, an accurate measurement of the EBL spectrum from optical to far-IR wavelengths, will address: What is the total energy released by stellar nucleosynthesis over cosmic history? Was significant energy released by non-stellar processes? Is there a diffuse component to the EBL anywhere from optical to sub-millimeter? When did first stars appear and how luminous was the reionization epoch? Absolute optical to mid-IR EBL spectrum to an astrophysically interesting accuracy can be established by wide field imagingat a distance of 5 AU or above the ecliptic plane where the zodiacal foreground is reduced by more than two orders of magnitude.

Attenuation of high--energy gamma rays by pair--production with UV, optical and IR background photons provides a link between the history of galaxy formation and high--energy astrophysics. We present results from our latest semi-analytic models (SAMs), based upon a $\Lambda$CDM hierarchical structural formation scenario and employing all ingredients thought to be important to galaxy formation and evolution, as well as reprocessing of starlight by dust to mid- and far-IR wavelengths. Our models also use results from recent hydrodynamic galaxy merger simulations. These latest SAMs are successful in reproducing a large variety of observational constraints such as number counts, luminosity and mass functions, and color bimodality. We have created 2 models that bracket the likely ranges of galaxy emissivities, and for each of these we show how the optical depth from pair--production is affected by redshift and gamma-ray energy. We conclude with a discussion of the implications of our work, and how the burgeoning science of gamma-ray astronomy will continue to help constrain cosmology.

We present two models for the cosmological UV background light, and calculate the opacity of GeV gamma--rays out to redshift 9. The contributors to the background include 2 possible quasar emissivities, and output from star--forming galaxies as determined by recent a semi--analytic model (SAM) of structure formation. The SAM used in this work is based upon a hierarchical build-up of structure in a $\Lambda$CDM universe and is highly successful in reproducing a variety of observational parameters. Above 1 Rydberg energy, ionizing radiation is subject to reprocessing by the IGM, which we treat using our radiative transfer code, CUBA. The two models for quasar emissivity differing above z = 2.3 are chosen to match the ionization rates observed using flux decrement analysis and the higher values of the line-of-sight proximity effect. We also investigate the possibility of a flat star formation rate density at z $>5$. We conclude that observations of gamma--rays from 10 to 100 GeV by Fermi (GLAST) and the next generation of ground based experiments should confirm a strongly evolving opacity from $1<$ z $<4$. Observation of attenuation in the spectra of gamma--ray bursts at higher redshift could constrain emission of UV radiation at these early times, either from a flat or increasing star-formation density or an unobserved population of sources.

A key obstacle to understanding the galaxy merger rate and its role in galaxy evolution is the difficulty in constraining the merger properties and time-scales from instantaneous snapshots of the real universe.The most common way to identify galaxy mergers is by morphology, yet current theoretical calculations of the time-scales for galaxy disturbances are quite crude. We present a morphological analysis of a large suite of GADGET N-Body/hydro-dynamical equal-mass gas-rich disc galaxy mergers which have been processed through the Monte-Carlo radiative transfer code SUNRISE. With the resulting images, we examine the dependence of quantitative morphology (G, M20, C, A) in the SDSS g-band on merger stage, dust, viewing angle, orbital parameters, gas properties, supernova feedback, and total mass. We find that mergers appear most disturbed in G-M20 and asymmetry at the first pass and at the final coalescence of their nuclei, but can have normal quantitative morphologies at other merger stages. The merger observability time-scales depend on the method used to identify the merger as well as the gas fraction, pericentric distance, and relative orientation of the merging galaxies. Enhanced star formation peaks after and lasts significantly longer than strong morphological disturbances. Despite their massive bulges, the majority of merger remnants appear disc-like and dusty in g-band light because of the presence of a low-mass star-forming disc.

We study the effects of dust in hydrodynamic simulations of spiral galaxies when different radial metallicity gradients are assumed. SUNRISE, a Monte-Carlo radiative-transfer code, is used to make detailed calculations of the internal extinction of disk galaxies caused by their dust content. SUNRISE is used on eight different Smooth Particle Hydrodynamics (SPH) simulations of isolated spiral galaxies. These galaxies vary mainly in mass and hence luminosity, spanning a range in luminosities from -16 to -22 magnitudes in the B band. We focus on the attenuation in different wavelength bands as a function of the disk inclination and the luminosity of the models, and compare this to observations. Observations suggest different metallicity gradients for galaxies of different luminosities. These metallicity gradients were explored in our different models, finding that the resulting dust attenuation matches observations for edge-on galaxies, but do not show a linear behaviour in log axis ratio as some observations have suggested. A quadratic law describing the dependence of attenuation on inclination, as proposed by more recent observations, reconciles the attenuation of the simulations at intermediate inclinations with observations. We also compare the total infrared-to-ultraviolet flux ratios for the simulated galaxies with those of the SINGS sample and find general agreement. Finally we compare our results with those from simpler models that do not take into account structure such as spiral arms, finding that the inclusion of sub-structure on the size scale of spiral arms does not change conclusions about the attenuation dependence on inclination or wavelength.

We explore the role of AGN in establishing and/or maintaining the bimodal colour distribution of galaxies by quenching their star-formation and hence, causing their transition from the blue to the red cloud. Important tests for this scenario include (i) the X-ray properties of galaxies in the transition zone between the two clouds and (ii) the incidence of AGN in post-starbursts, i.e. systems observed shortly after (<1Gyr) the termination of their star-formation. We perform these tests by combining deep Chandra observations with multiwavelength data from the AEGIS survey. Stacking the X-ray photons at the positions of galaxies (0.4<z<0.9) not individually detected at X-ray wavelengths suggests a population of obscured AGN among sources in the transition zone and in the red cloud. Their mean X-ray and mid-IR properties are consistent with moderately obscured low-luminosity AGN, Compton thick sources or a mix of both. Morphologies show that major mergers are unlikely to drive the evolution of this population but minor interactions may play a role. The incidence of obscured AGN in the red cloud (both direct detections and stacking results) suggests that BH accretion outlives the termination of the star-formation. This is also supported by our finding that post-starburst galaxies at z~0.8 and AGN are associated, in agreement with recent results at low-z. A large fraction of post-starbursts and red cloud galaxies show evidence for at least moderate levels of AGN obscuration. This implies that if AGN outflows cause the colour transformation of galaxies, then some nuclear gas and dust clouds either remain unaffected or relax to the central galaxy regions after the quenching their star-formation.

We construct a physically motivated model for predicting the properties of the remnants of gaseous galaxy mergers, given the properties of the progenitors and the orbit. The model is calibrated using a large suite of SPH merger simulations. It implements generalized energy conservation while accounting for dissipative energy losses and star formation. The dissipative effects are evaluated from the initial gas fractions and from the orbital parameters via an "impulse" parameter, which characterizes the strength of the encounter. Given the progenitor properties, the model predicts the remnant stellar mass, half-mass radius, and velocity dispersion to an accuracy of 25%. The model is valid for both major and minor mergers. We provide an explicit recipe for semi-analytic models of galaxy formation.

Utilizing Chandra X-ray observations in the All-wavelength Extended Groth Strip International Survey (AEGIS) we identify 241 X-ray selected Active Galactic Nuclei (AGNs, L > 10^{42} ergs/s) and study the properties of their host galaxies in the range 0.4 < z < 1.4. By making use of infrared photometry from Palomar Observatory and BRI imaging from the Canada-France-Hawaii Telescope, we estimate AGN host galaxy stellar masses and show that both stellar mass and photometric redshift estimates (where necessary) are robust to the possible contamination from AGNs in our X-ray selected sample. Accounting for the photometric and X-ray sensitivity limits of the survey, we construct the stellar mass function of X-ray selected AGN host galaxies and find that their abundance decreases by a factor of ~2 since z~1, but remains roughly flat as a function of stellar mass. We compare the abundance of AGN hosts to the rate of star formation quenching observed in the total galaxy population. If the timescale for X-ray detectable AGN activity is roughly 0.5-1 Gyr--as suggested by black hole demographics and recent simulations--then we deduce that the inferred AGN "trigger" rate matches the star formation quenching rate, suggesting a link between these phenomena. However, given the large range of nuclear accretion rates we infer for the most massive and red hosts, X-ray selected AGNs may not be directly responsible for quenching star formation.

Using hydrodynamic simulations of disc-galaxy major mergers, we investigate the star formation history and remnant properties when various parametrizations of a simple stellar feedback model are implemented. The simulations include radiative cooling, a density-dependent star formation recipe and a model for feedback from massive stars. The feedback model stores supernova feedback energy within individual gas particles and dissipates this energy on a time-scale specified by two free parameters; tau_fb, which sets the dissipative time-scale, and n, which sets the effective equation of state in star-forming regions. Using a self-consistent disc galaxy, modelled after a local Sbc spiral, in both isolated and major-merger simulations, we investigate parametrizations of the feedback model that are selected with respect to the quiescent disc stability. These models produce a range of star formation histories that are consistent with the star formation relation found by Kennicutt. All major mergers produce a population of new stars that is highly centrally concentrated, demonstrating a distinct break in the r1/4 surface density profile, consistent with previous findings. The half-mass radius and one-dimensional velocity dispersion are affected by the feedback model used. Finally, we compare our results to those of previous simulations of star formation in disc-galaxy major mergers, addressing the effects of star formation normalization, the version of smoothed particle hydrodynamics (SPH) employed and assumptions about the interstellar medium.

We employ numerical simulations of galaxy mergers to explore the effect of galaxy mass ratio on merger--driven starbursts. Our numerical simulations include radiative cooling of gas, star formation, and stellar feedback to follow the interaction and merger of four disk galaxies. The galaxy models span a factor of 23 in total mass and are designed to be representative of typical galaxies in the local Universe. We find that the merger--driven star formation is a strong function of merger mass ratio, with very little, if any, induced star formation for large mass ratio mergers. We define a burst efficiency that is useful to characterize the merger--driven star formation and test that it is insensitive to uncertainties in the feedback parameterization. In accord with previous work we find that the burst efficiency depends on the structure of the primary galaxy. In particular, the presence of a massive stellar bulge stabilizes the disk and suppresses merger--driven star formation for large mass ratio mergers. Direct, co--planar merging orbits produce the largest tidal disturbance and yield that most intense burst of star formation. Contrary to naive expectations, a more compact distribution of gas or an increased gas fraction both decrease the burst efficiency. Owing to the efficient feedback model and the newer version of SPH employed here, the burst efficiencies of the mergers presented here are smaller than in previous studies.

We present the quantitative rest-frame B morphological evolution and galaxy merger fractions at 0.2 < z < 1.2 as observed by the All-wavelength Extended Groth Strip International Survey (AEGIS). We use the Gini coefficent and M_20 to identify major mergers and classify galaxy morphology for a volume-limited sample of 3009 galaxies brighter than 0.4 L_B^*, assuming pure luminosity evolution of 1.3 M_B per unit redshift. We find that the merger fraction remains roughly constant at 10 +/- 2% for 0.2 < z < 1.2. The fraction of E/S0/Sa increases from 21+/- 3% at z ~ 1.1 to 44 +/- 9% at z ~ 0.3, while the fraction of Sb-Ir decreases from 64 +/- 6% at z ~ 1.1 to 47 +/- 9% at z ~ 0.3. The majority of z < 1.2 Spitzer MIPS 24 micron sources with L(IR) > 10^11 L_sun are disk galaxies, and only ~ 15% are classified as major merger candidates. Edge-on and dusty disk galaxies (Sb-Ir) are almost a third of the red sequence at z ~ 1.1, while E/S0/Sa makeup over 90% of the red sequence at z ~ 0.3. Approximately 2% of our full sample are red mergers. We conclude (1) the galaxy merger rate does not evolve strongly between 0.2 < z < 1.2; (2) the decrease in the volume-averaged star-formation rate density since z ~ 1 is a result of declining star-formation in disk galaxies rather than a disappearing population of major mergers; (3) the build-up of the red sequence at z < 1 can be explained by a doubling in the number of spheroidal galaxies since z ~ 1.2.

We analyze star formation (SF) as a function of stellar mass (M*) and redshift z in the All Wavelength Extended Groth Strip International Survey (AEGIS). For 2905 field galaxies, complete to 10^10(10^10.8) Msun at z<0.7(1), with Keck spectroscopic redshifts out to z=1.1, we compile SF rates (SFR) from emission lines, GALEX, and Spitzer MIPS 24 micron photometry, optical-NIR M* measurements, and HST morphologies. Galaxies with reliable signs of SF form a distinct "main sequence (MS)", with a limited range of SFR at a given M* and z (1 sigma < +-0.3 dex), and log(SFR) approximately proportional to log(M*). The range of log(SFR) remains constant to z>1, while the MS as a whole moves to higher SFR as z increases. The range of SFR along the MS constrains the amplitude of episodic variations of SF, and the effect of mergers on SFR. Typical galaxies spend ~67(95)% of their lifetime since z=1 within a factor of <~ 2(4) of their average SFR at a given M* and z. The dominant mode of the evolution of SF since z~1 is apparently a gradual decline of the average SFR in most individual galaxies, not a decreasing frequency of starburst episodes, or a decreasing factor by which SFR are enhanced in starbursts. LIRGs at z~1 seem to mostly reflect the high SFR typical for massive galaxies at that epoch. The smooth MS may reflect that the same set of few physical processes governs star formation prior to additional quenching processes. A gradual process like gas exhaustion may play a dominant role.

We analyze star formation (SF) as a function of stellar mass (M_*) and redshift z in the All Wavelength Extended Groth Strip International Survey (AEGIS), for star-forming field galaxies with M_* >~ 10^10 M_sun out to z=1.1. The data indicate that the high specific SF rates (SFR) of many less massive galaxies do not represent late, irregular or recurrent, starbursts in evolved galaxies. They rather seem to reflect the onset (initial burst) of the dominant SF episode of galaxies, after which SF gradually declines on Gyr timescales to z=0 and forms the bulk of a galaxy's M_*. With decreasing mass, this onset of major SF shifts to decreasing z for an increasing fraction of galaxies (staged galaxy formation). This process may be an important component of the ``downsizing'' phenomenon. We find that the predominantly gradual decline of SFR (Noeske et al. 2007, this volume) can be reproduced by exponential SF histories (tau models), if less massive galaxies have systematically longer e-folding times tau, and a later onset of SF (z_f). Our model can provide a first parametrization of SFR as a function of M_* and z, and quantify mass-dependences of tau and z_f, from direct observations of M_* and SFR up to z>1. The observed evolution of SF in galaxies can plausibly reflect the dominance of gradual gas exhaustion. The data are also consistent with the history of cosmological accretion onto Dark Matter halos.

After briefly reviewing the good agreement between large-scale observations and the predictions of the now-standard CDM theory and problems with the MOND alternative, I summarize several of the main areas of possible disagreement between theory and observation: galaxy centers, dark matter substructure, angular momentum, and the sequence of cosmogony, updating earlier reviews [1]. All of these issues are sufficiently complicated that it is not yet clear how serious they are, but there is at least some reason to think that the problems will be resolved through a deeper understanding of the complicated "gastrophysics" of star formation and feedback from supernovae and AGN.

We visually and quantitatively determine the host galaxy morphologies of 94 intermediate redshift (0.2 < z < 1.2) active galactic nuclei (AGN), selected using Chandra X-ray and Spitzer mid-infrared data in the Extended Groth Strip. Using recently developed morphology measures, the second-order moment of the brightest 20% of a galaxy's flux (M_{20}) and the Gini coefficient, we find that X-ray-selected AGN mostly reside in E/S0/Sa galaxies (53^{+11}_{-10}%), while IR-selected AGN show no clear preference for host morphology. X-ray-selected AGN hosts are members of close pairs more often than the field population by a factor of 3.3+/-1.4, but most of these pair members appear to be undisturbed early-type galaxies and do not tend to show direct evidence of gravitational perturbations or interactions. Thus, the activation mechanism for AGN activity remains unknown, even for pair members.

Using a sample of 89 snapshots from 58 hydrodynamic binary galaxy major merger simulations, we find that stellar remnants are mostly oblate while dark matter halos are mostly prolate or triaxial. The stellar minor axis and the halo major axis are almost always nearly perpendicular. This can be understood by considering the influence of angular momentum and dissipation during the merger. If binary mergers of spiral galaxies are responsible for the formation of elliptical galaxies or some subpopulation thereof, these galaxies can be expected to be oblate and inhabit their halos with the predicted shapes and orientations. These predictions are potentially relevant to observational studies of weak gravitational lensing, where one must stack many optically aligned galaxies in order to determine the shape of the resulting stacked mass distribution. The simple relationship between the dark and luminous matter presented here can be used to guide the stacking of galaxies to minimize the information lost.

We apply a new approach to quantifying galaxy morphology and identifying galaxy mergers to the rest-frame far-ultraviolet images of 82 z ~ 4 Lyman break galaxies (LBGs) and 55 1.2 < z < 1.8 emission-line galaxies in the GOODS and Ultra Deep Fields. We compare the distributions of the Gini coefficient (G), second-order moment of the brightest 20% of galaxy light (M20), and concentration (C) for high-redshift and low-redshift galaxies with average signal to noise per pixel > 2.5 and Petrosian radii > 0.3 arcsec. Ten of the 82 LBGs have M20 >= -1.1 and possess bright double or multiple nuclei, implying a major-merger fraction of star-forming galaxies ~ 10-25% at M_{FUV} < -20, depending on our incompleteness corrections. Galaxies with bulge-like morphologies (G >= 0.55, M20 < -1.6) make up ~ 30% of the z ~ 4 LBG sample, while the remaining ~ 50% have G and M20 values higher than expected for smooth bulges and disks and may be star-forming disks, minor-mergers or post-mergers. The star-forming z ~ 1.5 galaxy sample has a morphological distribution which is similar to the UDF z ~ 4 LBGs, with an identical fraction of major-merger candidates but fewer spheroids. The observed morphological distributions are roughly consistent with current hierarchical model predictions for the major-merger rates and minor-merger induced starbursts at z ~ 1.5 and ~4. We also examine the rest-frame FUV-NUV and FUV-B colors as a function of morphology and find no strong correlations at either epoch.

(Abridged) We study predictions for galaxy cluster observables that can test the statistics of dark matter halo shapes expected in a flat LCDM universe. We present a simple analytical model for the prediction of cluster-scale X-ray observations, approximating clusters as isothermal systems in hydrostatic equilibrium, and dark matter haloes as ellipsoids with uniform axial ratios. We test the model against high-resolution, hydrodynamic cluster simulations to gauge its reliability. We find that this simple prescription does a good job of predicting the distribution of cluster X-ray ellipticities compared to the simulations as long as one focuses on cluster regions that are less sensitive to recent mergers. Based on this simple model, the distribution of cluster-size halo shapes expected in the concordance LCDM cosmology implies an X-ray ellipticity distribution with a mean of 0.32 +- 0.01 and a scatter of 0.14 +- 0.01 for the mass range (1-4)x10^{14} Msun/h. We find it important to include the mass dependence of halo shape to make comparisons to observational samples that contain many, very massive clusters. We analyse the systematics of four observational samples of cluster ellipticities and find that our results are statistically compatible with observations. In particular, we find remarkably good agreement between two recent ROSAT samples and LCDM predictions that DO NOT include gas cooling. We also test how well our analytical model can predict Sunyaev-Zel'dovich decrement maps and find that it is less successful although still useful; the model does not perform as well as a function of flux level in this case because of the changing triaxiality of dark matter haloes as a function of radial distance. Both this effect and the changing alignment of isodensity shells of dark matter haloes leave an imprint on cluster gas...

Using six high resolution dissipationless simulations with a varying box size in a flat LCDM universe, we study the mass and redshift dependence of dark matter halo shapes for M_vir = 9.0e11 - 2.0e14, over the redshift range z=0-3, and for two values of sigma_8=0.75 and 0.9. Remarkably, we find that the redshift, mass, and sigma_8 dependence of the mean smallest-to-largest axis ratio of halos is well described by the simple power-law relation <s> = (0.54 +- 0.02)(M_vir/M_*)^(-0.050 +- 0.003), where s is measured at 0.3 R_vir and the z and sigma_8 dependences are governed by the characteristic nonlinear mass, M_*=M_*(z,sigma_8). We find that the scatter about the mean s is well described by a Gaussian with sigma ~ 0.1, for all masses and redshifts. We compare our results to a variety of previous works on halo shapes and find that reported differences between studies are primarily explained by differences in their methodologies. We address the evolutionary aspects of individual halo shapes by following the shapes of the halos through ~100 snapshots in time. We determine the formation scalefactor a_c as defined by Wechsler et al. (2002) and find that it can be related to the halo shape at z = 0 and its evolution over time.

There is strong evidence that the mass in the Universe is dominated by dark matter, which exerts gravitational attraction but whose exact nature is unknown. In particular, all galaxies are believed to be embedded in massive haloes of dark matter. This view has recently been challenged by surprisingly low random stellar velocities in the outskirts of ordinary elliptical galaxies, which were interpreted as indicating a lack of dark matter (Mendez et al. 2001; Romanowsky et al. 2003). Here we show that the low velocities are in fact compatible with galaxy formation in dark-matter haloes. Using numerical simulations of disc-galaxy mergers, we find that the stellar orbits in the outer regions of the resulting ellipticals are very elongated. These stars were torn by tidal forces from their original galaxies during the first close passage and put on outgoing trajectories. The elongated orbits, combined with the steeply falling density profile of the observed tracers, explain the observed low velocities even in the presence of large amounts of dark matter. Projection effects when viewing a triaxial elliptical can lead to even lower observed velocities along certain lines of sight.

A new Monte-Carlo radiative-transfer code, Sunrise, is used in conjunction with hydrodynamic simulations of major galaxy mergers to calculate the effects of dust in such systems. The simulations are in good agreement with observations of dust absorption in starburst galaxies, and the dust has a profound effect on their appearance. The dust attenuation increases with luminosity such that at peak luminosities ~90% of the bolometric luminosity is absorbed by dust. In general, the detailed appearance of the merging event depends on the stage of the merger and the geometry of the encounter. The fraction of bolometric energy absorbed by the dust, however, is a robust quantity that can be predicted from the intrinsic properties bolometric luminosity, baryonic mass, star-formation rate, and metallicity of the system. This paper presents fitting formulae, valid over a wide range of masses and metallicities, from which the absorbed fraction of luminosity (and consequently also the infrared dust luminosity) can be predicted. The attenuation of the luminosity at specific wavelengths can also be predicted, albeit with a larger scatter due to the variation with viewing angle. These formulae for dust attenuation appear to be valid for both isolated and interacting galaxies, are consistent with earlier studies, and would be suitable for inclusion in theoretical models, e.g. semi-analytic models of galaxy formation.

We discuss how measurements of the absorption of gamma-rays from GeV to TeV energies via pair production on the extragalactic background light (EBL) can probe important issues in galaxy formation. Semi-analytic models (SAMs) of galaxy formation, based on the flat LCDM hierarchical structure formation scenario, are used to make predictions of the EBL from 0.1 to 1000 microns. SAMs incorporate simplified physical treatments of the key processes of galaxy formation -- including gravitational collapse and merging of dark matter halos, gas cooling and dissipation, star formation, supernova feedback and metal production. We will summarize SAM successes and failures in accounting for observations at low and high redshift. New ground- and space-based gamma ray telescopes will help to determine the EBL, and also help to explain its origin by constraining some of the most uncertain features of galaxy formation theory, including the stellar initial mass function, the history of star formation, and the reprocessing of light by dust. On a separate topic concerning gamma ray cosmology, we discuss a new theoretical insight into the distribution of dark matter at the center of the Milky Way, and its implications concerning the high energy gamma rays observed from the Galactic center.

A new Monte-Carlo radiative-transfer code, Sunrise, is used to study the effects of dust in N-body/hydrodynamic simulations of interacting galaxies. Dust has a profound effect on the appearance of the simulated galaxies. At peak luminosities, about 90% of the bolometric luminosity is absorbed, and the dust obscuration scales with luminosity in such a way that the brightness at UV/visual wavelengths remains roughly constant. A general relationship between the fraction of energy absorbed and the ratio of bolometric luminosity to baryonic mass is found. Comparing to observations, the simulations are found to follow a relation similar to the observed IRX-Beta relation found by Meurer et al (1999) when similar luminosity objects are considered. The highest-luminosity simulated galaxies depart from this relation and occupy the region where local (U)LIRGs are found. This agreement is contingent on the presence of Milky-Way-like dust, while SMC-like dust results in far too red a UV continuum slope to match observations. The simulations are used to study the performance of star-formation indicators in the presence of dust. The far-infrared luminosity is found to be reliable. In contrast, the H-alpha and far-UV luminosity suffer severely from dust attenuation, and dust corrections can only partially remedy the situation.

We investigate the effect of dark energy on the density profiles of dark matter haloes with a suite of cosmological N-body simulations and use our results to test analytic models. We consider constant equation of state models, and allow both w>-1 and w<-1. Using five simulations with w ranging from -1.5 to -0.5, and with more than ~1600 well-resolved haloes each, we show that the halo concentration model of Bullock et al. (2001) accurately predicts the median concentrations of haloes over the range of w, halo masses, and redshifts that we are capable of probing. We find that the Bullock et al. (2001) model works best when halo masses and concentrations are defined relative to an outer radius set by a cosmology-dependent virial overdensity. For a fixed power spectrum normalization and fixed-mass haloes, larger values of w lead to higher concentrations and higher halo central densities, both because collapse occurs earlier and because haloes have higher virial densities. While precise predictions of halo densities are quite sensitive to various uncertainties, we make broad comparisons to galaxy rotation curve data. At fixed power spectrum normalization (fixed sigma_8), w>-1 quintessence models seem to exacerbate the central density problem relative to the standard w=-1 model. Meanwhile w<-1 models help to reduce the apparent discrepancy. We confirm that the Jenkins et al. (2001) halo mass function provides an excellent approximation to the abundance of haloes in our simulations and extend its region of validity to include models with w<-1.

The good agreement between large-scale observations and the predictions of the now-standard $\Lambda$CDM theory gives us hope that this will become a lasting foundation for cosmology. After briefly reviewing the current status of the key cosmological parameters, I summarize several of the main areas of possible disagreement between theory and observation: big bang nucleosynthesis, galaxy centers, dark matter substructure, and angular momentum, updating my earlier reviews [1]. The issues in all of these are sufficiently complicated that it is not yet clear how serious they are, but there is at least some reason to think that the problems will be resolved through a deeper understanding of the complicated astrophysics involved in such processes as gas cooling, star formation, and feedback from supernovae and AGN. Meanwhile, searches for dark matter are dramatically improving in sensitivity, and gamma rays from dark matter annihilation at the galactic center may have been detected by H.E.S.S.

We describe a quasi-equilibrium profile of dark matter particles in the inner parsec of the Galaxy, rho ~ r^-3/2. This "mini-cusp" profile is caused by scattering with the dense stellar cluster around the supermassive black hole in Sgr A* and is independent of the initial conditions. The implications for detection of gamma rays from WIMP dark matter annihilation in the Galactic center are a mild enhancement of the flux, and a characteristic central feature in the angular distribution which could be detectable by high-resolution Atmospheric Cherenkov Telescopes.

We use high-resolution dissipationless simulations of the concordance flat LCDM model to make predictions for the galaxy--mass (GM) correlations and compare them to the recent SDSS weak lensing measurements.We use a simple observationally motivated scheme to assign luminosities and colors to the halos.This allows us to closely match the selection criteria used to define observational samples.The simulations reproduce the observed GM correlation function and its observed dependencies on luminosity and color.The luminosity dependence of the correlation function is primarily determined by the changing relative contribution of central and satellite galaxies at different luminosities. The color dependence of the GM correlations reflects the difference in the typical environments of blue and red galaxies. We also find agreement between the predicted and observed cross-bias, b_x=b/r,at all probed scales.The GM correlation coefficient, r, is close to unity on scales >1/h Mpc.The cross bias is thus expected to measure the actual bias of galaxy clustering on these scales.The aperture mass-to-light ratio is independent of galaxy color.The aperture mass scales approximately linearly with luminosity at L_r>10^{10}h^{-2} Lsun, while at lower luminosities the scaling is shallower: L_r^{0.5}. We show that most of the luminous galaxies (M_r<-21) are near the centers of their halos and their GM correlation function at r<100/h kpc can therefore be interpreted as the average dark matter density profile of these galaxies. We find that for galaxies in a given narrow luminosity range, there is a broad and possibly non-gaussian distribution of halo virial masses. Therefore, the average relation between mass and luminosity derived from the weak lensing analyses should be interpreted with caution.

We analyze the halo occupation distribution (HOD), the probability for a halo of mass M to host a number of subhalos N, and two-point correlation function of galaxy-size dark matter halos using high-resolution dissipationless simulations of the concordance flat LCDM model. The halo samples include both the host halos and the subhalos, distinct gravitationally-bound halos within the virialized regions of larger host systems. We find that the first moment of the HOD, <N>(M), has a complicated shape consisting of a step, a shoulder, and a power law high-mass tail. The HOD can be described by a Poisson statistics at high halo masses but becomes sub-Poisson for <N><4. We show that the HOD can be understood as a combination of the probability for a halo of mass M to host a central galaxy and the probability to host a given number Ns of satellite galaxies. The former can be approximated by a step-like function, while the latter can be well approximated by a Poisson distribution, fully specified by its first moment <Ns>(M). We find that <Ns>~M^b with b~1 for a wide range of number densities, redshifts, and different power spectrum normalizations. This formulation provides a simple but accurate model for the halo occupation distribution found in simulations. At z=0, the two-point correlation function (CF) of galactic halos can be well fit by a power law down to ~100/h kpc with an amplitude and slope similar to those of observed galaxies. At redshifts z>~1, we find significant departures from the power-law shape of the CF at small scales. If the deviations are as strong as indicated by our results, the assumption of the single power law often used in observational analyses of high-redshift clustering is likely to bias the estimates of the correlation length and slope of the correlation function.

Sand, Treu, & Ellis (2002) have measured the central density profile of cluster MS2137-23 with gravitational lensing and velocity dispersion and removed the stellar contribution with a reasonable M/L. The resulting dark matter distribution within r<50 kpc was fitted by a density cusp of r^{-beta} with beta=0.35. This stands in an apparent contradiction to the CDM prediction of beta~1, and the disagreement worsens if adiabatic compression of the dark matter by the infalling baryons is considered. Following El-Zant, Shlosman & Hoffman (2001), we argue that dynamical friction acting on galaxies moving within the dark matter background counters the effect of adiabatic compression by transfering the orbital energy of galaxies to the dark matter, thus heating up and softening the central density cusp. Using N-body simulations of massive solid clumps moving in clusters we show that indeed the inner dark matter distribution flattens (with beta approx 0.35 for a cluster like MS2137-23) when the galaxies spiral inward. We find as a robust result that while the dark matter distribution becomes core-like, the overall mass distribution preserves its cuspy nature, in agreement with X-ray and lensing observations of clusters.

We present two new non-parametric methods for quantifying galaxy morphology: the relative distribution of the galaxy pixel flux values (the Gini coefficient or G) and the second-order moment of the brightest 20% of the galaxy's flux (M20). We test the robustness of G and M20 to decreasing signal-to-noise and spatial resolution, and find that both measures are reliable to within 10% at average signal-to-noise per pixel greater than 3 and resolutions better than 1000 pc and 500 pc, respectively. We have measured G and M20, as well as concentration (C), asymmetry (A), and clumpiness (S) in the rest-frame near-ultraviolet/optical wavelengths for 150 bright local "normal" Hubble type galaxies (E-Sd) galaxies and 104 0.05 < z < 0.25 ultra-luminous infrared galaxies (ULIRGs).We find that most local galaxies follow a tight sequence in G-M20-C, where early-types have high G and C and low M20 and late-type spirals have lower G and C and higher M20. The majority of ULIRGs lie above the normal galaxy G-M20 sequence, due to their high G and M20 values. Their high Gini coefficients arise from very bright nuclei, while the high second-order moments are produced by multiple nuclei and bright tidal tails. All of these features are signatures of recent and on-going mergers and interactions. We also find that in combination with A and S, G is more effective than C at distinguishing ULIRGs from the "normal" Hubble-types. Finally, we measure the morphologies of 45 1.7 < z < 3.8 galaxies from HST NICMOS observations of the Hubble Deep Field North. We find that many of the z $\sim$ 2 galaxies possess G and A higher than expected from degraded images of local elliptical and spiral galaxies, and have morphologies more like low-redshift single nucleus ULIRGs.

We report on the merger-induced generation of a shock-heated gas wind and formation of a remnant gas halo in simulations of colliding disk galaxies. The simulations use cosmologically motivated initial conditions and include the effects of radiative cooling, star formation, stellar feedback and the non-adiabatic heating of gas. The non-adiabatic heating, i.e. shocks, generated in the final merger forces gas out of the central region of the merger remnant and into the dark-matter halo. We demonstrate that the amount of heating depends on the size of the progenitor disk galaxy as well as the initial orbit the galaxies are placed on. Based upon these dependencies, we motivate a possible recipe for including this effect in semi-analytic models of galaxy formation.

This is a summary of the discussion in Panel I of IAU Symposium 220 on Dark Matter in Galaxies. The panel topics and panellists were as follows -- Cusps Theory: Julio Navarro; Cusps Observations: Erwin de Blok, Rob Swaters; Substructure Observations: Mario Mateo; Substructure Theory: James Taylor, Shude Mao. In addition to the talks by the panellists and the discussion among them, there was a good deal of discussion with the audience. Despite the differences of opinion expressed I think a consensus emerged, and I try to summarize it here.

Cold Dark Matter with a large cosmological constant (LambdaCDM) appears to fit large scale structure observations well. Of the possible small scale problems, the Central Cusps and Too Many Satellites problems now appear to be at least partly solved, so Angular Momentum has become the most serious remaining CDM problem. There are actually at least two different angular momentum problems: A. Too much transfer of angular momentum to the dark halo to make big disks, and B. Wrong distribution of spedific angular momentum to make spiral galaxies, if the baryonic material has the same angular momentum distribution as the dark matter. The angular momentum of dark matter halos, and presumably that of the galaxies they host, appears to arise largely from the orbital angular momentum of the satellites that they accrete. Since the dark and baryonic matter behave very differently in such accretion events, it is possible that the resulting angular momentum distribution of the baryons is different from that of the dark matter, as required to make the sort of galactic disks that are observed. The latest hydrodynamical simulations give some grounds for hope on this score, but much higher resolution simulations are needed.

We investigate a model for the high-ionization state gas associated with observed damped Lyman-alpha systems, based on a semi-analytic model of galaxy formation set within the paradigm of hierarchical structure formation. In our model, the hot gas in halos and sub-halos gives rise to CIV absorption, while the low-ionization state gas is associated with the cold gas in galaxies. The model matches the distribution of CIV column densities and leads naturally to kinematic properties that are in good agreement with the data. We examine the contribution of both hot and cold gas to sub-damped systems and suggest that the properties of these systems can be used as an important test of the model. We expect that sub-DLA systems will generally be composed of a single gas disk and thus predict that they should have markedly different kinematics than the damped systems. Finally, we find that hot halo gas produces less than one third of Lyman limit systems at redshift three. We model the contribution of mini-halos (halos with virial velocities < 35 km/s) to Lyman limit systems and find that they may contain as much gas as is observed in these systems. However, if we adopt realistic models of the gas density distribution we find that these systems are not a significant source of Lyman limit absorption. Instead we suggest that uncollapsed gas outside of virialized halos is responsible for most of the Lyman limit systems at high redshift.

We propose a new explanation for the origin of angular momentum in galaxies and their dark halos, in which the halos obtain their spin through the cumulative acquisition of angular momentum from satellite accretion. In our model, the build-up of angular momentum is a random walk process associated with the mass assembly history of the halo's major progenitor. We assume no correlation between the angular momenta of accreted objects. Using the extended Press-Schechter approximation, we calculate the growth of mass, angular momentum, and spin parameter $\lambda$ for many halos. Our random walk model reproduces the key features of the angular momentum of halos found in N-body simulations: a lognormal distribution in $\lambda$ with an average of $<\lambda> \approx 0.04$, independent of mass and redshift. The evolution of the spin parameter in individual halos in this model is quite different from the steady increase with time of angular momentum in the tidal torque picture. We find both in N-body simulations and in our random walk model that the value of $\lambda$ changes significantly with time for a halo's major progenitor. It typically has a sharp increase due to major mergers, and a steady decline during periods of gradual accretion of small satellites. The model predicts that on average the $\lambda$ of halos which had major mergers after redshift $z=2$ should be substantially larger than the $\lambda$ of those which did not. Perhaps surprisingly, this suggests that halos that host late-forming elliptical galaxies should rotate faster than halos of spiral galaxies.

Cold Dark Matter (CDM) has become the standard modern theory of cosmological structure formation. Its predictions appear to be in good agreement with data on large scales, and it naturally accounts for many properties of galaxies. But despite its many successes, there has been concern about CDM on small scales because of the possible contradiction between the linearly rising rotation curves observed in some dark-matter-dominated galaxies vs. the $1/r$ density cusps at the centers of simulated CDM halos. Other CDM issues on small scales include the very large number of small satellite halos in simulations, far more than the number of small galaxies observed locally, and problems concerning the angular momentum of the baryons in dark matter halos. The latest data and simulations have lessened, although not entirely resolved, these concerns. Meanwhile, the main alternatives to CDM that have been considered to solve these problems, self-interacting dark matter (SIDM) and warm dark matter (WDM), have been found to have serious drawbacks.

I discuss how measurements of the absorption of $\gamma$-rays from GeV to TeV energies via pair production on the extragalactic background light (EBL) can probe important issues in galaxy formation. We use semi-analytic models (SAMs) of galaxy formation, based on the flat LCDM hierarchical structure formation scenario with $\Omega_m=0.3$ and Hubble parameter $h=0.65$, to obtain predictions of the EBL from 0.1 to 1000$\mu$m. SAMs incorporate simplified physical treatments of the key processes of galaxy formation - including gravitational collapse and merging of dark matter halos, gas cooling and dissipation, star formation, supernova feedback and metal production -- and have been shown to reproduce key observations at low and high redshift. We investigate the consequences of variations in input assumptions such as the stellar initial mass function (IMF) and the efficiency of converting cold gas into stars. We also discuss recent attempts to determine the emitted spectrum of high energy gamma rays from blazars such as Mrk 501 from the observed X-rays using the synchrotron self-Compton model, and note that our favorite SAM EBL plus the observed spectrum of Mrk 501 do not imply unphysical upturns in the high energy emitted spectrum - thus undermining recent claims of a crisis with drastic possible consequences such as breaking of Lorentz invariance. We conclude that observational studies of the absorption of $\gamma$-rays will help to determine the EBL, and also help to explain its origin by constraining some of the most uncertain features of galaxy formation theory, including the IMF, the history of star formation, and the reprocessing of light by dust.

Data on galaxies at high redshift, identified by the Lyman-break photometric technique, can teach us about how galaxies form and evolve. The stellar masses and other properties of such Lyman break galaxies (LBGs) depend sensitively on the details of star formation. In this paper we consider three different star formation prescriptions, and use semi-analytic methods applied to the now-standard $\Lambda$CDM theory of hierarchical structure formation to show how these assumptions about star formation affect the predicted masses of the stars in these galaxies and the masses of the dark matter halos that host them. We find that, within the rather large uncertainties, recent estimates of the stellar masses of LBGs from multi-color photometry are consistent with the predictions of all three models. However, the estimated stellar masses are more consistent with the predictions of two of the models in which star formation is accelerated at high redshifts $z\gsim3$, and of these models the one in which many of the LBGs are merger-driven starbursts is also more consistent with indications that many high redshift galaxies are gas rich. The clustering properties of LBGs have put some constraints on the masses of their host halos, but due to similarities in the halo occupation of the three models we consider and degeneracies between model parameters, current constraints are not yet sufficient to distinguish between realistic models.

This article summarizes the possible roles of neutrinos in cosmology, from the first three minutes onward. The fact that primordial neutrinos are about as numerous as the photons of the cosmological background radiation means that neutrino (or ``hot'') dark matter can have a significant impact upon how matter is distributed in the universe. Over the past two decades, the liklihood of this has risen and fallen as more and more data has become available from laboratory experiments and the latest telescopes. Hot dark matter models were studied in the late 1970s and early 1980s, and there was considerable interest in cold plus hot dark matter models for a few years in the mid-1990s. We now know from Super-Kamiokande and Sudbury Neutrino Observatory data that neutrinos have mass, and that the total neutrino contribution to the cosmological density must be at least as large as $\Omega_\nu \sim 0.001$, a little less than that of the visible stars in the universe. But the upper limit on the hot dark matter cosmological density is $\Omega_\nu < 0.1$; correspondingly, the sum of the masses of all three active neutrino species is not more than about 5 eV.

Cold Dark Matter (CDM) has become the standard modern theory of cosmological structure formation. Its predictions appear to be in good agreement with data on large scales, and it naturally accounts for many properties of galaxies. But despite its many successes, there has been concern about CDM on small scales because of the possible contradiction between the linearly rising rotation curves observed in some dark-matter-dominated galaxies vs. the $1/r$ density cusps at the centers of simulated CDM halos. Other CDM issues on small scales include the very large number of small satellite halos in simulations, far more than the number of small galaxies observed locally, and problems concerning the angular momentum of the baryons in dark matter halos. The latest data and simulations have lessened, although not entirely resolved, these concerns. Meanwhile, the main alternatives to CDM that have been considered to solve these problems, self-interacting dark matter (SIDM) and warm dark matter (WDM), have been found to have serious drawbacks.

(abridged) We study the relation between the density profiles of dark matter halos and their mass assembly histories, using a statistical sample of halos in a high-resolution N-body simulation of the LCDM cosmology. For each halo at z=0, we identify its merger-history tree, and determine concentration parameters c_vir for all progenitors, thus providing a structural merger tree for each halo. We fit the mass accretion histories by a universal function with one parameter, the formation epoch a_c, defined when the log mass accretion rate dlogM/dloga falls below a critical value S. We find that late forming galaxies tend to be less concentrated, such that c_vir ``observed'' at any epoch a_o is strongly correlated with a_c via c_vir=c_1*a_o/a_c. Scatter about this relation is mostly due to measurement errors in c_v and a_c, implying that the actual spread in c_vir for halos of a given mass can be mostly attributed to scatter in a_c. We demonstrate that this relation can also be used to predict the mass and redshift dependence of c_v, and the scatter about the median c_vir(M,z), using accretion histories derived from the Extended Press-Schechter (EPS) formalism, after adjusting for a constant offset between the formation times as predicted by EPS and as measured in the simulations;this new ingredient can thus be easily incorporated into semi-analytic models of galaxy formation. The correlation found between halo concentration and mass accretion rate suggests a physical interpretation: for high mass infall rates the central density is related to the background density; when the mass infall rate slows, the central density stays approximately constant and the halo concentration just grows as R_vir. The tight correlation demonstrated here provides an essential new ingredient for galaxy formation modeling.

We study the relation between the density profiles of dark matter halos and their mass assembly histories, using a statistical sample of halos in a high-resolution N-body simulation of the LCDM cosmology. For each halo at z=0, we identify its merger-history tree, and determine concentration parameters c_vir for all progenitors, thus providing a structural merger tree for each halo. We fit the mass accretion histories by a universal function with one parameter, the formation epoch a_c, defined when the log mass accretion rate dlog M/doga falls below a critical value. We find that late forming galaxies tend to be less concentrated, such that c_vir ``observed'' at any epoch a_obs is strongly correlated with a_c via c_vir=c_1 a_obs/a_c. Scatter about this relation is mostly due to measurement errors in c_vir and a_c, implying that the actual spread in c_vir for halos of a given mass can be mostly attributed to scatter in a_c. Because of the direct connection between halo concentration and velocity rotation curves, and because of probable connections between halo mass assembly history and star formation history, the tight correlation between these properties provides an essential new ingredient for galaxy formation modeling.

We present the three-dimensional power spectrum of galaxy clustering using measurements of the angular correlation function from the second incremental release of the Two-Micron All Sky Survey (2MASS). The angular positions of galaxies out to a limiting magnitude of K_s <= 14.0 (508,054 galaxies) determine the correlation function. We consider a variety of estimators of the catalog's angular correlation function, and we find the Landy and Szalay (LS) estimator to be the most reliable choice. Using the LS estimator, we find that the angular correlation function takes the form of a power law with a break. The angular correlation function is then used in combination with the selection function to determine the power spectrum by inverting Limber's equation. We find good agreement with previous measurements of the power spectrum for 0.04 < k < 0.1 h Mpc^-1, and evidence for a peak in the power spectrum around k = 0.03 h Mpc^-1. For k > 0.1 h Mpc^-1 we find that the power spectrum is biased higher than other surveys, which can be understood as a color-based selection effect.

We study the effects of mass and force resolution on the density profiles of galaxy-size Cold Dark Matter (CDM) halos in a flat, low-density cosmological model with vacuum energy. We show thatalthough increasing the mass and force resolution allows us to probe deeper into the inner halo regions, it does not lead to steeper inner density profiles. Instead, the halo profiles converge at scales larger than four times the formal resolution or the radius containing more than 200 particles, whichever is larger. In the simulations presented in this paper, we are able to probe density profile of a relaxed isolated galaxy-size halo at scales r=(0.005-1)r_vir. We find that the density distribution can be well approximated by the profile suggested by Moore etal (1998): rho= x^{-1.5}(1+x^{1.5})^-1, where x=r/r_s and r_s is the characteristic radius. The analytical profile proposed by Navarro et al. (1996) rho= x^{-1}(1+x)^-2, also provides a good fit, with the same relative errors of about 10% for radii larger than 1% of the virial radius. Both analytical profiles fit equally well because for high-concentration galaxy-size halos the differences between these profiles become significant only at scales well below 0.01r_vir. We also find that halos of similar mass may have somewhat different parameters (characteristic radius, maximum rotation velocity, etc.) and shapes of their density profiles. We associate this scatter in properties with differences in halo merger histories and the amount of substructure present in the analyzed halos.

We use N-body simulations combined with semi-analytic models to compute the clustering properties of modeled galaxies at z~3, and confront these predictions with the clustering properties of the observed population of Lyman-break galaxies (LBGs). Several scenarios for the nature of LBGs are explored, which may be broadly categorized into models in which high-redshift star formation is driven by collisional starbursts and those in which quiescent star formation dominates. For each model, we make predictions for the LBG overdensity distribution, the variance of counts-in-cells, the correlation length, and close pair statistics. Models which assume a one-to-one relationship between massive dark-matter halos and galaxies are disfavored by close pair statistics, as are models in which colliding halos are associated with galaxies in a simplified way. However, when modeling of gas consumption and star formation is included using a semi-analytic treatment, the quiescent and collisional starburst models predict similar clustering properties and none of these models can be ruled out based on the available clustering data. None of the ``realistic'' models predict a strong dependence of clustering amplitude on the luminosity threshold of the sample, in apparent conflict with some observational results.

We discuss how measurements of the absorption of $\gamma$-rays from GeV to TeV energies via pair production on the extragalactic background light (EBL) can probe important issues in galaxy formation. We use semi-analytic models (SAMs) of galaxy formation, set within the hierarchical structure formation scenario, to obtain predictions of the EBL for 0.1-1000$\mu$m. SAMs incorporate simplified physical treatments of the key processes of galaxy formation --- including gravitational collapse and merging of dark matter halos, gas cooling and dissipation, star formation, supernova feedback and metal production --- and have been shown to reproduce key observations at low and high redshift. Here we also introduce improved modelling of the spectral energy distributions in the mid-to-far-IR arising from emission by dust grains. Assuming a flat \lcdm cosmology with $\Omega_m=0.3$ and Hubble parameter $h=0.65$, we investigate the consequences of variations in input assumptions such as the stellar initial mass function (IMF) and the efficiency of converting cold gas into stars. We conclude that observational studies of the absorption of $\gamma$-rays with energies from 10s of Gev to 10s of TeV will help to determine the EBL, and also help to explain its origin by constraining some of the most uncertain features of galaxy formation theory, including the IMF, the history of star formation, and the reprocessing of light by dust.

[Abridged] We study the angular-momentum profiles of a statistical sample of halos drawn from a high-resolution N-body simulation of the LCDM cosmology. We find that the cumulative mass distribution of specific angular momentum, j, in a halo of mass Mv is well fit by a universal function, M(<j) = Mv \mu j/(j_0+j). This profile is defined by one shape parameter (\mu or j_0) in addition to the global spin parameter \lambda. It follows a power-law over most of the mass, and flattens at large j, with the flattening more pronounced for small values of \mu. Compared to a uniform sphere in solid-body rotation, most halos have a higher fraction of their mass in the low- and high-j tails of the distribution. The spatial distribution of angular momentum in halos tends to be cylindrical and is well-aligned within each halo for ~80% of the halos. We investigate two ideas for the origin of this profile. The first is based on a revised version of linear tidal-torque theory combined with extended Press-Schechter mass accretion, and the second focuses on j transport in minor mergers. Finally, we briefly explore implications of the M(<j) profile on the formation of galactic disks assuming that j is conserved during an adiabatic baryonic infall. The implied gas density profile deviates from an exponential disk, with a higher density at small radii and a tail extending to large radii. The steep central density profiles may imply disk scale lengths that are smaller than observed. This is reminiscent of the "angular-momentum problem" seen in hydrodynamic simulations, even though we have assumed perfect j conservation. A possible solution is to associate the central excesses with bulge components and the outer regions with extended gaseous disks.

We summarize recent developments in the study of the origin of halo spin profiles and preliminary implications on disk formation. The specific angular-momentum distributions within halos in N-body simulations match a universal shape, M(<j) \propto j/(j_0+j). It is characterized by a power law over most of the mass, and one shape parameter in addition to the spin parameter \lambda. The angular momentum tends to be aligned throughout the halo and of cylindrical symmetry. Even if angular momentum is conserved during baryonic infall, the resultant disk density profile is predicted to deviate from exponential, with a denser core and an extended tail. A slightly corrected version of the scaling relation due to linear tidal-torque theory is used to explain the origin of a typical power-law profile in shells, j(M) \propto M^s with s \gsim 1. While linear theory crudely predicts the amplitudes of halo spins, it is not a good predictor of their directions. Independently, mergers of halos are found to produce a similar profile due to j transfer from the orbit to the product halo via dynamical friction and tidal stripping. The halo spin is correlated with having a recent major merger, though this correlation is weakened by mass loss. These two effects are due to a correlation between the spins of neighboring halos and their orbit, leading to prograde mergers.

On the basis of a new convergence study of high-resolution N-body simulations, my colleagues and I now agree that the Navarro, Frenk, & White (1996) density profile $\rho_{NFW}(r) \propto r^{-1} (r+r_s)^{-2}$ is a good representation of typical dark matter halos of galactic mass. Comparing simulations of the same halo with numbers of particles ranging from $\sim10^3$ to $\sim10^6$, we have also shown that $r_s$, the radius where the log-slope is -2, can be determined accurately for halos with as few as $\sim10^3$ particles. Based on a study of thousands of halos at many redshifts in an Adaptive Refinement Tree (ART) simulation of a cosmological volume in a $\Lambda$CDM cosmology, we have found that the concentration $\cvir \equiv \Rvir/\rs$ has a log-normal distribution, with $1\sigma$ $\Delta (\log \cvir) = 0.18$ at a given mass, corresponding to a scatter in maximum rotation velocities of $\Delta \Vmax/\Vmax = 0.12$. The average concentration declines with redshift at fixed mass as $\cvir(z) \propto (1+z)^{-1}$. This may have important implications for galaxy rotation curves. Finally, we have found that the velocity function determined from galaxy luminosity functions plus luminosity-velocity relations agrees with the predictions from our $\Lambda$CDM simulations. But we also note that the very limited evolution with redshift of the velocity function predicted by $\Lambda$CDM conflicts with the data that is becoming available on the number density of bright galaxies unless there is significant evolution of the luminosity-velocity relation at $z>1$.

[Abridged] We analyze the exchange of dark matter between halos, subhalos, and their environments in a high-resolution cosmological N-body simulation of a Lambda CDM cosmology. At each analyzed redshift z we divide the dark matter particles into 4 components: (i) isolated galactic halos, (ii) subhalos, (iii) the diffuse medium of group and cluster halos, and (iv) the background outside of virialized halos. We follow the time evolution of the mass distribution and flows between these components and provide fitting functions for the exchange rates. We use our derived exchange rates to gauge the importance of metal redistribution in the universe due solely to gravity-induced interactions. The diffuse metallicity in clusters is predicted to be ~40% that in isolated galaxies (~55% of groups) at z=0, and should be lower only slightly by z=1, consistent with observations. The metallicity of the diffuse media in poor groups is expected to be lower by a factor of 5 by z~2, in agreement with the observed metallicity of damped Ly$\alpha$ systems. The metallicity of the background IGM is predicted to be (1-3)x10^{-4} that of z=0 clusters, also consistent with observations. The agreement of predicted and observed trends indicates that gravitational interaction alone may play an important role in metal enrichment of the intra-cluster and intergalactic media.

[Abridged] The interaction rates of dark-matter halos and subhalos, are computed using high-resolution cosmological N-body simulations of the Lambda-CDM model. Although the number fraction of subhalos of mass $>2x10^{11}\hsolmass$ is only ~10%, we find that the interaction rate of such subhalos is relatively high because they reside in high density environments. At low redshift, the subhalo collisions dominate the total collision rate, and even at z=3 they are involved in more than 30% of all collisions. About 40% of the "major" mergers are between subhalos. Therefore subhalo interactions must be incorporated in models of structure formation. We find that the collision rate between halos in physical density units, is $\propto (1+z)^\delta $, with $\delta = 3-4$, in agreement with earlier simulations and most observational data. We test previous analytic estimates of the interaction rates of subhalos based on statistical models, which could be very inaccurate because of the small number of subhalos and the variation of conditions within small host halos. We find that, while such statistical estimates may severely overestimate the rate within hosts $< 10^{13} \hsolmass$, typical of high redshifts, they are valid for larger hosts regardless of the number of subhalos in them. We find the Makino & Hut (1997) estimate of the subhalo merger rate to be valid for hosts $ > 6x10^{11}\hsolmass$ at all redshifts. The collision rate between subhalos and the central object of their host halo is approximated relatively well using the timescale for dynamical friction in circular orbits. This approximation fails in ~40% of the cases, partly because of deviations from the assumption of circular orbits and partly because of the invalidity of the assumption that the subhalo mass is negligible.

We study dark-matter halo density profiles in a high-resolution N-body simulation of an LCDM cosmology. Our statistical sample contains ~5000 haloes in the range 10^{11}-10^{14} M_sun. The profiles are parameterized by an NFW form with two parameters, an inner radius r_s and a virial radius r_v, and we define the halo concentration c_v = r_v/r_s. First, we find that, for a given halo mass, the redshift dependence of the median concentration is c_v ~ 1/(1+z), corresponding to a roughly constant r_s with redshift. We present an improved analytic treatment of halo formation that fits the measured relations between halo parameters and their redshift dependence. The implications are that high-redshift galaxies are predicted to be more extended and dimmer than expected before. Second, we find that the scatter in log(c_v) is ~0.18, corresponding to a scatter in maximum rotation velocities of dV/V ~ 0.12. We discuss implications for modelling the Tully-Fisher relation, which has a smaller reported intrinsic scatter. Third, haloes in dense environments tend to be more concentrated than isolated haloes. These results suggest that c_v is an essential parameter for the theory of galaxy modelling, and we briefly discuss implications for the universality of the Tully-Fisher relation, the formation of low surface brightness galaxies, and the origin of the Hubble sequence.

The cosmological parameters that I emphasize are the age of the universe $t_0$, the Hubble parameter $H_0 \equiv 100 h$ km s$^{-1}$ Mpc$^{-1}$, the average matter density $\Omega_m$, the baryonic matter density $\Omega_b$, the neutrino density $\Omega_\nu$, and the cosmological constant $\Omega_\Lambda$. The evidence currently favors $t_0 \approx 13$ Gyr, $h \approx 0.65$, $\Omega_m \approx 0.4\pm0.1$, $\Omega_b = 0.02h^{-2}$, $0.001 < \Omega_\nu < 0.1$, and $\Omega_\Lambda \approx 0.7$.

Cosmological dark matter in the form of neutrinos with masses of up to a few electron volts is known as hot dark matter. After an historical review of the subject, this article considers constraints on hot dark matter from current data on neutrino oscillations and on cosmology. The atmospheric neutrino oscillation data imply a lower limit on the HDM contribution to the cosmological density $\Omega_\nu \gsim 0.001$. The possible improvement of low-$\Omega_m$ flat ($\Lambda$CDM) cosmological models with the addition of light neutrinos appears to be rather limited, but$\Lambda$CDM models with $\Omega_\nu \lsim 0.1$ may be consistent with presently available data. Data expected soon may permit detection of such a hot dark matter contribution, or alternatively provide stronger upper limits on $\Omega_\nu$ and neutrino masses.

This paper has been withdrawn because the models on which it was based have undergone significant changes and improvements. A new paper with the same title, based on the improved models, is accepted for publication in MNRAS and is available at astro-ph/0006364. We have withdrawn this paper to avoid confusion.

Using semi-analytic models of galaxy formation set within the Cold Dark Matter (CDM) merging hierarchy, we investigate several scenarios for the nature of the high-redshift ($z \ga 2$) Lyman-break galaxies (LBGs). We consider a ``collisional starburst'' model in which bursts of star formation are triggered by galaxy-galaxy mergers, and find that a significant fraction of LBGs are predicted to be starbursts. This model reproduces the observed comoving number density of bright LBGs as a function of redshift and the observed luminosity function at $z\sim 3$ and $z\sim 4$, with a reasonable amount of dust extinction. Model galaxies at $z \sim 3$ have star formation rates, half-light radii, I-K colours, and internal velocity dispersions that are in good agreement with the data. Global quantities such as the star formation rate density and cold gas and metal content of the Universe as a function of redshift also agree well. Two ``quiescent'' models without starbursts are also investigated. In one, the star formation efficiency in galaxies remains constant with redshift, while in the other, it scales inversely with disc dynamical time, and thus increases rapidly with redshift. The first quiescent model is strongly ruled out as it does not produce enough high redshift galaxies once realistic dust extinction is accounted for. The second quiescent model fits marginally, but underproduces cold gas and very bright galaxies at high redshift. A general conclusion is that star formation at high redshift must be more efficient than locally. The collisional starburst model appears to accomplish this naturally without violating other observational constraints.

We study the distribution functions of mass and circular velocity for dark matter halos in N-body simulations of the $\Lambda$CDM cosmology, addressing redshift and environmental dependence. The dynamical range enables us to resolve subhalos and distinguish them from "distinct" halos. The mass function is compared to analytic models, and is used to derive the more observationally relevant circular velocity function. The distribution functions in the velocity range 100--500 km/s are well fit by a power-law with two parameters, slope and amplitude. We present the parameter dependence on redshift and provide useful fitting formulae. The amplitudes of the mass functions decrease with z, but, contrary to naive expectation, the comoving density of halos of a fixed velocity ~200 km/s actually increases out to z=5. This is because high-z halos are denser, so a fixed velocity corresponds to a smaller mass. The slope of the velocity function at z=0 is as steep as ~ -4, and the mass and velocity functions of distinct halos steepen with increasing z, while the functions of subhalos do not steepen with $z$, and become even flatter at z>2. A simple observable prediction is that the slope of the velocity function of isolated galaxies is steeper than that of galaxies in groups by as much as unity, reflecting the density biasing of high-velocity halos. We confirm that the Press-Schechter approximation typically overestimates the halo mass function by a factor of ~2, while modified approximations provide improved predictions.

We present a method for constraining the evolution of the galaxy luminosity-velocity (LV) relation in hierarchical scenarios of structure formation. The comoving number density of dark-matter halos with circular velocity of 200 km/s is predicted in favored CDM cosmologies to be nearly constant over the redshift range 0<z<5. Any observed evolution in the density of bright galaxies implies in turn a corresponding evolution in the LV relation. We consider several possible forms of evolution for the zero-point of the LV relation and predict the corresponding evolution in galaxy number density. The Hubble Deep Field suggests a large deficit of bright (M_V < -19) galaxies at 1.4 < z < 2. If taken at face value, this implies a dimming of the LV zero-point by roughly 2 magnitudes. Deep, wide-field, near-IR selected surveys will provide more secure measurements to compare with our predictions.

We investigate the properties of damped Lyman-alpha systems (DLAS) in semi-analytic models, focusing on whether the models can reproduce the kinematic properties of low-ionization metal lines described by Prochaska & Wolfe (1997b, 1998). We explore a variety of approaches for modelling the radial distribution of the cold neutral gas associated with the galaxies in our models, and find that our results are very sensitive to this ingredient. If we use an approach based on Fall & Efstathiou (1980), in which the sizes of the discs are determined by conservation of angular momentum, we find that the majority of the DLAS correspond to a single galactic disc. These models generically fail to reproduce the observed distribution of velocity widths. In alternative models in which the gas discs are considerably more extended, a significant fraction of DLAS arise from lines of sight intersecting multiple gas discs in a common halo. These models produce kinematics that fit the observational data, and also seem to agree well with the results of recent hydrodynamical simulations. Thus we conclude that Cold Dark Matter based models of galaxy formation can be reconciled with the kinematic data, but only at the expense of the standard assumption that DLAS are produced by rotationally supported gas discs whose sizes are determined by conservation of angular momentum. We suggest that the distribution of cold gas at high redshift may be dominated by another process, such as tidal streaming due to mergers.

Previously we have shown (Maller et al, 1998) that the kinematics of Damped Lyman Alpha Systems (DLAS) as measured by Prochaska and Wolfe (1998) can be reproduced in a multiple disk model (MDM) if the gaseous disks are of sufficient radial extent. Here we discuss this model's predictions for the relationship between DLAS and Lyman break galaxies (LBGs), which we here take to be objects at z~3 brighter than R=25.5. We expect that future observations of the correlations between DLAS and LBGs will provide a new data set able to discriminate between different theoretical models of the DLAS. Djorgovski (1997) has already detected a few optical counterparts and more studies are underway.

Previously we have shown that it is possible to account for the kinematic properties of damped Lyman alpha systems (DLAS) in the context of semi-analytic models. In these models, hierarchical structure formation is approximated by constructing a merger tree for each dark matter halo. A natural consequence is that every virialized halo may contain not only a central galaxy, but also a number of satellite galaxies as determined by its merging history. Thus the kinematics of the DLAS arise from the combined effects of the internal rotation of gas disks and the motions between gas disks within a common halo. Here we investigate the sensitivity of this model to some of the assumptions made previously, including the modeling of satellite dynamics, the scale height of the gas, and the cosmology.

The cosmological parameters that I will emphasize are the Hubble parameter $H_0 \equiv 100 h$ km s$^{-1}$ Mpc$^{-1}$, the age of the universe $t_0$, the average matter density $\Omega_m$, the baryonic matter density $\Omega_b$, the neutrino density $\Omega_\nu$, and the cosmological constant $\Omega_\Lambda$. The evidence currently favors $t_0 \approx 13$ Gyr, $h \approx 0.65$, $\Omega_m \approx 0.3$, $\Omega_\Lambda \approx 0.7$.

The degeneracy between the disk and the dark matter contribution to galaxy rotation curves remains an important uncertainty in our understanding of disk galaxies. Here we discuss a new method for breaking this degeneracy using gravitational lensing by spiral galaxies, and apply this method to the spiral lens B1600+434 as an example. The combined image and lens photometry constraints allow models for B1600+434 with either a nearly singular dark matter halo, or a halo with a sizable core. A maximum disk model is ruled out with high confidence. Further information, such as the circular velocity of this galaxy, will help break the degeneracies. Future studies of spiral galaxy lenses will be able to determine the relative contribution of disk, bulge, and halo to the mass in the inner parts of galaxies.

We investigate the clustering properties of high-redshift galaxies within three competing scenarios for assigning luminous galaxies to dark matter halos from N-body simulations: a one galaxy per massive halo model, a quiescent star formation model, and a collisional starburst model. We compare these models to observations of Lyman-Break galaxies at z~3$ With current data and the simple statistic used here, one cannot rule out any of these models, but we see potential for finding distinguishing features using statistics that are sensitive to the tails of the distribution, and statistics based on the number of multiple galaxies per halo, which we explore in an ongoing study.

Using high resolution N-body simulations with hydrodynamics and star formation, we investigate interactions and the resulting starbursts in galaxies with properties typical of $z\sim 3$. We apply spectral population models to produce mock-HST images, and discuss the observed magnitude, color, and morphological appearance of our simulated galaxies in both the rest-UV and rest-visual bands.

This paper represents a step toward constraining galaxy formation models via TeV gamm a ray observations. We use semi-analytic models of galaxy formation to predict a spectral distribution for the intergalactic infrared photon field, which in turn yields information about the absorption of TeV gamma rays from extra-galactic sources. By making predictions for integral flux observations at >200 GeV for several known EGRE T sources, we directly compare our models with current observational upper limits obtained by Whipple. In addition, our predictions may offer a guide to the observing programs for the current population of TeV gamma ray observatories.

The frequency with which background galaxies appear as long arcs as a result of gravitational lensing by foreground clusters of galaxies has recently been found to be a very sensitive probe of cosmological models by Bartelmann et al. (1998). They have found that such arcs would be expected far less frequently than observed (by an order of magnitude) in the currently favored model for the universe, with a large cosmological constant $\Omega_\Lambda \sim 0.7$. Here we analyze whether including the effect of cluster galaxies on the likelihood of clusters to generate long-arc images of background galaxies can change the statistics. Taking into account a variety of constraints on the properties of cluster galaxies, we find that there are not enough sufficiently massive galaxies in a cluster for them to significantly enhance the cross section of clusters to generate long arcs. We find that cluster galaxies typically enhance the cross section by only $\lesssim 15%$.

We present a galaxy circular velocity function, Psi(log v), derived from existing luminosity functions and luminosity-velocity relations. Such a velocity function is desirable for several reasons. First, it enables an objective comparison of luminosity functions obtained in different bands and for different galaxy morphologies, with a statistical correction for dust extinction. In addition, the velocity function simplifies comparison of observations with predictions from high-resolution cosmological N-body simulations. We derive velocity functions from five different data sets and find rough agreement among them, but about a factor of 2 variation in amplitude. These velocity functions are then compared with N-body simulations of a LCDM model (corrected for baryonic infall) in order to demonstrate both the utility and current limitations of this approach. The number density of dark matter halos and the slope of the velocity function near v_*, the circular velocity corresponding to an ~L_* spiral galaxy, are found to be comparable to that of observed galaxies. The primary sources of uncertainty in construction of Psi(log v) from observations and N-body simulations are discussed and explanations are suggected to account for these discrepancies.

Using semi-analytic models of galaxy formation, we investigate galaxy properties such as the Tully-Fisher relation, the B and K-band luminosity functions, cold gas contents, sizes, metallicities, and colours, and compare our results with observations of local galaxies. We investigate several different recipes for star formation and supernova feedback, including choices that are similar to the treatment in Kauffmann, White & Guiderdoni (1993) and Cole et al. (1994) as well as some new recipes. We obtain good agreement with all of the key local observations mentioned above. In particular, in our best models, we simultaneously produce good agreement with both the observed B and K-band luminosity functions and the I-band Tully-Fisher relation. Improved cooling and supernova feedback modelling, inclusion of dust extinction, and an improved Press-Schechter model all contribute to this success. We present results for several variants of the CDM family of cosmologies, and find that models with values of $\Omega_0 \simeq 0.3$--0.5 give the best agreement with observations.

Lyman break galaxies (LBGs) at $z \sim 3$ exhibit number densities and clustering similar to local $L_*$ galaxies; however, their star formation rates (SFRs) are much higher. We explore the scenario in which LBGs are starburst galaxies triggered by collisions, and thus provide an explanation for these key properties. The relative paucity of starburst galaxies at low redshift can be attributed to a much higher collision rate in the past. We use high-resolution cosmological $N$-body simulations and a hierarchical halo finder to estimate the galaxy collision rate as a function of time in the popular $\Lambda$CDM cosmological model. We find that bright collisional starbursts are frequent enough to account for most of the high-$z$ ($2.5-4.5$) LBGs. Moreover, many of the objects are of relatively small mass, but they cluster about large-mass halos. They therefore exhibit strong clustering, similar to that observed and stronger than that of the relevant massive halos.

We present here the extragalactic background light (EBL) predicted by semi-analytic models of galaxy formation, and show how measurements of the absorption of gamma rays of $\sim$ TeV energies via pair production on the EBL can probe cosmology and the formation of galaxies. Semi-analytic models permit a physical treatment of the key processes of galaxy formation -- including gravitational collapse and merging of dark matter halos, gas cooling and dissipation, star formation, supernova feedback and metal production -- and have been shown to reproduce key observations at low and high redshift. Using this approach, we investigate the consequences of variations in input assumptions such as the stellar initial mass function and the underlying cosmology. We conclude that observational studies of the absorption of $\sim 10^{-2}-10^{2}$ TeV gamma rays will help to constrain the star formation history of the universe, and the nature and extent of the extinction of starlight due to dust and reradiation of the absorbed energy at infrared wavelengths.

We review our analysis of the clustering properties of ``Lyman-break'' galaxies (LBGs) at redshift z~3, previously discussed in Wechsler et al (1998). We examine the likelihood of spikes found by Steidel et al (1998) in the redshift distribution of LBGs, within a suite of models for the evolution of structure in the Universe. Using high-resolution dissipationless N-body simulations, we analyze deep pencil-beam surveys from these models in the same way that they are actually observed, identifying LBGs with the most massive dark matter halos. We find that all the models (with SCDM as a marginal exception) have a substantial probability of producing spikes similar to those observed, because the massive halos are much more clumped than the underlying matter -- i.e., they are biased. Therefore, the likelihood of such a spike is not a good discriminator among these models. The LBG correlation functions are less steep than galaxies today (gamma~1.4), but show similar or slightly longer correlation lengths. We have extened this analysis and include a preliminary comparison to the new data presented in Adelberger et al (1998). We also discuss work in progress, in which we use semi-analytic models to identify Lyman-break galaxies within dark-matter halos.

The velocity profiles of weak metal absorption lines can be used to observationally probe the kinematic state of gas in damped Lyman-alpha systems. Prochaska and Wolfe have argued that the flat distribution of velocity widths combined with the asymmetric line profiles indicate that the DLAS are disks with large rotation velocities. An alternative explanation has been proposed by Haehnelt, Steinmetz, and Rauch, in which the observed large velocity widths and asymmetric profiles can be produced by lines of sight passing through two or more clumps each having relatively small internal velocity dispersions. We investigate the plausibility of this scenario in the context of semi-analytic models based on hierarchical merging trees and including simple treatments of gas dynamics, star formation, supernova feedback, and chemical evolution. We find that all the observed properties of the metal-line systems including the distribution of velocities and the asymmetric profiles, can be reproduced by lines of sight passing through sub-clumps that are bound within larger virialized dark matter halos. In order to produce enough multiple hits, we find that the cold gas must be considerably more extended than the optical radius of the proto-galaxies, perhaps even beyond the tidal radius of the sub-halo. This could occur due to tidal stripping or supernova-driven outflows.

Gamma rays of ~TeV energies from distant sources suffer attenuation due to pair production off of ~1 micron EBL photons. We may exploit this process in order to indirectly measure the EBL and constrain models of galaxy formation. Here, using semi-analytic models of galaxy formation, we examine how gamma ray absorption may be used as an indirect probe of the stellar initial mass function (IMF), although there is a degeneracy with dust modeling. We point out that with the new generation of gamma ray telescopes including STACEE, MAGIC, HESS, VERITAS, and Milagro, we should soon possess a wealth of new data and a new method for probing the nature of the IMF.

The recent atmospheric neutrino data from Super-Kamiokande provide strong evidence of neutrino oscillations and therefore of non-zero neutrino mass. These data imply a lower limit on the hot dark matter (i.e., light neutrino) contribution to the cosmological density $\Omega_\nu >~ 0.001$ --- almost as much as that of all the stars in the universe --- and permit higher $\Omega_\nu$. The ``standard'' COBE-normalized critical-matter-density (i.e., $\Omega_m=1$) Cold Dark Matter (CDM) model has too much power on small scales. But adding to CDM neutrinos with mass of about 5 eV, corresponding to $\Omega_\nu \approx 0.2$, results in a much improved fit to data on the nearby galaxy and cluster distribution. Indeed, the resulting Cold + Hot Dark Matter (CHDM) cosmological model is arguably the most successful $\Omega_m=1$ model for structure formation. However, other recent data has begun to make the case for $\Omega_m <~ 0.6$ fairly convincing. In light of all this new data, we reconsider whether cosmology still provides evidence favoring neutrino mass of a few eV in flat models with cosmological constant $\Omega_\Lambda = 1 - \Omega_m$. We find that the possible improvement of the low-$\Omega_m$ flat (LCDM) cosmological models with the addition of light neutrinos appears to be rather limited.

This article concerns the formation and structure of dark matter halos, including (1) their radial density profiles, (2) their abundance, and (3) their merger rates. The last topic may be relevant to the nature of the small, bright, high-redshift galaxies discovered by the Lyman break technique. (1) Study of a statistical sample of galaxy-mass dark halos in high-resolution Adaptive Refinement Tree simulations shows that they have a central density profile p(r) \propto 1/r^g with g \approx 0.2, in agreement with data on dark-matter-dominated disk galaxies. We present recent, higher resolution results on this. (2) Another important new result is that the Press-Schechter approximation predicts about twice as many galaxy-mass halos at z=0 as are present in large dissipationless N-body simulations; more generally, PS overpredicts the abundance of M ~< 0.1M_* halos at all redshifts. (3) Finally, we discuss the assembly of these halos, in particular the merger rate of (sub-)halos at high redshift and the distribution of the starbursts that these mergers are likely to trigger. If most of the Lyman-break galaxies are such starbursts, this perhaps resolves the apparent paradox that these galaxies appear to cluster like massive halos (~10^{12} M_\odot), while their relatively low linewidths and their spectral energy distributions suggest that they have relatively low mass (few X 10^{10} M_\odot) and young ages (few X 10^8 yr). It also predicts much more star formation at high redshift in CDM-type hierarchical models for structure formation than if only quiescent star formation is included.

Observations now probe the star formation history of the Universe back to a redshift of $z\sim5$. We investigate whether the predictions of semi-analytic models of galaxy formation based on hierarchical Cold Dark Matter (CDM) type models are in agreement with these direct observations and also with the ``fossil'' evidence contained in constraints on the ages of present day early-type galaxies. Previous models predicted that the star formation rate density falls off rather steeply at $z\ga 2$, and correspondingly that the majority of the stars in the Universe formed at relatively low redshift. We investigate the effect of including a bursting mode of star formation, assuming that galaxy-galaxy mergers trigger starbursts and using the merger rate that arises naturally in the CDM merging hierarchy. The resulting starbursts substantially increase the global star formation rate at high redshift, leading to predictions that are in good agreement with the star formation rate density at $z\sim3$ obtained from sub-millimeter observations (SCUBA) and optical/UV estimates after correction for dust extinction. The mass of stars formed at $z \ge 3$ is correspondingly in better agreement with the fossil evidence. We also investigate complementary global quantities such as the mass of cold gas and the average metallicity of cold gas as a function of redshift, and the integrated extra-galactic background light. We argue that these observations, taken together, provide strong constraints on the star formation history of the Universe, and that hierarchical models of the CDM type are in reasonable agreement with these observations when starbursts are included.

We review the basic arguments for the likelihood of non-Gaussian density perturbations in inflation models with primordial black hole (PBH) production. We discuss our derived distributions of field fluctuations and their implications, specifically commenting on the fine-tuning problem. We also discuss how the derived distributions may be affected when linked to metric perturbations. While linking the metric perturbations to field fluctuations in a nonlinear way may be important for determining exact probability distributions, the correct mapping is not self-evident. The calculation of P. Ivanov, which yields skew positive distribution, is based on an ansatz for the behavior of the nonlinear metric perturbation. We note that the ``natural'' generalization of the gauge-invariant formalism favored by Bond and Salopek yields an effective linear link between the distribution of field fluctuations and metric perturbations during inflation.

Using semi-analytic models of galaxy formation, we investigate the properties of $z\sim3$ galaxies and compare them with the observed population of Lyman-break galaxies (LBGs). In addition to the usual quiescent mode of star formation, we introduce a physical model for starbursts triggered by galaxy-galaxy interactions. We find that with the merger rate that arises naturally in the CDM-based merging hierarchy, a significant fraction of bright galaxies identified at high redshift ($z \ga 2$) are likely to be low-mass, bursting satellite galaxies. The abundance of LBGs as a function of redshift and the luminosity function of LBGs both appear to be in better agreement with the data when the starburst mode is included, especially when effects of dust extinction are considered. The objects that we identify as LBGs have observable properties including low velocity dispersions that are in good agreement with the available data. In this ``Bursting Satellite'' scenario, quiescent star formation at $z\ga2$ is relatively inefficient and most of the observed LBGs are starbursts triggered by satellite mergers within massive halos. In high-resolution N-body simulations, we find that the most massive dark matter halos cluster at redshift $z\sim 3$ much as the LBGs are observed to do. This is true for both the $\Omega=1$ CHDM model and low-$\Omega$ \LCDM and OCDM models, all of which have fluctuation power spectra $P(k)$ consistent with the distribution of low-redshift galaxies. The Bursting Satellite scenario can resolve the apparent paradox of LBGs that cluster like massive dark matter halos but have narrow linewidths and small stellar masses.

We use the rotation curves of a sample of dark matter dominated dwarf and low-surface brightness (LSB) late-type galaxies to study their radial mass distributions. We find that the shape of the rotation curves is remarkably similar for all (both dwarf and LSB) galaxies in the sample, suggesting a self-similar density distribution of their dark matter (DM) halos. This shape can be reproduced well by a density profile with a shallow central cusp (rho(r)\propto 1/r^gamma, gamma\approx 0.2-0.4), corresponding to a steeply rising velocity curve (v(r)\propto r^g, g\approx 0.9-0.8). We further show that the observed shape of the rotation curves is well matched by the density profiles of DM halos formed in very high resolution simulations of the CDM, LCDM, and CHDM models of structure formation. This is surprising in light of several previous studies which suggested that the structure of simulated DM halos is inconsistent with the dynamics of dwarf galaxies. We discuss possible explanations for this discrepancy. We show that although the mass distribution in the hierarchically formed halos is on average consistent with the shape of rotation curves of dark matter dominated galaxies, the scatter of the individual profiles around the average is substantial and should not be neglected in comparisons with the data. Finally, we show that the DM halos in our hierarchical simulations and the real galaxies in our sample exhibit a similar decrease in their characteristic densities with increasing characteristic radial scales, and increase in their maximum rotation velocities with increasing radii at which their maximum velocity occurs. (Abridged)

We present two matched sets of five simulations each, covering five presently favored simple modifications to the standard cold dark matter (CDM) scenario. One simulation suite, with a linear box size of 75 Mpc/h, is designed for high resolution and good statistics on the group/poor cluster scale, and the other, with a box size of 300 Mpc/h, is designed for good rich cluster statistics. All runs had 57 million cold particles, and models with massive neutrinos had an additional 113 million hot particles. We consider separately models with massive neutrinos, tilt, curvature, and a nonzero cosmological constant in addition to the standard CDM model. We find that our tilted Omega+Omega_Lambda=1 (TLCDM) model produces too much small-scale power by a factor of ~3, and our open Lambda=0 (OCDM) model also exceeds observed small-scale power by a factor of 2. In addition, we take advantage of the large dynamic range in detectable halo masses our simulations allow to check the shape of the Press-Schechter approximation. We find good fits at cluster masses for delta_c=1.27--1.35 for a Gaussian filter and delta_c=1.57--1.73 for a tophat filter. However, Press-Schechter overpredicts the number density of halos compared to the simulations in the high resolution suite by a weakly cosmology-dependent factor of 1.5--2 at galaxy and group masses, which cannot be fixed by adjusting delta_c within reasonable bounds. An appendix generalizes the spherical collapse model to any isotropic cosmology.

We address the high peaks found by Steidel et al (1997) in the redshift distribution of ``Lyman-break'' objects (LBOs) at redshift z~3. The highest spike represents a relative overdensity of 2.6 in the distribution of LBOs in pixels of comoving size ~10Mpc/h. We examine the likelihood of such a spike in the redshift distribution within a suite of models for the evolution of structure in the Universe, including models with Omega=1 (SCDM and CHDM) and with Omega=0.4-0.5 (LCDM and OCDM). Using high-resolution dissipationless N-body simulations, we analyze deep pencil-beam surveys from these models in the same way that they are actually observed, identifying LBOs with the most massive dark matter halos. We find that all the models (with SCDM as a marginal exception) have a substantial probability of producing spikes similar to those observed, because the massive halos are much more clumped than the underlying matter -- i.e., they are biased. Therefore, the likelihood of such a spike is not a good discriminator among these models. We find in these models that the mean biasing parameter b of LBOs with respect to dark matter varies within a range b ~2-5 on a scale of ~10Mpc/h. We also compute the two-body correlation functions of LBOs predicted in these models. The LBO correlation functions are less steep than galaxies today (gamma ~1.4), but show similar or slightly longer correlation lengths.

The combination of measurements of the local abundance of rich clusters of galaxies and its evolution to higher redshift offers the possibility of a direct measurement of $\Omega_0$ with little contribution from other cosmological parameters. We investigate the significance of recent claims that this evolution indicates that $\Omega_0$ must be small. The most recent cluster velocity dispersion function from a compilation including the ESO Northern Abell Cluster Survey (ENACS) results in a significantly higher normalization for models, corresponding to $\sigma_8\approx 0.6$ for $\Omega_0=1$, compared to the Eke, Cole, & Frenk result of $\sigma_8=0.52\pm 0.04$. Using the ENACS data for a $z=0$ calibration results in strong evolution in the abundance of clusters, and we find that the velocity dispersion function is consistent with $\Omega_0=1$. The results are dependent upon the choice and analysis of low-redshift and high-redshift data, so at present, the data is not good enough to determine $\Omega_0$ unambiguously.

This chapter aims to present an introduction to current research on the nature of the cosmological dark matter and the origin of galaxies and large scale structure within the standard theoretical framework: gravitational collapse of fluctuations as the origin of structure in the expanding universe. General relativistic cosmology is summarized, and the data on the basic cosmological parameters ($t_o$ and $H_0 \equiv 100 h \kmsmpc$, $\Omega_0$, $\Omega_\Lambda$ and $\Omega_b$) are reviewed. Various particle physics candidates for hot, warm, and cold dark matter are briefly reviewed, together with current constraints and experiments that could detect or eliminate them. Also included is a very brief summary of the theory of cosmic defects, and a somewhat more extended exposition of the idea of cosmological inflation with a summary of some current models of inflation. The remainder is a discussion of observational constraints on cosmological model building, emphasizing models in which most of the dark matter is cold and the primordial fluctuations are the sort predicted by inflation. It is argued that the simplest models that have a hope of working are Cold Dark Matter with a cosmological constant ($\Lambda$CDM) if the Hubble parameter is high ($h \gsim 0.7$), and Cold + Hot Dark Matter (CHDM) if the Hubble parameter and age permit an $\Omega=1$ cosmology, as seems plausible in light of the data from the Hipparcos astrometric satellite. The most attractive variants of these models and the critical tests for each are discussed.

Spheroidal components of spiral galaxies have been considered the only dynamically important component in gravitational lensing studies thus far. Here we point out that including the disk component can have a significant effect, depending on the disk inclination, on a variety of lensing properties that are relevant to present studies and future surveys. As an example, we look at the multiple image system B1600+434, recently identified as being lensed by a spiral galaxy. We find that including the disk component one can understand the fairly large image separation as being due to the inclination of a typical spiral, rather than the presence of a very massive halo. The fairly low magnification ratio can also be readily understood if the disk is included. We also discuss how such lensed systems might allow one to constrain parameters of spiral galaxies such as a disk-to-halo mass ratio, and disk mass scale length. Another example we consider is the quasar multiple-lensing cross section, which we find can increase many-fold at high inclination for a typical spiral. Finally, we discuss the changes in the gravitational lensing effects on damped Lyman alpha systems (DLAS) when disk lensing is included.

Reconstruction of the linear power spectrum from observational data provides a way to compare cosmological models to a large amount of data, as Peacock & Dodds (1994, 1996) have shown. By applying the appropriate corrections to the observational power spectrum it is possible to recover the underlying linear power spectrum for any cosmological model. Using extensive N-body simulations we demonstrate that the method is applicable to a wide range of cosmological models. However, we find that the recovery of the linear power spectrum from observations following PD94 is misleading because the corrections are model- dependent. When we apply the proper corrections for a given model to the observational power spectrum, we find that no model in our test group recovers the linear power spectrum well for the bias suggested by PD94 between Abell, Radio, Optical, and IRAS catalogs 4.5:1.9:1.3:1, with b_IRAS=1. When we allow b_IRAS to vary we find that: (i)CHDM models give very good fits to observations if optically-selected galaxies are slightly biased b_Opt=1.1 (ii) Most LCDM models give worse but acceptable fits if blue galaxies are considerably antibiased: 0.6<b_Opt<0.9 and fail if optical galaxies are biased. (iii)There is a universal shape of the recovered linear power spectrum of all LCDM models over their entire range of explored wavenumbers,0.01<k<0.6h\Mpc. Recovered spectra of CDM and CHDM models are nearly the same as that of LCDM in the region 0.01<k<0.2h/Mpc but diverge from this spectrum at higher k.

We discuss variants of Cold Dark Matter (CDM) dominated cosmological models that give good agreement with a range of observations. We consider models with hot dark matter, tilt, $\Omega < 1$, or a cosmological constant. We also discuss the sensitivity of the results to other parameters, such as the Hubble parameter and the baryon fraction. We obtain constraints by combining the COBE data, cluster abundances, abundance of damped Lyman-$\alpha$ systems at $z\sim3$, the small-angle Cosmic Microwave Background anisotropy, and the small-scale non-linear power spectrum. We present non-linear power spectra from a new suite of N-body simulations for the ``best-bet'' models from each category.

We explore the role of non-Gaussian fluctuations in primordial black hole (PBH) formation and show that the standard Gaussian assumption, used in all PBH formation papers to date, is not justified. Since large spikes in power are usually associated with flat regions of the inflaton potential, quantum fluctuations become more important in the field dynamics, leading to mode-mode coupling and non-Gaussian statistics. Moreover, PBH production requires several sigma (rare) fluctuations in order to prevent premature matter dominance of the universe, so we are necessarily concerned with distribution tails, where any intrinsic skewness will be especially important. We quantify this argument by using the stochastic slow-roll equation and a relatively simple analytic method to obtain the final distribution of fluctuations. We work out several examples with toy models that produce PBH's, and show that the naive Gaussian assumption can result in errors of many orders of magnitude. For models with spikes in power, our calculations give sharp cut-offs in the probability of large positive fluctuations, meaning that Gaussian distributions would vastly over-produce PBH's. The standard results that link inflation-produced power spectra and PBH number densities must then be reconsidered, since they rely quite heavily on the Gaussian assumption. We point out that since the probability distributions depend on the nature of the potential, it is impossible to obtain results for general models. However, calculating the distribution of fluctuations for any specific model seems to be relatively straightforward, at least in the single inflaton case.

(Shortened Abstract). We analyze a redshift sample of Abell/ACO clusters and compare them with numerical simulations based on the truncated Zel'dovich approximation (TZA), for a list of eleven dark matter (DM) models. For each model we run several realizations, on which we estimate cosmic variance effects. We analyse correlation statistics, the probability density function, and supercluster properties from percolation analysis. As a general result, we find that the distribution of galaxy clusters provides a constraint only on the shape of the power spectrum, but not on its amplitude: a shape parameter 0.18 < \Gamma < 0.25 and an effective spectral index at 20Mpc/h in the range [-1.1,-0.9] are required by the Abell/ACO data. In order to obtain complementary constraints on the spectrum amplitude, we consider the cluster abundance as estimated using the Press--Schechter approach, whose reliability is explicitly tested against N--body simulations. We conclude that, of the cosmological models considered here, the only viable models are either Cold+Hot DM ones with \Omega_\nu = [0.2-0.3], better if shared between two massive neutrinos, and flat low-density CDM models with \Omega_0 = [0.3-0.5].

We apply several statistical estimators to high-resolution N-body simulations of two currently viable cosmological models: a mixed dark matter model, having $\Omega_\nu=0.2$ contributed by two massive neutrinos (C+2\nuDM), and a Cold Dark Matter model with Cosmological Constant (\LambdaCDM) with $\Omega_0=0.3$ and h=0.7. Our aim is to compare simulated galaxy samples with the Perseus-Pisces redshift survey (PPS). We consider the n-point correlation functions (n=2-4), the N-count probability functions P_N, including the void probability function P_0, and the underdensity probability function U_\epsilon (where \epsilon fixes the underdensity threshold in percentage of the average). We find that P_0 (for which PPS and CfA2 data agree) and P_1 distinguish efficiently between the models, while U_\epsilon is only marginally discriminatory. On the contrary, the reduced skewness and kurtosis are, respectively, S_3\simeq 2.2 and S_4\simeq 6-7 in all cases, quite independent of the scale, in agreement with hierarchical scaling predictions and estimates based on redshift surveys. Among our results, we emphasize the remarkable agreement between PPS data and C+2\nuDM in all the tests performed. In contrast, the above \LambdaCDM model has serious difficulties in reproducing observational data if galaxies and matter overdensities are related in a simple way.

I have been asked to make the case here that CHDM is the best theory of cosmic structure formation, and indeed I believe that it is the best of all those I have considered if the cosmological matter density is near critical and if the expansion rate is not too large (i.e. $h \equiv H_0/(100 \kmsMpc) \lsim 0.6$). I discuss CHDM together with its chief competitor among CDM variants, low-$\Omega_0$ CDM with a cosmological constant ($\Lambda$CDM). $\Lambda$CDM with $\Omega_0 \sim 0.3$ has the possible virtue of allowing a higher expansion rate H_0 for a given cosmic age t_0, but the defect of predicting too much fluctuation power on small scales. Also, except for the $H_0-t_0$ problem, there is not a shred of evidence in favor of a nonzero cosmological constant, only increasingly stringent observational upper bounds on it. CHDM has less power on small scales, in good agreement with data, although it remains to be seen whether it predicts early enough galaxy formation to be compatible with the latest high-redshift data. I also briefly compare CHDM to other CDM variants such as Warm Dark Matter (WDM) and tilted CDM. CHDM has the advantage among $\Omega=1$ CDM-type models of requiring little or no tilt, which appears to be an advantage in fitting recent small-angle CMB anisotropy data. The presence of a hot component that clusters less than cold dark matter lowers the effective $\Omega_0$ that would be measured on small scales, which appears to be in accord with observations, and it may also avoid the discrepancy between the high central density of dark matter halos from CDM simulations compared to evidence from rotation curves of dwarf spiral galaxies.

A new class of geometric statistics for analyzing galaxy catalogs is presented. Filament statistics quantify filamentarity and planarity in large scale structure in a manner consistent with catalog visualizations. These statistics are based on sequences of spatial links which follow local high-density structures. From these link sequences we compute the discrete curvature, planarity, and torsion. Filament statistics are applied to CDM and CHDM ($\Omega_\nu = 0.3$) simulations of Klypin \etal (1996), the CfA1-like mock redshift catalogs of Nolthenius, Klypin and Primack (1994, 1996), and the CfA1 catalog. We also apply the moment-based shape statistics developed by Babul \& Starkman (1992), Luo \& Vishniac (1995), and Robinson \& Albrecht (1996) to these same catalogs, and compare their robustness and discriminatory power versus filament statistics. For 100 Mpc periodic simulation boxes ($H_0 = 50$ km s$^{-1}$ Mpc$^{-1}$), we find discrimination of $\sim 4\sigma$ (where $\sigma$ represents resampling errors) between CHDM and CDM for selected filament statistics and shape statistics, including variations in the galaxy identification scheme. Comparing the CfA1 data versus the models does not yield a conclusively favored model; no model is excluded at more than a $\sim 2\sigma$ level for any statistic, not including cosmic variance which could further degrade the discriminatory power. We find that CfA1 discriminates between models poorly mainly due to its sparseness and small number of galaxies, not due to redshift distortion, magnitude limiting, or geometrical effects. We anticipate that the proliferation of large redshift surveys and simulations will enable the statistics presented here to provide robust discrimination between large-scale structure in various cosmological models.

Here we discuss what are perhaps the two most popular variants of CDM that might agree with the data: \lcdm\ and CHDM. While the predictions of COBE-normalized \lcdm\ and CHDM both agree well with the available data on scales of $\sim 10$ to $100 \hMpc$, each has potential virtues and defects. \lcdm\ with $\Omega_0 \sim 0.3$ has the possible virtue of allowing a higher expansion rate $H_0$ for a given cosmic age $t_0$, but the defect of predicting too much fluctuation power on small scales. CHDM has less power on small scales, so its predictions appear to be in good agreement with data on the galaxy distribution, but it remains to be seen whether it predicts early enough galaxy formation to be compatible with the latest high-redshift data. Also, two very recent observational results favor high cosmic density, and thus favor $\Omega=1$ models such as CHDM over \lcdm\ --- (1) the positive deceleration parameter $q_0>0$ measured using high-redshift Type Ia supernovae, and (2) the low primordial deuterium/hydrogen ratio measured in two different quasar absorption spectra. We try to identify ``best'' variants of both \lcdm\ and CHDM, and discuss critical observational tests for both models.

The cosmological parameters that I discuss are the Hubble parameter $H_0 \equiv 100 h$ km s$^{-1}$ Mpc$^{-1}$, the age of the universe $t_0$, the average density $\Omega_0$, and the cosmological constant $\Lambda$. To focus the discussion, I concentrate on the the value of $\Omega_0$ in currently popular models in which most of the dark matter is cold, especially Cold + Hot Dark Matter (CHDM) and flat ($\Omega_0 + \Omega_\Lambda=1$) low-$\Omega$ CDM with a Cosmological Constant ($\Lambda$CDM). The evidence would favor small $\Omega_0 \approx 0.3$ if (1) the Hubble parameter actually has the high value $h \approx 0.75$ favored by many observers, and $t_0 \geq 13$ Gy; or (2) the baryonic/total mass ratio in clusters of galaxies is actually $\sim 15$\%, about 3 times larger than expected for standard BBN in an $\Omega=1$ universe, $\Omega_b \approx 0.0125 h^{-2}$, despite the recent measurement by Tytler of $D/H=2.4\times 10^{-5}$ in two high-redshift Lyman limit systems, implying $\Omega_b\approx 0.024 h^{-2}$. The evidence would favor $\Omega=1$ if (1) the POTENT analysis of galaxy peculiar velocity data is right, in particular regarding outflows from voids or the inability to obtain the present-epoch non- Gaussian density distribution from Gaussian initial fluctuations in a low- $\Omega$ universe; or (2) the preliminary LSND report indicating neutrino mass $\gsim 2.4$ eV is right, since that would be too much hot dark matter to allow significant structure formation in a low-$\Omega_0$ $\Lambda$CDM model. Statistics on gravitational lensing of quasars provide an upper limit on $\Lambda$, and the preliminary results on the deceleration parameter $q_0=\Omega_0/2-\Omega_\Lambda$ on very large scales from high-redshift Type Ia supernovae suggest that $\Omega_0 \sim 1$ and $\Omega_\Lambda$ is small.

The velocity dispersion of galaxies on small scales ($r\sim1h^{-1}$ Mpc), $\sigma_{12}(r)$, can be estimated from the anisotropy of the galaxy-galaxy correlation function in redshift space. We apply this technique to ``mock-catalogs'' extracted from N-body simulations of several different variants of Cold Dark Matter dominated cosmological models to obtain results which may be consistently compared to similar results from observations. We find a large variation in the value of $\sigma_{12}(1 h^{-1} Mpc)$ in different regions of the same simulation. We conclude that this statistic should not be considered to conclusively rule out any of the cosmological models we have studied. We attempt to make the statistic more robust by removing clusters from the simulations using an automated cluster-removing routine, but this appears to reduce the discriminatory power of the statistic. However, studying $\sigma_{12}$ as clusters with different internal velocity dispersions are removed leads to interesting information about the amount of power on cluster and subcluster scales. We also compute the pairwise velocity dispersion directly and compare this to the values obtained using the Davis-Peebles method, and find that the agreement is fairly good. We evaluate the models used for the mean streaming velocity and the pairwise peculiar velocity distribution in the original Davis-Peebles method by comparing the models with the results from the simulations.

The velocity dispersion of galaxies on scales of $r\sim1h^{-1}$ Mpc, $\sigma_{12}(r)$, may be estimated from the anisotropy of the galaxy-galaxy correlation function in redshift space. We present a reanalysis of the CfA1 survey, correct an error in the original analysis of Davis and Peebles (1983), and find that $\sigma_{12}(r)$ is extremely sensitive to the details of how corrections for infall into the Virgo cluster are applied. We conclude that a robust value of $\sigma_{12}$ cannot be obtained from this survey. We also discuss results from other redshift surveys, including the effect of removing clusters.

Many multiply--imaged quasars have been found over the years, but none so far with image separation in excess of $8\arcsec$. The absence of such large splittings has been used as a test of cosmological models: the standard Cold Dark Matter model has been excluded on the basis that it predicts far too many large--separation double images. These studies assume that the lensing structure has the mass profile of a singular isothermal sphere. However, such large splittings would be produced by very massive systems such as clusters of galaxies, for which other gravitational lensing data suggest less singular mass profiles. Here we analyze two cases of mass profiles for lenses: an isothermal sphere with a finite core radius (density $\rho \propto (r^2+r_{core}^2)^{-1})$, and a Hernquist profile ($\rho \propto r^{-1}(r+a)^{-3}$). We find that small core radii $r_{core} \sim 30 h^{-1}$ kpc, as suggested by the cluster data, or large $a \gsim 300 h^{-1}$ kpc, as needed for compatibility with gravitational distortion data, would reduce the number of large--angle splittings by an order of magnitude or more. Thus, it appears that these tests are sensitive both to the cosmological model (number density of lenses) and to the inner lens structure, which is unlikely to depend sensitively on the cosmology, making it difficult to test the cosmological models by large--separation quasar lensing until we reliably know the structure of the lenses themselves.

Cosmological models with a positive cosmological constant and $\Omega_0<1$ have a number of attractive features. A larger Hubble constant, which can be compatible with the recent HST estimate, and a large fraction of baryon density in galaxy clusters make them current favorites. Early galaxy formation also is considered as a welcome feature of these models But early galaxy formation implies that fluctuations on few megaparsec scales spent more time in the nonlinearregime, as compared with standard Cold Dark Matter (CDM) or Cold+Hot Dark Matter models. This results in excessive clustering on small scales. We show that a typical LCDM model with $H_0=70$ km/s/Mpc, $\Omega_0=0.3$ normalized to COBE on large scales and compatible with the number-density of galaxy clusters, predicts a power spectrum of galaxy clustering in real space which is too high: {\it at least} twice larger than CfA estimates (Park \etal 1994) and 3 times larger than APM estimates (Baugh \& Efstathiou 1994) for wavenumbers $k=(0.4-1)h/{\rm Mpc}$. This conclusion holds if we assume either that galaxies trace the dark matter ($\sigma_8\approx 1.1$ for this model) or just that a region with higher density produces more galaxies than a region with lower density. The only way to reconcile the model with the observed power spectrum is to assume that regions with high dark matter density produce fewer galaxies than regions with low density. Theoretically this is possible, but it seems very unlikely: X-ray emission from groups and clusters indicates that places with a large density of dark matter produce a large number of galaxies. Since it follows that the low- $\Omega$ LCDM models are in serious trouble, we discuss which LCDM models have the best hope of surviving the confrontation with available observational data.

Theoretical models for structure formation with Gaussian initial fluctuations have been worked out in considerable detail and compared with observations on various scales. It is on nonlinear scales $\lsim 10 \ h^{-1}\ {\rm Mpc}$ that the greatest differences exist between $\Omega=1$ models that have been normalized to agree on the largest scales with the COBE data; here especially there is a need for better statistical tests which are simultaneously {\it robust}, {\it discriminatory}, and {\it interpretable}. The era at which galaxy and cluster formation occurs is also a critical test of some models. Needs for the future include faster and cleverer codes, better control of cosmic variance in simulations, better understanding of processes leading to galaxy formation, better ways of comparing observational data with models, and better access to observational and simulation data.

A new class of geometric statistics for analyzing galaxy catalogs is presented. {\it Filament statistics} quantify filamentarity and planarity in large scale structure in a manner consistent with catalog visualizations. These statistics are based on sequences of spatial links which follow local high-density structures. From these link sequences we compute the discrete curvature, planarity, and torsion. Filament statistics are applied to CDM and CHDM ($\Omega_\nu = 0.3$) simulations of Klypin \etal (1995), the CfA1-like sky catalogs of Nolthenius, Klypin and Primack (1994, 1995), and the CfA1 catalog. For 100 Mpc periodic simulation boxes ($H_0 = 50$ km s$^{-1}$ Mpc$^{-1}$), we find robust discrimination of over 4$\sigma$ (where $\sigma$ represents resampling errors) between CHDM and CDM. The {\it reduced filament statistics} show that CfA1 data is intermediate between CHDM and CDM, but more consistent with the CHDM models. Filament statistics provide robust and discriminatory shape statistics with which to test cosmological simulations of various models against present and future redshift survey data.

We present models of the extragalactic background light (EBL) based on several scenarios of galaxy formation and evolution. We have treated galaxy formation with the Press-Schecter approximation for both cold dark matter (CDM) and cold+hot dark matter (CHDM) models, representing a moderate ($z_f \sim 3$) and a late ($z_f \sim 1$) era of galaxy formation respectively. Galaxy evolution has been treated by considering a variety of stellar types, different initial mass functions and star formation histories, and with an accounting of dust absorption and emission. We find that the dominant factor influencing the EBL is the epoch of galaxy formation. A recently proposed method for observing the EBL utilizing the absorption of $\sim 0.1$ to 10 TeV gamma-rays from active galactic nuclei (AGN) is shown to be capable of discriminating between different galaxy formation epochs. The one AGN viewed in TeV light, Mrk 421, does show some evidence for a cutoff above 3 TeV; based on the EBL models presented here, we suggest that this is due to extinction in the source. The large absorption predicted at energies $>200$ GeV for sources at $z>0.5$ indicates that observations of TeV gamma-ray bursts (GRB) would constrain or eliminate models in which the GRB sources lie at cosmological distances.

I discuss the Hubble parameter $H_0$, the age of the universe $t_0$, the density $\Omega_0$, and the cosmological constant $\Lambda$. To focus the discussion, I concentrate on the issue of the value of the density $\Omega_0$ in currently popular cosmological models in which most of the dark matter is cold, especially Cold + Hot Dark Matter (CHDM) and flat low- $\Omega$ CDM with a Cosmological Constant ($\Lambda$CDM). The evidence would favor a small $\Omega_0 \approx 0.3$ if (1) the Hubble parameter actually has the high value $h \approx 0.8$ favored by many observers, and the age of the universe $t_0 \geq 13$ Gy; or (2) the baryonic/total mass ratio in clusters is actually $\sim 20$\%, about 3-4 times larger than expected for standard Big Bang Nucleosynthesis in an $\Omega=1$ universe, and standard BBN is actually right in predicting that the density of ordinary matter $\Omega_b$ lies in the range $0.009 \leq \Omega_b h^2 \leq 0.02$. The evidence would favor $\Omega=1$ if (1) the POTENT analysis of galaxy peculiar velocity data is right, in particular regarding outflows from voids or the inability to obtain the present-epoch non-Gaussian density distribution from Gaussian initial fluctuations in a low-$\Omega$ universe; or (2) the preliminary report from LSND indicating a neutrino mass $\geq 2.4$ eV is right, since that would be too much hot dark matter to allow significant structure formation in a low-$\Omega$ $\Lambda$CDM model. Statistics on gravitational lensing of quasars provide a strong upper limit on $\Lambda$. Reliable data on all of these issues is becoming available so rapidly today that a clear decision between these alternatives may be possible within the next few years.

The recent results from the Los Alamos neutrino oscillation experiment, together with assumptions of neutrino oscillation solutions for the solar and atmospheric neutrino deficit problems, may place powerful constraints on any putative scheme for neutrino masses and mixings. Assuming the validity of these experiments and assumptions, we argue that a nearly unique spectrum of neutrino masses emerges as a fit, if two additional astrophysical arguments are adopted: (1) the sum of the light neutrino masses is $\sim 5\ {\rm eV}$, as large scale structure simulations with mixed cold plus hot dark matter seem to suggest; and (2) $r$-process nucleosynthesis originates in neutrino-heated ejecta from Type II supernovae. In this fit, the masses of the neutrinos must satisfy $m_{{\nu}_e} \approx m_{{\nu}_s} \approx 2.7\ {\rm eV}$ (where ${\nu}_e$ is split from a sterile species, ${\nu}_s$, by $\sim {10}^{-6} \ {\rm eV}$) and $m_{{\nu}_{\tau}} \approx m_{{\nu}_{\mu}} \approx 1.1\ {\rm eV}$ (where these species are split by $\sim {10}^{-2} \ {\rm eV}$). We discuss alternative neutrino mass spectra that are allowed if we decline to adopt certain experiments or astrophysical models.

We present a series of four simulations of Cold Dark Matter (CDM) and Cold + Hot Dark Matter (CHDM) cosmologies. We discuss the power spectrum and correlation functions in real and redshift space, with comparisons to the CfA2 and IRAS redshift data, the pairwise velocity of galaxies, and the distribution of hot and cold particles in CHDM simulations. We confirm that CHDM with cold/hot/baryon density ratios $\Omega_c/\Omega_\nu/\Omega_b= 0.6/0.3/0.1$ is a good fit to a wide variety of present-epoch data, much better than CDM. In particular, with reasonable assumptions about identification of galaxies and biasing, we find that the power spectrum agrees rather well with both the CfA2 and IRAS power specta in both the nonlinear and linear regimes. New CHDM models (e.g., with $\Omega_\nu=0.2$ or with two massive neutrinos) predict a significantly larger rate of formation of galaxies at high redshift, which may be needed to explain some observational data. At the same time, the difference between the variants is rather small at $z=0$. The results presented in this paper are interesting for two purposes: (i) For a comparison with other classes of models (like CDM and $\Lambda$CDM) at $z=0$.(ii) As a reference point for comparison between different variants of the CHDM model.

We present models of the extragalactic background light (EBL) based on several scenarios of galaxy formation and evolution. We have treated galaxy formation with the Press-Schecter approximation for both cold dark matter (CDM) and cold+hot dark matter (CHDM) models. Galaxy evolution has been treated by considering a variety of stellar types, different initial mass functions and star formation histories, and with an accounting of dust absorption and emission. We find that the dominant factor influencing the EBL is the epoch of galaxy formation. A recently proposed method for observing the EBL utilizing the absorption of $\sim 0.1$ to 10 TeV gamma-rays from active galactic nuclei (AGN) is shown to be capable of discriminating between different galaxy formation epochs. Observations of TeV gamma-ray bursts (GRB) would constrain or eliminate models in which the GRB sources lie at cosmological distances.

Cold + Hot Dark Matter (CHDM) $\Omega=1$ cosmological models require a total neutrino mass $\sim 5$ eV. Because recent data support the $\nu_\mu \to \nu_\tau$ oscillation explanation of the cosmic ray $\nu_\mu$ deficit, which requires that $m(\nu_\mu) \approx m(\nu_\tau)$, this suggests that $m(\nu_\mu) \approx m(\nu_\tau) \approx 2.4$ eV. The linear calculations and N-body simulation reported here indicate that an $\Omega=1$ CHDM model with two 2.4 eV neutrinos (designated C$\nu^2$DM) agrees remarkably well with all available observations, but only if the Hubble parameter $h \approx 0.5$. We also show that even one 2.4 eV neutrino raises serious difficulties for low-$\Omega$ flat CDM models.

We present techniques for identifying and analyzing galaxy groups and show that they provide a powerful discriminator between cosmological models. We apply these to high-resolution PM N-body simulations of structure formation in CHDM (cold/hot/baryon fractions=.6/.3/.1, b=1.5 (COBE norm) and two CDM models; b=1.5 and b=1.0 (COBE norm). Groups are identified with the Nolthenius (1993) algorithm. We find: (1) Properties of groups are a powerful and robust discriminator between these Gaussian models whose spatial and velocity properties differ on Mpc scales. We test robustness against several methods for assigning luminosity to dark matter halos, for merging CfA1 galaxies, and for breaking up massive over-mergers (2) When allowance is made for the high large scale power P(k) present in the CfA1 data, CHDM at our $\Omnu=0.30$ produces slightly too many groups and groups too high a fraction of the galaxies, while the fraction grouped in CDM is far too low. Lowering $\Omnu$ would appear to produce close agreement with all measures, save one: for all simulations, median group sizes are x1.4-1.7 larger than equivalently selected CfA1 groups. (3) The standard group M/L method gives $\Omega\simeq 0.1$ for CHDM and CfA1, and $\Omega\simeq 0.35$ for CDM. A conspiracy of three effects produce these very low values, and we conclude that low observed $\Omega$'s need not argue for a low $\Omega$ universe. When overmergers are broken up, the median virial- to-DM mass of 3D selected groups is $\sim1$ for all simulations. $M_{DM}>10^{14}\Msun$ groups appear virialized in all simulations. Global velocity biases $b_v$ are similar to previous studies. Within 3D-selected groups, CHDM and CDM b=1.5 show a stronger bias of $b_v=0.7-0.8$, while CDM b=1.0 shows group $b_v$'s $\simeq 1$.

Damped Ly$\alpha$ systems provide possibly the most significant evidence for early structure formation, and thus a stringent constraint on the Cold + Hot Dark Matter (CHDM) cosmology. Using the numbers of halos in N-body simulations to normalize Press-Schechter (PS) estimates of the number densities of protogalaxies as a function of redshift, we find that CHDM with $\Omega_c/\Omega_\nu/\Omega_b = 0.6/0.3/0.1$ is compatible with the damped Ly$\alpha$ data only at $z<2.5$, but that it is probably incompatible with data at $z>3$. The predictions of CHDM are quite sensitive to the neutrino fraction. We find that $\Omega_c/\Omega_\nu/\Omega_b = 0.725/0.20/0.075$ is compatible with the $z>3$ data. With one massive neutrino species, this corresponds to lowering the neutrino mass from 7.0 to 4.7 eV, for $H_0=50\kmsMpc$ and $T=2.726$ K. By analysing our numerical simulations with different resolutions and box sizes as well as those of Ma \& Bertchinger (1994), we show that for the CHDM models with $\Omega_\nu$=0.2--0.3 the PS approximation should be used with Gaussian filter with $\delta_c=1.3-1.4$ if one tries to recover the total mass of a collapsed halo and to include nonlinear effects, due to waves both longer and shorter than those within the simulation box.

We use the void probability statistics to study the redshift-space galaxy distribution as described by a volume-limited subsample of the Perseus-Pisces survey. We compare the results with the same analysis realized on artificial samples, extracted from high-resolution N-body simulations by reproducing the observational biases of the real data set. Simulations are run for the Cold+HotDM model (CHDM) and for unbiased and biased (b=1.5) CDM models in a 50 Mpc/h box. We identify galaxies as residing in peaks of the evolved density field. We fragment overmerged structures into individual galaxies so as to reproduce both the correct luminosity function (after assuming M/ L values for the resulting galaxy groups) and the two-point correlation function. Our main result is that a void-probability function (VPF) from the standard CHDM model with fractions 60% cold, 30% hot, 10% barions, exceeds the observational VPF with a high confidence level. CDM models produce smaller VPF independent of the biasing parameter. We verify the robustness of this result against changing the observer position in the simulations and the galaxy identification in the evolved density field.

We discuss the effects of redshift--space distortions on the estimate of high order galaxy correlation functions. We use both the Perseus--Pisces redshift survey and the results of high--resolution N--body simulations to explore the consequences of working in redshift space on the detection of deviations from hierarchical clustering. Both for real and simulated data, significant deviations from hierarchical clustering seem to be present in real space. Their behaviour is coherent with that expected from an initially biased galaxy field, once the displacement from the sites where galaxies formed, due to the nonlinear gravitational evolution, is taken into account. The passage to redshift space has the net effect of filtering out the higher powers required to fit the distribution in real space. This magnifies the distortions operated by nonlinear evolution on the initial distribution, erasing the residual narrow scale range (2--5~$h^{-1}$Mpc) where deviations from hierarchical clustering are still detectable in real space. We conclude that such deviations can hardly be estimated using data in redshift space.

Although the Cold + Hot Dark Matter (CHDM) cosmology provides perhaps the best fit of any model to all the available data at the current epoch, CHDM produces structure at relatively low redshifts and thus could be ruled out if there were evidence for formation of massive objects at high redshifts. Damped Ly$\alpha$ systems are abundant in quasar absorption spectra and thus provide possibly the most significant evidence for early structure formation, and thus perhaps the most stringent constraint on CHDM. Using the numbers of halos in N-body simulations to normalize Press-Schechter estimates of the number densities of protogalaxies as a function of redshift, we find that CHDM with $\Omega_c/\Omega_\nu/\Omega_b = 0.6/0.3/0.1$ is compatible with the damped Ly$\alpha$ data at $\le 2.5$, but that it is probably incompatible with the limited $z>3$ damped Ly$\alpha$ data. The situation is uncertain because there is very little data for $z>3$, and also it is unclear whether all damped Ly$\alpha$ systems are associated with collapsed protogalaxies. The predictions of CHDM are quite sensitive to the hot (neutrino) fraction, and we find that $\Omega_c/\Omega_\nu/\Omega_b = 0.675/0.25/0.075$ is compatible even with the $z>3$ data. This corresponds to lowering the neutrino mass from 6.8 to 5.7 eV, for $H_0=50\kmsMpc$. In CHDM, the higher redshift damped Ly$\alpha$ systems are predicted to have lower masses, which can be checked by measuring the velocity widths of the associated metal line systems.

The distribution of dark matter around galactic or cluster halos has usually been assumed to be approximately isothermal with a non-zero core radius, which is expected to be of the order of the size of the visible matter distribution. Recently, the possibility has been raised that dark matter halos might be singular in the sense that the dark matter density $\rho$ could increase monotonically with radius $r$ down to a very small distance from the center of galaxies or clusters. Such central cusps in the dark matter density could lead to a high flux of gamma rays from WIMP dark matter annihilation. Here we analyze two possibilities that have been discussed in the literature, $\rho \propto r^{-n}$ with $n \approx 1\ {\rm or}\ 2$, and point out that such density profiles are excluded by gravitional lensing analyses on cluster scales and by the rotation curves of gas-rich, halo-dominated dwarf spirals on small scales. We also point out that if spiral galaxies form by gas infall inside dark matter halos, as they are expected to do in any hierarchical clustering model, such profiles almost always lead to falling rotation curves after infall, contrary to observations.

This letter presents results of new high resolution $\Omega=1$ Cold + Hot Dark Matter (CHDM) and Cold Dark Matter (CDM) simulations. Properties of groups in these simulations reflect the lower small-scale velocities and the greater tendency to form distinct filaments on both small and large scales in CHDM as compared to CDM. The fraction of galaxies in groups and the median group rms velocity are found to be powerful discriminators between models. We combine these two features into a very robust statistic, median group rms velocity $v_{\rm gr}(f_{\rm gr})$ as a function of the fraction $f_{\rm gr}$ of galaxies in groups. Using this statistic, we compare ``observed'' simulations to CfA data in redshift space in a careful and consistent way. We find that CHDM remains a promising model, with for example $v_{\rmgr}(0.45) \approx 125 \pm 25 \kms$ in agreement with the CfA data, while CDM with bias b=1.0 (COBE-compatible) or b=1.5, both giving $v_{\rm gr}(0.45) \approx 400 \pm 25 \kms$, can be virtually ruled out. Using median $M/L$, the observed value of $\Omega$ is $0.10$ (CHDM) to $0.38$ (CDM).

We use high resolution PM N-body simulations to follow the development of non-linear clustering in a flat Universe, dominated by Cold + Hot Dark Matter (CHDM) with 60% of CDM, 30% of HDM and 10% of baryons; a simulation box of 100 Mpc a side ($h=0.5$) is used. We analyze two CHDM simulations with $b =1.5$ (COBE normalization). We also compare them with CDM simulations with $b=1.5$ and $b=1$. We evaluate high-order correlation functions and the void-probability-function (VPF). Correlation functions are obtained both from counts in cells and counts of neighbors. The analysis is made for DM particles and for galaxies, identified as massive halos in the evolved density field. We also check the effects of dynamical evolution and redshift space distortions. We find that clustering of DM particles exhibits deviations from the hierarchical scaling, which decrease somewhat in redshift space. Galaxies follow hierarchical scaling far more closely, with coefficients S_3=2.5 and S_4=7.5, in general agreement with observational results, just marginally affected by redshift distortions and on the choices for the initial spectra. The hierarchical scaling of galaxy clustering is confirmed by the VPF analysis and in all the cases considered a good agreement with observational results is obtained.

We report results from high-resolution particle-mesh (PM) N-body simulations of structure formation in an $\Omega=1$ cosmological model with a mixture of Cold plus Hot Dark Matter (C+HDM) having $\Omega_{\rm cold}=0.6$, $\Omega_\nu=0.3$, and $\Omega_{\rm baryon}=0.1$. We present analytic fits to the C+HDM power spectra for both cold and hot ($\nu$) components, which provide initial conditions for our nonlinear simulations. In order to sample the neutrino velocities adequately, these simulations included six times as many neutrino particles as cold particles. Our simulation boxes were 14, 50, and 200~Mpc cubes (with $H_0=50$ km s$^{-1}$ Mpc$^{-1}$); we also did comparison simulations for Cold Dark Matter (CDM) in a 50~Mpc box. C+HDM with linear bias factor $b=1.5$ is consistent both with the COBE data and with the galaxy correlations we calculate. We find the number of halos as a function of mass and redshift in our simulations; our results for both CDM and C+HDM are well fit by a Press-Schechter model. The number density of galaxy-mass halos is smaller than for CDM, especially at redshift $z>2$, but the numbers of cluster-mass halos are comparable. We also find that on galaxy scales the neutrino velocities and flatter power spectrum in C+HDM result in galaxy pairwise velocities that are in good agreement with the data, and about 30\% smaller than in CDM with the same biasing factor. On scales of several tens of Mpc, the C+HDM streaming velocities are considerably larger than CDM. Thus C+HDM looks promising as a model of structure formation.

The IV Rencontres de Blois, on Particle Astrophysics, held at the Ch\^ateau de Blois, June 15-20, 1992, was a meeting well-timed for a reconsideration of the issues in particle astrophysics in the light of the COBE discovery of the cosmic microwave background (CMB) fluctuations. This is a summary of what I thought were the most interesting things discussed at Blois: (1) The near-success of Cold Dark Matter (CDM) in predicting the COBE fluctuation amplitude, which favors the hypothesis that structure formed in the universe through gravitational collapse. (2) The indications that $\Omega\approx1$ and that the power spectrum has a little more power on supercluster and larger scales than CDM. These are suggested by the IRAS and CfA redshift surveys and POTENT galaxy peculiar velocity analysis, and also by the COBE data. (3) The consequent demise of CDM and the rise of hybrid schemes such as Cold+Hot Dark Matter (C+HDM). (4) The possible implications for neutrino masses and mixings, and for cosmology, of the recent results on solar neutrinos. (5) The first discovery of TeV $\gamma$ rays from an extragalactic source, which was announced at Blois. I also summarize here a number of the exciting ongoing and planned experiments and observations discussed at Blois: CERN experiments on $\nu_\mu \nu_\tau$ oscillations, which may be sufficiently sensitive to detect the $\nu_\tau$ if