Search

We investigate the biases and uncertainties in estimates of physical parameters of high-redshift Lyman break galaxies (LBGs), such as stellar mass, mean stellar population age, and star formation rate (SFR), obtained from broad-band photometry. By combining LCDM hierarchical structure formation theory, semi-analytic treatments of baryonic physics, and stellar population synthesis models, we construct model galaxy catalogs from which we select LBGs at redshifts z ~ 3.4, 4.0, and 5.0. The broad-band spectral energy distributions (SEDs) of these model LBGs are then analysed by fitting galaxy template SEDs derived from stellar population synthesis models with smoothly declining SFRs. We compare the statistical properties of LBGs' physical parameters -- such as stellar mass, SFR, and stellar population age -- as derived from the best-fit galaxy templates with the intrinsic values from the semi-analytic model. We find some trends in these distributions: first, when the redshift is known, SED-fitting methods reproduce the input distributions of LBGs' stellar masses relatively well, with a minor tendency to underestimate the masses overall, but with substantial scatter. Second, there are large systematic biases in the distributions of best-fit SFRs and mean ages, in the sense that single-component SED-fitting methods underestimate SFRs and overestimate ages. We attribute these trends to the different star formation histories predicted by the semi-analytic models and assumed in the galaxy templates used in SED-fitting procedure, and to the fact that light from the current generation of star-formation can hide older generations of stars. These biases, which arise from the SED-fitting procedure, can significantly affect inferences about galaxy evolution from broadband photometry., Comment: 85 pages, 34 figures, submittted to ApJS

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., Comment: 12 pages, 8 figures, to be published in the Proceedings of the 4th Heidelberg International Symposium on High Energy Gamma-Ray Astronomy, held July 2008 in Heidelberg, Germany

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., Comment: 4 pages, 7 figures, To be published in the Proceedings of the 4th Heidelberg International Symposium on High Energy Gamma-Ray Astronomy, held July 2008 in Heidelberg, Germany

The treatment of dust attenuation is crucial in order to compare the predictions of galaxy formation models with multiwavelength observations. Most past studies have either used simple analytic prescriptions or else full radiative transfer (RT) calculations. Here, we couple star formation histories and morphologies predicted by the semi-analytic galaxy formation model MORGANA with RT calculations from the spectrophotometric and dust code GRASIL to create a library of galaxy SEDs from the UV/optical through the far Infrared, and compare the predictions of the RT calculations with analytic prescriptions. We consider a low and high redshift sample, as well as an additional library constructed with empirical, non-cosmological star formation histories and simple (pure bulge or disc) morphologies. Based on these libraries, we derive fitting formulae for the effective dust optical depth as a function of galaxy physical properties such as metallicity, gas mass, and radius. We show that such fitting formulae can predict the V-band optical depth with a scatter smaller than 0.4 dex for both the low and high redshift samples, but that there is a large galaxy-to-galaxy scatter in the shapes of attenuation curves, probably due to geometrical variations, which our simple recipe does not capture well. However, our new recipe provides a better approximation to the GRASIL results at optical wavelength than standard analytic prescriptions from the literature, particularly at high redshift., Comment: 20 pages, 13 figures, 4 tables; accepted for publication by MNRAS

We investigate the correlation of star formation quenching with internal galaxy properties and large scale environment (halo mass) in empirical data and theoretical models. We make use of the halo-based Group Catalog of Yang and collaborators, which is based on the Sloan Digital Sky Survey. Data from the Galaxy Evolution Explorer (GALEX) are also used to extract the recent star formation rate. In order to investigate the environmental effects, we examine the properties of "central" and "satellite" galaxies separately. For central galaxies, we are unable to conclude whether star formation quenching is primarily connected with halo mass or stellar mass, because these two quantities are themselves strongly correlated. For satellite galaxies, a nearly equally strong dependence on halo mass and stellar mass is seen. We make the same comparison for five different semi-analytic models based on three independently developed codes. We find that the models with AGN feedback reproduce reasonably well the dependence of the fraction of central red and passive galaxies on halo mass and stellar mass. However, for satellite galaxies, the same models badly overproduce the fraction of red/passive galaxies and do not reproduce the empirical trends with stellar mass or halo mass. This {\em satellite overquenching problem} is caused by the too-rapid stripping of the satellites' hot gas halos, which leads to rapid strangulation of star formation., Comment: 18 pages, 11 figures, Accepted for publication in MNRAS

The observational link between Supermassive Black Holes (SMBH) and galaxies at low redshift seems to be very tight, and statistically the global evolution of star formation activity and BH accretion activity also seem to trace each other closely. However, pinning down the co-evolution of galaxies and BH on an object-by-object basis remains elusive. I present results from new models for the joint evolution of galaxies, SMBH, and AGN, which may be able to help resolve some of the observational puzzles. A unique aspect of these models is our treatment of self-regulated BH growth based on hydrodynamic simulations of galaxy-galaxy mergers. Although these models do quite well at reproducing the observed evolution of galaxies, they do not reproduce the observed history of BH accretion, predicting too much early accretion and not enough at late times. I suggest two possible resolutions to this problem., Comment: to appear in "Panoramic Views of Galaxy Formation and Evolution", First Subaru International Conference, Hayama, Japan, December 2007

We present a new semi-analytic model that self-consistently traces the growth of supermassive black holes (BH) and their host galaxies within the context of the LCDM cosmological framework. In our model, the energy emitted by accreting black holes regulates the growth of the black holes themselves, drives galactic scale winds that can remove cold gas from galaxies, and produces powerful jets that heat the hot gas atmospheres surrounding groups and clusters. We present a comprehensive comparison of our model predictions with observational measurements of key physical properties of low-redshift galaxies, such as cold gas fractions, stellar metallicities and ages, and specific star formation rates. We find that our new models successfully reproduce the exponential cutoff in the stellar mass function and the stellar and cold gas mass densities at z~0, and predict that star formation should be largely, but not entirely, quenched in massive galaxies at the present day. We also find that our model of self-regulated BH growth naturally reproduces the observed relation between BH mass and bulge mass. We explore the global formation history of galaxies in our models, presenting predictions for the cosmic histories of star formation, stellar mass assembly, cold gas, and metals. We find that models assuming the "concordance" LCDM cosmology overproduce star formation and stellar mass at high redshift (z>2). A model with less small-scale power predicts less star formation at high redshift, and excellent agreement with the observed stellar mass assembly history, but may have difficulty accounting for the cold gas in quasar absorption systems at high redshift (z~3-4)., Comment: MNRAS accepted

The QSO luminosity function at z>5 provides strong constraints on models of joint evolution of QSO and their hosts. However, these observations are challenging because the low space densities of these objects necessitate surveying of large areas, in order to obtain statistically meaningful samples, while at the same time cosmological redshifting and dimming means that rather deep Near Infrared (NIR) imaging must be carried out. Several upcoming and proposed facilities with wide-field NIR imaging capabilities will open up this new region of parameter space. In this paper we present predictions for the expected number counts of z>5 QSOs, based on simple empirical and semi-empirical models of QSO evolution, as a function of redshift, depth and surveyed area. We compute the evolution of observed-frame QSO magnitudes and colors in a representative photometric system covering the wavelength range 550nm

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

We present the clustering properties of 151 Lyman alpha emitting galaxies at z ~ 4.5 selected from the Large Area Lyman Alpha (LALA) survey. Our catalog covers an area of 36' x 36' observed with five narrowband filters. We assume that the angular correlation function w(theta) is well represented by a power law A_w = Theta^(-beta) with slope beta = 0.8, and we find A_w = 6.73 +/- 1.80. We then calculate the correlation length r_0 of the real-space two-point correlation function xi(r) = (r/r_0)^(-1.8) from A_w through the Limber transformation, assuming a flat, Lambda-dominated universe. Neglecting contamination, we find r_0 = 3.20 +/- 0.42 Mpc/h. Taking into account a possible 28% contamination by randomly distributed sources, we find r_0 = 4.61 +/- 0.6 Mpc/h. We compare these results with the expectations for the clustering of dark matter halos at this redshift in a Cold Dark Matter model, and find that the measured clustering strength can be reproduced if these objects reside in halos with a minimum mass of 1-2 times 10^11 Solar masses/h. Our estimated correlation length implies a bias of b ~ 3.7, similar to that of Lyman-break galaxies (LBG) at z ~ 3.8-4.9. However, Lyman alpha emitters are a factor of ~ 2-16 rarer than LBGs with a similar bias value and implied host halo mass. Therefore, one plausible scenario seems to be that Lyman alpha emitters occupy host halos of roughly the same mass as LBGs, but shine with a relatively low duty cycle of 6-50%., Comment: 23 pages in preprint format, 4 figures, ApJ accepted

In the context of measuring structure and morphology of intermediate redshift galaxies with recent HST/ACS surveys, we tune, test, and compare two widely used fitting codes (GALFIT and GIM2D) for fitting single-component Sersic models to the light profiles of both simulated and real galaxy data. We find that fitting accuracy depends sensitively on galaxy profile shape. Exponential disks are well fit with Sersic models and have small measurement errors, whereas fits to de Vaucouleurs profiles show larger uncertainties owing to the large amount of light at large radii. We find that both codes provide reliable fits and little systematic error, when the effective surface brightness is above that of the sky. Moreover, both codes return errors that significantly underestimate the true fitting uncertainties, which are best estimated with simulations. We find that GIM2D suffers significant systematic errors for spheroids with close companions owing to the difficulty of effectively masking out neighboring galaxy light; there appears to be no work around to this important systematic in GIM2D's current implementation. While this crowding error affects only a small fraction of galaxies in GEMS, it must be accounted for in the analysis of deeper cosmological images or of more crowded fields with GIM2D. In contrast, GALFIT results are robust to the presence of neighbors because it can simultaneously fit the profiles of multiple companions thereby deblending their effect on the fit to the galaxy of interest. We find GALFIT's robustness to nearby companions and factor of >~20 faster runtime speed are important advantages over GIM2D for analyzing large HST/ACS datasets. Finally we include our final catalog of fit results for all 41,495 objects detected in GEMS., Comment: Accepted for publication in ApJS October 2007, v172n2; 25 pages, 16 Figures, 9 Tables; for hi-resolution version, see http://www.mpia.de/homes/bhaeussl/galaxy_fitting.pdf. For results, catalogues and files for code-testing, see http://www.mpia.de/GEMS/fitting_paper.html

Surveys of distant galaxies with the Hubble Space Telescope and from the ground have shown that there is only mild evolution in the relationship between radial size and stellar mass for galactic disks from z~1 to the present day. Using a sample of nearby disk-dominated galaxies from the Sloan Digital Sky Survey (SDSS), and high redshift data from the GEMS (Galaxy Evolution from Morphology and SEDs) survey, we investigate whether this result is consistent with theoretical expectations within the hierarchical paradigm of structure formation. The relationship between virial radius and mass for dark matter halos in the LCDM model evolves by about a factor of two over this interval. However, N-body simulations have shown that halos of a given mass have less centrally concentrated mass profiles at high redshift. When we compute the expected disk size-stellar mass distribution, accounting for this evolution in the internal structure of dark matter halos and the adiabatic contraction of the dark matter by the self-gravity of the collapsing baryons, we find that the predicted evolution in the mean size at fixed stellar mass since z~1 is about 15-20 percent, in good agreement with the observational constraints from GEMS. At redshift z~2, the model predicts that disks at fixed stellar mass were on average only 60% as large as they are today. Similarly, we predict that the rotation velocity at a given stellar mass (essentially the zero-point of the Tully-Fisher relation) is only about 10 percent larger at z~1 (20 percent at z~2) than at the present day., Comment: 13 pages, 6 figures, accepted for publication in ApJ. Revised in response to referee's comments to improve clariry. Results are unchanged

This paper presents the Hubble Ultra Deep Field (HUDF), a one million second exposure of an 11 square minute-of-arc region in the southern sky with the Hubble Space Telescope. The exposure time was divided among four filters, F435W (B435), F606W (V606), F775W (i775), and F850LP (z850), to give approximately uniform limiting magnitudes mAB~29 for point sources. The image contains at least 10,000 objects presented here as a catalog. Few if any galaxies at redshifts greater than ~4 resemble present day spiral or elliptical galaxies. Using the Lyman break dropout method, we find 504 B-dropouts, 204 V-dropouts, and 54 i-dropouts. Using these samples that are at different redshifts but derived from the same data, we find no evidence for a change in the characteristic luminosity of galaxies but some evidence for a decrease in their number densities between redshifts of 4 and 7. The ultraviolet luminosity density of these samples is dominated by galaxies fainter than the characteristic luminosity, and the HUDF reveals considerably more luminosity than shallower surveys. The apparent ultraviolet luminosity density of galaxies appears to decrease from redshifts of a few to redshifts greater than 6. The highest redshift samples show that star formation was already vigorous at the earliest epochs that galaxies have been observed, less than one billion years after the Big Bang., Comment: 44 pages, 18 figures, to appear in the Astronomical Journal October 2006

We present deep Hubble Space Telescope Advanced Camera for Surveys / High Resolution Channel U, narrow-V, and I images of nine "ultracompact" blue dwarf galaxies (UCBDs) selected from the Sloan Digital Sky Survey. We define UCBDs as local (z < 0.01) star-forming galaxies having angular diameters < 6" and physical sizes < 1 kpc. They are also among the most metal-poor galaxies known, and are found to reside within voids. Both the HST images and the objects' optical spectra reveal that they are composites of young (~1 Myr) populations that dominate their light, and older (~10 Gyr) populations that dominate their stellar masses, which we estimate to be ~10^7 - 10^8 Msol. An intermediate-age population is also indicated in most cases. The objects are not as dynamically disturbed as the prototype UCBD, POX 186, but the structure of several of them suggests that their current starbursts have been triggered by the collisions/mergers of smaller clumps of stars. In one case, HS 0822+3542, the ACS/HRC images resolve the object into two small (~100 pc) components which appear to have recently collided, supporting this interpretation. In six of the objects much of their star formation is concentrated in Young Massive Star clusters. The evidence that the galaxies consist mainly of ~10 Gyr old stars establishes that they are not protogalaxies; their low metallicities are more likely to be the result of the escape of supernova ejecta, as opposed to youth. These results are consistent with recent galaxy formation simulations which predict that cosmic re-ionization at z ~ 6 significantly limited the subsequent star formation of dwarf galaxies in voids due to the photo-evaporation of baryons from their cold dark matter halos (Abridged)., Comment: 43 pages, 6 figures, accepted for publication in ApJ; figures available at ftp://ftp.nofs.navy.mil/pub/outgoing/mcorbin/

Using a model for self-regulated growth of black holes (BHs) in mergers involving gas-rich galaxies, we study the relationship between quasars and the population of merging galaxies and predict the merger-induced star formation rate density of the Universe. Mergers drive nuclear gas inflows, fueling starbursts and 'buried quasars' until accretion feedback expels the gas, rendering a briefly visible optical quasar. Star formation is shut down and accretion declines, leaving a passively evolving remnant with properties typical of red, elliptical galaxies. Based on evolution of these events in our simulations, we demonstrate that the observed statistics of merger rates, luminosity functions (LFs) and mass functions, SFR distributions, specific SFRs, quasar and quasar host galaxy LFs, and elliptical/red galaxy LFs are self-consistent and follow from one another as predicted by the merger hypothesis. We use our simulations to de-convolve both quasar and merging galaxy LFs to determine the birthrate of black holes of a given final mass and merger rates as a function of stellar mass. We use this to predict the merging galaxy LF in several observed wavebands, color-magnitude relations, mass functions, absolute and specific SFR distributions and SFR density, and quasar host galaxy LFs, as a function of redshift from z=0-6. We invert this and predict e.g. quasar LFs from observed merger LFs or SFR distributions. Our results agree well with observations, but idealized models of quasar lightcurves are ruled out by comparison of merger and quasar observations at >99.9% confidence. Using only observations of quasars, we estimate the contribution of mergers to the SFR density of the Universe even to high redshifts z~4., Comment: 26 pages, 15 figures, matches version accepted to ApJ

We calculate the projected two point correlation function for samples of luminous and massive galaxies in the COMBO-17 photometric redshift survey, focusing particularly on the amplitude of the correlation function at small projected radii and exploring the constraints such measurements can place on the galaxy merger rate. For nearly volume-limited samples with 0.42.5e10 M_sun is 5+/-1%. Incorporating close pair fractions from the literature, the 2dFGRS and the SDSS, we find a fairly rapid evolution of the merger fraction of massive galaxies between z=0.8 and the present day. Assuming that the major merger timescale is of order the dynamical timescale for close massive galaxy pairs, we tentatively infer that ~50% (70%) of all galaxies with present-day masses M_*>5e10 M_sun (remnants of mergers between galaxies with M_*>2.5e10 M_sun) have undergone a major merger since z=0.8(1): major mergers between massive galaxies are a significant driver of galaxy evolution over the last eight billion years., Comment: ApJ, in press. 8 pages, 3 figures. Expanded discussion section with explicit discussion of merger fraction vs. close pair fraction. Change of typical close pair timescale results in increased inferred merger rate

We report on the angular correlation function of Lyman-break galaxies (LBGs) at z~4 and 5 from deep samples obtained from the Great Observatories Deep Origins Survey (GOODS). Similar to LBGs at z~3, the shape of w(theta) of the GOODS LBGs is well approximated by a power-law with slope beta~0.6 at angular separation theta > 10 arcsec. The clustering strength of z~4, 5 LBGs also depends on the rest-frame UV luminosity, with brighter galaxies more strongly clustered than fainter ones, implying a general correlation between halos' mass and LBGs' star-formation rate. At smaller separations, w(theta) of deep samples significantly exceeds the extrapolation of the large-scale power-law fit, implying enhanced spatial clustering at scales r < 1 Mpc. We also find that bright LBGs statistically have more faint companions on scales theta < 20 arcsec than fainter ones, showing that the enhanced small-scale clustering is very likely due to sub-structure, namely the fact that massive halos can host multiple galaxies. A simple model for the halo occupation distribution and the CDM halo mass function reproduce well the observed w(theta). The scaling relationship of the clustering strength with volume density and with redshift is quantitatively consistent with that of CDM halos. A comparison of the clustering strength of three samples of equal luminosity limit at z ~ 3, 4 and 5 shows that the LBGs at z~5 are hosted in halos about one order of magnitude less massive than those in the lower redshift bins, suggesting that star-formation was more efficient at higher-redshift., Comment: replaced with the version accepted for publication in ApJ. 46 pages, 10 figures; minor changes to text, one subsection added

We have used the 28'x 28' HST image mosaic from the GEMS (Galaxy Evolution from Morphology and SEDs) survey in conjunction with the COMBO-17 photometric redshift survey to constrain the incidence of major mergers between spheroid-dominated galaxies with little cold gas (dry mergers) since z = 0.7. A set of N-body merger simulations was used to explore the morphological signatures of such interactions: they are recognizable either as < 5kpc separation close pairs or because of broad, low surface brightness tidal features and asymmetries. Data with the depth and resolution of GEMS are sensitive to dry mergers between galaxies with M_V < -20.5 for z < 0.7; dry mergers at higher redshifts are not easily recovered in single-orbit HST imaging. Six dry mergers (12 galaxies) with luminosity ratios between 1:1 and 4:1 were found from a sample of 379 red early-type galaxies with M_V < -20.5 and 0.1 < z < 0.7. The simulations suggest that the morphological signatures of dry merging are visible for ~250Myr and we use this timescale to convert the observed merger incidence into a rate. On this basis we find that present day spheroidal galaxies with M_V < -20.5 on average have undergone between 0.5 and 2 major dry mergers since z ~ 0.7. We have compared this result with the predictions of a Cold Dark Matter based semi-analytic galaxy formation model. The model reproduces the observed declining major merger fraction of bright galaxies and the space density of luminous early-type galaxies reasonably well. The predicted dry merger fraction is consistent with our observational result. Hence, hierarchical models predict and observations now show that major dry mergers are an important driver of the evolution of massive early-type galaxies in recent epochs., Comment: ApJ, in press. The paper has been extensively modified, detailing the automated+visual selection and dry merger classification. 11 pages emulateapj with 9 reduced-quality figures. A high quality copy is available at http://www.mpia-hd.mpg.de/homes/bell/papers/dry.ps.gz

Using Hubble Space Telescope and ground-based U through K- band photometry from the Great Observatories Origins Deep Survey (GOODS), we measure the evolution of the luminosity function and luminosity density in the rest-frame optical (UBR) to z ~ 2, bridging the poorly explored ``redshift desert'' between z~1 and z~2. We also use deep near-infrared observations to measure the evolution in the rest-frame J-band to z~1. Compared to local measurements from the SDSS, we find a brightening of the characteristic magnitude, (M*), by ~2.1, \~0.8 and ~0.7 mag between z=0.1 and z=1.9, in U, B, and R bands, respectively. The evolution of M* in the J-band is in the opposite sense, showing a dimming between redshifts z=0.4 and z=0.9. This is consistent with a scenario in which the mean star formation rate in galaxies was higher in the past, while the mean stellar mass was lower, in qualitative agreement with hierarchical galaxy formation models. We find that the shape of the luminosity function is strongly dependent on spectral type and that there is strong evolution with redshift in the relative contribution from the different spectral types to the luminosity density. We find good agreement in the luminosity function derived from an R-selected and a K-selected sample at z~1, suggesting that optically selected surveys of similar depth (R < 24) are not missing a significant fraction of objects at this redshift relative to a near-infrared-selected sample. We compare the rest-frame B-band luminosity functions from z~0--2 with the predictions of a semi-analytic hierarchical model of galaxy formation, and find qualitatively good agreement. In particular, the model predicts at least as many optically luminous galaxies at z~1--2 as are implied by our observations., Comment: 43 pages; 15 Figures; 5 Tables, Accepted for publication in Ap.J

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., Comment: update to match published version

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., Comment: 12 pages, 10 figures, submitted to ApJ

We combine HST imaging from the GEMS survey with photometric redshifts from COMBO-17 to explore the evolution of disk-dominated galaxies since z<1.1. The sample is comprised of all GEMS galaxies with Sersic indices n<2.5, derived from fits to the galaxy images. We account fully for selection effects through careful analysis of image simulations; we are limited by the depth of the redshift and HST data to the study of galaxies with absolute magnitudes M(V)<-20, or equivalently stellar masses log(M)>10. We find strong evolution in the magnitude-size scaling relation for galaxies with M(V)<-20, corresponding to a brightening of 1 mag per sqarcsec in rest-frame V-band by z=1. Yet, disks at a given absolute magnitude are bluer and have lower stellar mass-to-light ratios at z=1 than at the present day. As a result, our findings indicate weak or no evolution in the relation between stellar mass and effective disk size for galaxies with log(M)>10 over the same time interval. This is strongly inconsistent with the most naive theoretical expectation, in which disk size scales in proportion to the halo virial radius, which would predict that disks are a factor of two denser at fixed mass at z=1. The lack of evolution in the stellar mass-size relation is consistent with an ``inside-out'' growth of galaxy disks on average (galaxies increasing in size as they grow more massive), although we cannot rule out more complex evolutionary scenarios., Comment: 22 pages, 16 figures, submitted to ApJ

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., Comment: 11 pages, 7 figures, to appear in the proceedings of the Gamma 2004 Symposium on High Energy Gamma Ray Astronomy, Heidelberg, July 2004 (AIP Proceedings Series)

We combine HST/ACS imaging from the GEMS survey with redshifts and rest-frame quantities from COMBO-17 to study the evolution of morphologically early-type galaxies with red colors since z=1. We use a new large sample of 728 galaxies with centrally-concentrated radial profiles (Sersic n>2.5) and rest-frame U-V colors on the red sequence. By appropriate comparison with the local relations from SDSS, we find that the luminosity-size (L-R) and stellar mass-size (M-R) relations evolve in a manner that is consistent with the passive aging of ancient stars. By itself, this result is consistent with a completely passive evolution of the red early-type galaxy population. If instead, as demonstrated by a number of recent surveys, the early-type galaxy population builds up in mass by a factor of 2 since z=1, our results imply that new additions to the early-type galaxy population follow similar L-R and M-R correlations, compared to the older subset of early-type galaxies. Adding early-type galaxies to the red sequence through disk fading appears to be consistent with the data. Through comparison with models, the role of dissipationless merging is limited to <1 major merger on average since z=1 for the most massive galaxies. Predictions from models of gas-rich mergers are not yet mature enough to allow a detailed comparison to our observations. We find tentative evidence that the amount of luminosity evolution depends on galaxy stellar mass, such that the least massive galaxies show stronger luminosity evolution compared to more massive early types. This could reflect a different origin of low-mass early-type galaxies and/or younger stellar populations; the present data is insufficient to discriminate between these possibilities. (abridged), Comment: Submitted to ApJ, 23 pages, Latex using emulateapj5.sty and onecolfloat.sty (included), 10 figures, version with full resolution figures at http://www.astro.umass.edu/~dmac/Papers/ETevol.hires.ps

Pan-chromatic galaxy surveys are providing tightening constraints on the global mass assembly history, and high-resolution imaging of large fields is telling us when and where stars were formed. How well are state-of-the-art hierarchical galaxy formation models currently doing at reproducing these observations? I present results here that suggest that hierarchical models are doing quite well at reproducing the global star formation and stellar mass assembly history obtained from galaxies selected in the optical and Near IR. However, the same models fail to reproduce two very important populations at high redshift: quiescent red spheroids and vigorously star-forming, dust-enshrouded starbursts. This mismatch carries important lessons about how star formation is triggered and regulated in early galaxies, and may force us to consider new ideas about the formation of massive spheroids., Comment: 8 pages, 3 figures, to appear in the Proceedings of the ESO/USM/MPE Workshop on "Multiwavelength Mapping of Galaxy Formation and Evolution", eds. R. Bender and A. Renzini

GEMS, Galaxy Evolution from Morphologies and SEDs, is a large-area (800 arcmin2) two-color (F606W and F850LP) imaging survey with the Advanced Camera for Surveys on HST. Centered on the Chandra Deep Field South, it covers an area of ~28'x28', or about 120 Hubble Deep Field areas, to a depth of m_AB(F606W)=28.3 (5sigma and m_AB(F850LP)=27.1 (5sigma) for compact sources. In its central ~1/4, GEMS incorporates ACS imaging from the GOODS project. Focusing on the redshift range 0.2<=z<=1.1, GEMS provides morphologies and structural parameters for nearly 10,000 galaxies where redshift estimates, luminosities and SEDs exist from COMBO-17. At the same time, GEMS contains detectable host galaxy images for several hundred faint AGN. This paper provides an overview of the science goals, the experiment design, the data reduction and the science analysis plan for GEMS., Comment: 24 pages, TeX with 6 eps Figures; to appear in ApJ Supplement. Low resolution figures only. Full resolution at http://zwicky.as.arizona.edu/~rix/Misc/GEMS.ps.gz

The integrated colors of distant galaxies provide a means for interpreting the properties of their stellar content. Here, we use rest-frame UV-to-optical colors to constrain the spectral-energy distributions and stellar populations of color-selected, B-dropout galaxies at z ~ 4 in the Great Observatories Origins Deep Survey. We combine the ACS data with ground-based near-infrared images, which extend the coverage of galaxies at z ~ 4 to the rest-frame B-band. We observe a color-magnitude trend in the rest-frame m(UV) - B versus B diagram for the z ~ 4 galaxies that has a fairly well-defined "blue-envelope", and is strikingly similar to that of color-selected, U-dropout galaxies at z ~ 3. We also find that although the co-moving luminosity density at rest-frame UV wavelengths (1600 Angstroms) is roughly comparable at z ~ 3 and z ~ 4, the luminosity density at rest-frame optical wavelengths increases by about one-third from z ~ 4 to z ~ 3. Although the star-formation histories of individual galaxies may involve complex and stochastic events, the evolution in the global luminosity density of the UV-bright galaxy population corresponds to an average star-formation history with a star-formation rate that is constant or increasing over these redshifts. This suggests that the evolution in the luminosity density corresponds to an increase in the stellar-mass density of more than 33%., Comment: Accepted for publication in ApJ Letters, 4 pages, 3 figures, uses emulateapj style. GOODS data and the related papers can be found at http://www.stsci.edu/science/goods/

Using U'- through Ks-band imaging data in the GOODS-South field, we construct a large, complete sample of 275 ``extremely red objects'' (EROs; K_s<22.0, R-K_s>3.35; AB), all with deep HST/ACS imaging in B_435, V_606, i_775, and z_850, and well-calibrated photometric redshifts. Quantitative concentration and asymmetry measurements fail to separate EROs into distinct morphological classes. We therefore visually classify the morphologies of all EROs into four broad types: ``Early'' (elliptical-like), ``Late'' (disk galaxies), ``Irregular'' and ``Other'' (chain galaxies and low surface brightness galaxies), and calculate their relative fractions and comoving space densities. For a broad range of limiting magnitudes and color thresholds, the relative number of early-type EROs is approximately constant at 33-44%, and the comoving space densities of Early- and Late-type EROs are comparable. Mean rest-frame spectral energy distributions (SEDs) at wavelengths between 0.1 and 1.2 um are constructed for all EROs. The SEDs are extremely similar in their range of shapes, independent of morphological type. The implication is that any differences between the broad-band SEDs of Early-type EROs and the other types are relatively subtle, and there is no robust way of photometrically distinguishing between different morphological types with usual optical/near-infrared photometry., Comment: Submitted to the ApJL. A version with full-resolution figures, all GOODS data and all GOODS collaboration papers may be found at http://www.stsci.edu/science/goods/

Cosmic variance is the uncertainty in observational estimates of the volume density of extragalactic objects such as galaxies or quasars arising from the underlying large-scale density fluctuations. This is often a significant source of uncertainty, especially in deep galaxy surveys, which tend to cover relatively small areas. We present estimates of the relative cosmic variance for one-point statistics (i.e. number densities) for typical scales and volumes sampled by the Great Observatories Origins Deep Survey (GOODS). We use two approaches: for objects with a known two-point correlation function that is well-approximated by a power law, one can use the standard analytic formalism to calculate the cosmic variance (in excess of shot noise). We use this approach to estimate the cosmic variance for several populations that are being studied in the GOODS program: Extremely Red Objects (ERO) at z~1, and Lyman Break Galaxies (LBG) at z~3 and z~4, using clustering information for similar populations in the literature. For populations with unknown clustering, one can use predictions from Cold Dark Matter theory to obtain a rough estimate of the variance as a function of number density. We present a convenient plot which allows one to use this approach to read off the cosmic variance for a population with a known mean redshift and estimated number density. We conclude that for the volumes sampled by GOODS, cosmic variance is a significant source of uncertainty for strongly clustered objects (~40-60% for EROs) and less serious for less clustered objects, ~10-20% for LBGs., Comment: 5 pages, 3 figures, to appear in the special GOODS issue of ApJ Letters. GOODS data and all GOODS collaboration papers may be found at http://www.stsci.edu/science/goods/

The redshift distribution of near-IR selected galaxies is often used to attempt to discriminate between the classical view of galaxy formation, in which present-day luminous galaxies were assembled at early times and evolve due to the passive aging of their stellar populations, and that of hierarchical structure formation, in which galaxies were assembled more recently via the merging of smaller objects. We carry out such a test here, by computing the distribution of photometric redshifts of K_AB<22 galaxies in the Great Observatories Origins Deep Field Survey GOODS Southern Field, and comparing the results with predictions from a semi-analytic model based on hierarchical structure formation, and a classical `passive evolution' model. We find that the redshift distributions at z <~ 1.5 of both the hierarchical and passive models are very similar to the observed one. At z>~1.5, the hierarchical model shows a deficit of galaxies, while the passive model predicts an excess. We investigate the nature of the observed galaxies in the redshift range where the models diverge, and find that the majority have highly disturbed morphologies, suggesting that they may be merger-induced starbursts. While the hierarchical model used here does not produce these objects in great enough numbers, the appearance of this population is clearly in better qualitative agreement with the hierarchical picture than with the classical passive evolution scenario. We conclude that the observations support the general framework of hierarchical formation, but suggest the need for new or modified physics in the models., Comment: 5 pages, 2 figures. To appear in the GOODS special edition of ApJ Letters. GOODS data and all GOODS collaboration papers may be found at http://www.stsci.edu/science/goods/

Hubble Space Telescope images of high-redshift galaxies selected via color and photometric redshifts are used to examine the size and axial-ratio distribution of galaxies as a function of redshift at lookback times t > 8 Gyr. These parameters are measured at rest-frame UV wavelengths (1200-2000 Angstroms) on images with a rest-frame resolution of less than 0.8 kpc. Galaxy radii are found to scale with redshift approximately as the inverse of the Hubble parameter H(z). This is in accord with the theoretical expectation that the typical sizes of the luminous parts of galaxies should track the expected evolution in the virial radius of dark-matter halos. The mean ratio of semi-major to semi-minor axis for a bright well-resolved sample of galaxies at z ~ 4 is b/a = 0.65, suggesting that these Lyman break galaxies are not drawn from a spheroidal population. However the median concentration index of this sample is C = 3.5, which is closer to the typical concentration indices, C ~ 4, of nearby elliptical galaxies than to the values, C<2 for local disk galaxies of type Sb and later., Comment: 5 pages, including 3 figures; uses emulateapj style; accepted by ApJL for publication in a special issue of early GOODS results

Reducing the power on small scales relative to the `standard' LCDM model alleviates a number of possible discrepancies with observations, and is favored by the recent analysis of WMAP plus galaxy and Lyman-alpha forest data. Here, we investigate the epoch of reionization in several models normalized to WMAP on large scales, and with sufficiently reduced power on small scales to solve the halo concentration and substructure problems. These include a tilted model, the WMAP running-index model, and a Warm Dark Matter model. We assume that the Universe was reionized by stellar sources comprised of a combination of supermassive (~200 M_sun) Pop III stars and Pop II stars with a `normal' IMF. We find that in all of these models, structure formation and hence reionization occurs late, certainly at redshifts below ten, and more probably at z < 6. This is inconsistent (at 2 sigma) with the determination of z_reion ~ 17 from the WMAP temperature-polarization data and is only marginally consistent with SDSS quasar observations. The tension between the galactic-scale observations, which favor low-power models, and the early reionization favored by WMAP can only be resolved if the efficiency of Pop III star formation is dramatically higher than any current estimate, or if there is an exotic population of ionizing sources such as mini-quasars. Otherwise, we may have to live with the standard LCDM power spectrum, and solve the small-scale problems in some other way., Comment: submitted to ApJL

We calculate the global star formation rate density (SFRD) from z ~ 30-3 using a semi-analytic model incorporating the hierarchical assembly of dark matter halos, gas cooling via atomic hydrogen, star formation, supernova feedback, and suppression of gas collapse in small halos due to the presence of a photoionizing background. We compare the results with the predictions of simpler models based on the rate of dark matter halo growth and a fixed ratio of stellar-to-dark mass, and with observational constraints on the SFRD at 3 < z < 6. We also estimate the star formation rate due to very massive, metal-free Pop III stars using a simple model based on the halo formation rate, calibrated against detailed hydrodynamic simulations of Pop III star formation. We find that the total production rate of hydrogen-ionizing photons during the probable epoch of reionization (15 < z < 20) is approximately equally divided between Pop II and Pop III stars, and that if reionization is late (less than about 15, close to the lower limit of the range allowed by the WMAP results), then Pop II stars alone may be able to reionize the Universe., Comment: submitted to ApJL

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

We investigate the properties of a galaxy similar to the Milky Way within the context of standard disk formation theory in a LCDM universe. Using the standard assumption that baryons conserve specific angular momentum when they collapse, we conclude that the mean properties of the model galaxies are in good agreement with the Milky Way and other similar spiral galaxies, but the predicted scatter in disk scale lengths may be too large. A model in which half of the initial specific angular momentum is transfered to the dark matter may produce a smaller scatter, if very compact disks are unstable and evolve into spheroids or early type galaxies., Comment: 3 pages, 2 figures, Latex, to appear in the proceedings of "Galaxy Evolution: Theory and Observations", Cozumel, Mexico, April 2002. Eds. V. Avila-Reese, C. Firmani, C. Frenk, & C. Allen, RevMexAA SC (2002)

In this paper, I discuss the capabilities and limitations of an 8-10 meter ultraviolet/optical telescope in space, the proposed successor to the Hubble Space Telescope, in the context of galaxy studies. The exquisite spatial resolution and excellent sensitivity of such a facility would open up new possibilities for the study of nearby dwarf galaxies ($z \la 0.5$), and for studying the internal structure and kinematics of more luminous galaxies at high redshift ($z \ga 2$). These applications are of particular importance because they would address areas in which the popular Cold Dark Matter theory is in potential conflict with observations., Comment: 8 pages, 3 figures, Latex, uses newpasp.sty. To appear in "Hubble's Science Legacy: Future Optical-Ultraviolet Astronomy from Space", eds. K.R. Sembach, J.C. Blades, G.D. Illingworth, & R.C. Kennicutt, ASP Conference Series

A new approach to estimating photometric redshifts - using Artificial Neural Networks (ANNs) - is investigated. Unlike the standard template-fitting photometric redshift technique, a large spectroscopically-identified training set is required but, where one is available, ANNs produce photometric redshift accuracies at least as good as and often better than the template-fitting method. The Bayesian priors on the underlying redshift distribution are automatically taken into account. Furthermore, inputs other than galaxy colours - such as morphology, angular size and surface brightness - may be easily incorporated, and their utility assessed. Different ANN architectures are tested on a semi-analytic model galaxy catalogue and the results are compared with the template-fitting method. Finally the method is tested on a sample of ~ 20000 galaxies from the Sloan Digital Sky Survey. The r.m.s. redshift error in the range z < 0.35 is ~ 0.021., Comment: Submitted to MNRAS, 9 pages, 9 figures, substantial improvements to paper structure

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., Comment: 6 pages, 4 figures, for Proceedings, Mass of Galaxies at Low and High Redshift, Venice, October 2001, eds. R. Bender and A. Renzini (ESO Astrophysics Symposia, Springer-Verlag)

Wide field near-infrared (IR) surveys have revealed a population of galaxies with very red optical$-$IR colors, which have been termed ``Extremely Red Objects'' (EROs). Modeling suggests that such red colors (R-K > 5) could be produced by galaxies at z>~1 with either very old stellar populations or very high dust extinction. Recently it has been discovered that EROs are strongly clustered. Are these objects the high-redshift progenitors of present day giant ellipticals (gEs)? Are they already massive at this epoch? Are they the descendents of the $z\sim3$ Lyman Break Galaxies (LBG), which have also been identified as possible high redshift progenitors of giant ellipticals? We address these questions within the framework of the Cold Dark Matter paradigm using an analytic model that connects the number density and clustering or bias of an observed population with the halo occupation function (the number of observed galaxies per halo of a given mass). We find that EROs reside in massive dark matter halos, with average mass > 1E13/h100 Msun. The occupation function that we derive for EROs is very similar to the one we derive for z=0, L>L* early type galaxies, whereas the occupation function for LBGs is skewed towards much smaller host halo masses ( ~ 1E11 to 1E12/h100 Msun. We then use the derived occupation function parameters to explore the possible evolutionary connections between these three populations., Comment: Replaced to match version accepted for publication in ApJ. New model added; appendix added with dark matter correlation function fits. 12 pages, uses emulateapj5.sty

We apply standard disk formation theory with adiabatic contraction within cuspy halo models predicted by the standard LCDM cosmology. The resulting models score remarkably well when confronted with the broad range of observational data available for the Milky Way and M31 galaxies, giving a Milky Way virial mass of 1-2x10^12Msun and concentration C=12-17. We consider two types of models, in which: (A) baryons conserve angular momentum and (B) some of the angular momentum of the baryons is transferred to the dark matter. Type-A models produce good agreement with observed rotation curves and obey constraints in the solar neighborhood, but may have too much dark matter in the center to allow a fast rotating bar. The type-B models with angular momentum transport have a slightly more massive disk and less dark matter in the central part, allowing a fast rotating bar to persist. Both classes of models probably have sufficient baryonic mass in the central 3.5kpc to reproduce recent observational values of the optical depth to microlensing events towards the Galactic center. All models require that about 1/2 of all baryons expected inside the virial radius must not be in the disk or bulge. We investigate whether the range of virial masses allowed by our dynamical models is compatible with constraints from the galaxy luminosity function, and find a range of parameter space that is allowed by this constraint. We conclude that rotation curves and dynamical properties of ``normal'' high surface brightness spiral galaxies appear to be consistent with standard LCDM., Comment: 26 pages, 11 figures, submited to ApJ

There is accumulating evidence that the faint end of the galaxy luminosity function might be very different in different locations. The luminosity function might be rising in rich clusters and flat or declining in regions of low density. If galaxies form according to the model of hierarchical clustering then there should be many small halos compared to the number of big halos. If this theory is valid then there must be a mechanism that eliminates at least the visible component of galaxies in low density regions. A plausible mechanism is photoionization of the intergalactic medium at a time before the epoch that most dwarf galaxies form in low density regions but after the epoch of formation for similar systems that ultimately end up in rich clusters. The dynamical timescales are found to accommodate this hypothesis in a flat universe with Omega_m < 0.4. If small halos exist but simply cannot be located because they have never become the sites of significant star formation, they still might have dynamical manifestations. These manifestations are hard to identify in normal groups of galaxies because small halos do not make a significant contribution to the global mass budget. However, it could be entertained that there are clusters of halos where there are only small systems, clusters that are at the low mass end of the hierarchical tree. There may be places where only a few small galaxies managed to form, enough for us to identify and use as test probes of the potential. It turns out that such environments might be common. Four probable groups of dwarfs are identified within 5 Mpc and the assumption they are gravitationally bound suggests M/L_B ~ 300 - 1200 M_sun/L_sun, 6 +/- factor 2 times higher than typical values for groups with luminous galaxies., Comment: Accepted ApJ 569, (April 20), 2002, 12 pages, 6 figures, 1 table

Cold Dark Matter theory predicts that the Local Group should contain many more dwarf-sized objects than the observed number of dwarf galaxies --- the so-called sub-structure problem. We investigate whether the suppression of star formation in these small objects due to the presence of a photoionizing background can resolve the problem. We make use of results from recent hydrodynamic simulations to build a recipe for the suppression of gas infall into semi-analytic galaxy formation models, and use these to predict the luminosity function of dwarf galaxies in the Local Group. In the models without photoionization ``squelching'', we predict a large excess of faint dwarf galaxies compared with the observed number in the Local Group --- thus, the usual recipe for supernovae feedback used in semi-analytic models does not solve the sub-structure problem on its own. When we include photoionization squelching, we find good agreement with the observations. We have neglected tidal destruction, which probably further reduces the number of dwarf galaxies. We conclude that photoionizing squelching easily solves the sub-structure problem. In fact, it is likely that once this effect is taken into account, models with reduced small-scale power (e.g. Warm Dark Matter) would underproduce dwarf galaxies., Comment: LaTeX, 5 pages, 2 figures, Submitted to ApJ Letters

We discuss how current and future data on the clustering and number density of z~3 Lyman-break galaxies (LBGs) can be used to constrain their relationship to dark matter haloes. We explore a three-parameter model in which the number of LBGs per dark halo scales like a power-law in the halo mass: N(M) = (M/M_1)^S for M>M_m. Here, M_m is the minimum mass halo that can host an LBG, M_1 is a normalization parameter, associated with the mass above which haloes host more than one observed LBG, and S determines the strength of the mass dependence. We show how these three parameters are constrained by three observable properties of LBGs: the number density, the large-scale bias, and the fraction of objects in close pairs. Given these three quantities, the three unknown model parameters may be estimated analytically, allowing a full exploration of parameter space. As an example, we assume an LCDM cosmology and consider the observed properties of a recent sample of spectroscopically confirmed LBGs. We find that the favored range for our model parameters is M_m = (0.4-8)x10^10 Msun/h, M_1 = (6-10)x10^12 Msun/h, and 0.9 < S < 1.1. The preferred region in M_m expands by an order of magnitude and slightly shallower slopes are acceptable if the allowed range of bias is permitted to span all recent observational estimates. We also discuss how the observed clustering of LBGs as a function of luminosity can be used to constrain halo occupation, although due to current observational uncertainties we are unable to reach any strong conclusions. Our methods and results can be used to constrain more realistic models that aim to derive the occupation function N(M) from first principles, and offer insight into how basic physical properties affect the observed properties of LBGs., Comment: 13 pages, 5 figures, Submitted to MNRAS, minor corrections

We model the acquisition of spin by dark-matter halos in semi-analytic merger trees. We explore two different algorithms; one in which halo spin is acquired from the orbital angular momentum of merging satellites, and another in which halo spin is gained via tidal torquing on shells of material while still in the linear regime. We find that both scenarios produce the characteristic spin distribution of halos found in N-body simulations, namely, a log-normal distribution with mean ~0.04 and standard deviation ~0.5 in the log. A perfect match requires fine-tuning of two free parameters. Both algorithms also reproduce the general insensitivity of the spin distribution to halo mass, redshift and cosmology seen in N-body simulations. The spin distribution can be made strictly constant by physically motivated scalings of the free parameters. In addition, both schemes predict that halos which have had recent major mergers have systematically larger spin values. These algorithms can be implemented within semi-analytic models of galaxy formation based on merger trees. They yield detailed predictions of galaxy properties that strongly depend on angular momentum (such as size and surface brightness) as a function of merger history and environment., Comment: 11 pages, 6 figures, submitted to MNRAS. revised version adds additional reference to Vitvitska et al (2001) and fixes minor typos. Revised version because of bug found in calculations which leads to changes in the values of free parameters introduced in the paper. The content of the paper remains unchanged

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., Comment: 16 pages, 7 figures, to appear in Proceedings of the International Symposium on Gamma-Ray Astronomy, Heidelberg, June 2000, eds F. Aharonian and H. V\"olk, AIP Conf. Proc

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., Comment: 20 pages, Latex, uses emulateapj5.sty, minor changes to match ApJ in press version

Do disk galaxies exist at redshifts much greater than unity, and how might they look different from local disks? How does the morphological mix of galaxies change with redshift? What can we learn from current observations about the properties of high redshift galactic disks? I present theoretical predictions based on semi-analytic hierarchical models, focussing on the role played by interactions, mergers, and starbursts in determining the observable properties of disk galaxies at high redshift, and discuss the interpretation of high redshift observations of possible proto-disks (damped Lyman-$\alpha$ systems and Lyman-break galaxies) in light of these predictions., Comment: 8 pages, 3 figures, Invited review to appear in the proceedings of "Galaxy Disks and Disk Galaxies", J.G. Funes and E.M. Corsini, eds

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

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., Comment: 4 pages, 2 figures to appear in The Hy-Redshift Universe: Galaxy Formation and Evolution at High Redshift Edited by A.J. Bunker and W.J.M. van Breugel, A.S.P. Conference Seriers, Vol. 193

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., Comment: 2 pages, 1 figure to appear in Proceedings of Rencontres Internationales de l'IGRAP, Clustering at High Redshift, Edited by A. Mazure, O. Le Fevre and V. Le Brun