Different methodologies lead to order-of-magnitude variations in predicted galaxy merger rates. We examine and quantify the dominant uncertainties. Different halo merger rates and subhalo 'destruction' rates agree to within a factor ~2 given proper care in definitions. If however (sub)halo masses are not appropriately defined or are under-resolved, the major merger rate can be dramatically suppressed. The dominant differences in galaxy merger rates owe to baryonic physics. Hydrodynamic simulations without feedback and older models that do not agree with the observed galaxy mass function propagate factor ~5 bias in the resulting merger rates. However, if the model matches the galaxy mass function, properties of central galaxies are sufficiently converged to give small differences in merger rates. But variations in baryonic physics of satellites also have dramatic effects. The known problem of satellite 'over-quenching' in most semi-analytic models (SAMs), whereby SAM satellites are too efficiently stripped of gas, could lead to order-of-magnitude under-estimates of merger rates for low-mass, gas-rich galaxies. Fixing the satellite properties to observations tends to predict higher merger rates, but with factor ~2 empirical uncertainties. Choice of mass ratio definition matters: at low masses, most true major mergers (in baryonic/dynamical galaxy mass) will appear to be minor mergers in their stellar or luminosity mass ratio. Observations and models using these criteria may underestimate major merger rates by factors ~5. Orbital parameters and gas fractions also introduce factor ~3 differences in amount of bulge formed by mergers, even for fixed mass ratio encounters.

We present high-resolution interferometric Submillimeter Array (SMA) imaging at 890 microns (~2" resolution) of two millimeter selected galaxies -- MMJ100015+021549 and MMJ100047+021021 -- discovered with the Max-Planck Millimeter Bolometer (MAMBO) on the IRAM 30 m telescope and also detected with Bolocam on the CSO, in the COSMOS field. The first source is significantly detected at the ~11 sigma level, while the second source is tentatively detected at the ~4 sigma level, leading to a positional accuracy of ~0.2-0.3". MM100015+021549 is identified with a faint radio and K-band source. MMJ100047+021021 shows no radio emission and is tentatively identified with a very faint K-band peak which lies at ~1.2" from a clumpy optical source. The submillimeter-to-radio flux ratio for MM100015+021549 yields a redshift of ~4.8, consistent with the redshift implied by the UV-to-submillimeter photometry, z~3.0-5.0. We find evidence for warm dust in this source with an infrared luminosity in the range ~0.9-2.5x10^{13} L_sun, supporting the increasing evidence for a population luminous submillimeter galaxies at z>3. Finally, the lack of photometric data for MMJ100047+021021 does not allow us to investigate its properties in detail, however its submillimeter-to-radio ratio implies z>3.5.

We present continent-scale VLBI - obtained with the European VLBI Network (EVN) at a wavelength of 18cm - of six distant, luminous submm-selected galaxies (SMGs). Our images have a synthesized beam width of ~30 milliarcsec FWHM - three orders of magnitude smaller in area than the highest resolution VLA imaging at this wavelength - and are capable of separating radio emission from ultra-compact radio cores (associated with active super-massive black holes - SMBHs) from that due to starburst activity. Despite targeting compact sources - as judged by earlier observations with the VLA and MERLIN - we identify ultra-compact cores in only two of our targets. This suggests that the radio emission from SMGs is produced primarily on larger scales than those probed by the EVN, and therefore is generated by star formation rather than an AGN - a result consistent with other methods used to identify the presence of SMBHs in these systems.

We present a physical model for the origin of z~2 Dust-Obscured Galaxies (DOGs), a class of high-redshift ULIRGs selected at 24 micron which are particularly optically faint (24/R>1000). By combining N-body/SPH simulations of high redshift galaxy evolution with 3D polychromatic dust radiative transfer models, we find that luminous DOGs (with F24 > 0.3 mJy at z~2 are well-modeled as extreme gas-rich mergers in massive (~5x10^12-10^13 Msun) halos, with elevated star formation rates (~500-1000 Msun/yr) and/or significant AGN growth (Mdot > 0.5 Msun/yr), whereas less luminous DOGs are more diverse in nature. At final coalescence, merger-driven DOGs transition from being starburst dominated to AGN dominated, evolving from a "bump" to a power-law shaped mid-IR (IRAC) spectral energy distribution (SED). After the DOG phase, the galaxy settles back to exhibiting a "bump" SED with bluer colors and lower star formation rates. While canonically power-law galaxies are associated with being AGN-dominated, we find that the power-law mid-IR SED can owe both to direct AGN contribution, as well as to a heavily dust obscured stellar bump at times that the galaxy is starburst dominated. Thus power-law galaxies can be either starburst or AGN dominated. Less luminous DOGs can be well-represented either by mergers, or by massive ($M_{\rm baryon} ~5x10^11 Msun) secularly evolving gas-rich disc galaxies (with SFR > 50 Msun/yr). By utilising similar models as those employed in the SMG formation study of Narayanan et al. (2010), we investigate the connection between DOGs and SMGs. We find that the most heavily star-forming merger driven DOGs can be selected as Submillimetre Galaxies (SMGs), while both merger-driven and secularly evolving DOGs typically satisfy the BzK selection criteria.

We present a physical model for origin of the cosmic diffuse infrared background (CDIRB). By utilizing the observed stellar mass function and its evolution as input to a semi-empirical model of galaxy formation, we isolate the physics driving diffuse IR emission. The model includes contributions from three primary sources of IR emission: steady-state star formation owing to isolated disk galaxies, interaction-driven bursts of star formation owing to close encounters and mergers, and obscured active galactic nuclei (AGN). We find that most of the CDIRB is produced by equal contributions from objects at z=0.5-1 and z>1, as suggested by recent observations. Of those sources, the vast majority of the emission originates in systems with low to moderate IR luminosities (L_{IR}<10^{12} $L_sun); the most luminous objects contribute significant flux only at high-redshifts (z>2). All star formation in ongoing mergers accounts for <10% of the total at all wavelengths and redshifts, while emission directly attributable to the interaction-driven burst itself accounts for <5%. We furthermore find that obscured AGN contribute <1-2% of the CDIRB at all wavelengths and redshifts, with a strong upper limit of less than 4% of the total emission. Finally, since electron-positron pair production interactions with the CDIRB represent the primary source of opacity to very high energy (VHE: E_\gamma > 1 TeV) \gamma-rays, the model provides predictions for the optical depth of the Universe to the most energetic photons. We find that these predictions agree with observations of high-energy cutoffs at TeV energies in nearby blazars, and suggest that while the Universe is extremely optically thick at >10 TeV, the next generation of VHE \gamma-ray telescopes can reasonably expect detections from out to 50-150 Mpc.

We present high-resolution 345 GHz interferometric observations of two extreme luminous (L_{IR}>10^{13} L_sun), submillimetre-selected galaxies (SMGs) in the COSMOS field with the Submillimeter Array (SMA). Both targets were previously detected as unresolved point-sources by the SMA in its compact configuration, also at 345 GHz. These new data, which provide a factor of ~3 improvement in resolution, allow us to measure the physical scale of the far-infrared in the submillimetre directly. The visibility functions of both targets show significant evidence for structure on 0.5-1 arcsec scales, which at z=1.5 translates into a physical scale of 5-8 kpc. Our results are consistent with the angular and physical scales of two comparably luminous objects with high-resolution SMA followup, as well as radio continuum and CO sizes. These relatively compact sizes (<5-10 kpc) argue strongly for merger-driven starbursts, rather than extended gas-rich disks, as the preferred channel for forming SMGs. For the most luminous objects, the derived sizes may also have important physical consequences; under a series of simplifying assumptions, we find that these two objects in particular are forming stars close to or at the Eddington limit for a starburst.

Massive galaxies in the early Universe have been shown to be forming stars at surprisingly high rates. Prominent examples are dust-obscured galaxies which are luminous when observed at sub-millimeter (sub-mm) wavelengths and which may be forming stars at rates upto 1,000Mo/yr. These intense bursts of star formation are believed to be driven by mergers between gas rich galaxies. However, probing the properties of individual star-forming regions within these galaxies is beyond the spatial resolution and sensitivity of even the largest telescopes at present. Here, we report observations of the sub-mm galaxy SMMJ2135-0102 at redshift z=2.3259 which has been gravitationally magnified by a factor of 32 by a massive foreground galaxy cluster lens. This cosmic magnification, when combined with high-resolution sub-mm imaging, resolves the star-forming regions at a linear scale of just ~100 parsecs. We find that the luminosity densities of these star-forming regions are comparable to the dense cores of giant molecular clouds in the local Universe, but they are ~100x larger and 10^7 times more luminous. Although vigorously star-forming, the underlying physics of the star formation processes at z~2 appears to be similar to that seen in local galaxies even though the energetics are unlike anything found in the present-day Universe.

We use a suite of semi-empirical models to predict galaxy merger rates and contributions to bulge growth as functions of merger mass, redshift, and mass ratio. The models use empirical halo occupation constraints to identify mergers, together with high-resolution simulations to quantify how mergers with different properties contribute to the bulge population. We find good agreement with a variety of observational constraints, and provide fitting functions for merger rates and contributions to bulge growth. We identify several robust conclusions. (1) Major mergers dominate formation and assembly of L* bulges and the spheroid mass density, minor mergers contribute ~30%. (2) This is mass-dependent: bulge formation is dominated by more minor mergers in lower-mass systems. At higher masses, bulges form in major mergers near L*, but subsequently assemble in minor mergers. (3) The minor/major contribution is also morphology-dependent: higher B/T systems form in more major mergers, lower B/T systems form in situ from minor mergers. (4) Low-mass galaxies, being gas-rich, require more major mergers to reach the same B/T as high-mass systems. (5) Absolute merger rates increase with galaxy mass. (6) Predicted rates agree well with observations, but suggest that some morphology-selected samples include contamination from minor mergers. (7) Predicted rates agree with integrated growth in bulge mass with cosmic time, but with factor ~2 uncertainty - half the bulge mass density could come from non-mergers. We consider ~1000 model variations and quantify resulting uncertainties. Conclusions regarding the major/minor contribution to bulge growth are very robust, absolute merger rates have systematic factor ~2 uncertainties.

We use a novel method to predict the contribution of normal star-forming galaxies, merger-induced bursts, and obscured AGN, to IR luminosity functions (LFs) and global SFR densities. We use empirical halo occupation constraints to populate halos with galaxies and determine the distribution of normal and merging galaxies. Each system can then be associated with high-resolution hydrodynamic simulations. We predict the distribution of observed luminosities and SFRs, from different galaxy classes, as a function of redshift from z=0-6. We provide fitting functions for the predicted LFs, quantify the uncertainties, and compare with observations. At all redshifts, 'normal' galaxies dominate the LF at moderate luminosities ~L* (the 'knee'). Merger-induced bursts increasingly dominate at L>>L*; at the most extreme luminosities, AGN are important. However, all populations increase in luminosity at higher redshifts, owing to increasing gas fractions. Thus the 'transition' between normal and merger-dominated sources increases from the LIRG-ULIRG threshold at z~0 to bright Hyper-LIRG thresholds at z~2. The transition to dominance by obscured AGN evolves similarly, at factor of several higher L_IR. At all redshifts, non-merging systems dominate the total luminosity/SFR density, with merger-induced bursts constituting ~5-10% and AGN ~1-5%. Bursts contribute little to scatter in the SFR-stellar mass relation. In fact, many systems identified as 'ongoing' mergers will be forming stars in their 'normal' (non-burst) mode. Counting this as 'merger-induced' star formation leads to a stronger apparent redshift evolution in the contribution of mergers to the SFR density.

We describe a model for the formation of \zsim 2 Submillimeter Galaxies (SMGs) which simultaneously accounts for both average and bright SMGs while providing a reasonable match to their mean observed spectral energy distributions (SEDs). By coupling hydrodynamic simulations of galaxy mergers with the high resolution 3D polychromatic radiative transfer code Sunrise, we find that a mass sequence of merger models which use observational constraints as physical input naturally yield objects which exhibit black hole, bulge, and H2 gas masses similar to those observed in SMGs. The dominant drivers behind the 850 micron flux are the masses of the merging galaxies and the stellar birthcloud covering fraction. The most luminous (S850 ~ 15 mJy) sources are recovered by ~10^13 Msun 1:1 major mergers with a birthcloud covering fraction close to unity, whereas more average SMGs ~5-7 mJy) may be formed in lower mass halos ~5x10^12 Msun. These models demonstrate the need for high spatial resolution hydrodynamic and radiative transfer simulations in matching both the most luminous sources as well as the full SEDs of SMGs. While these models suggest a natural formation mechanism for SMGs, they do not attempt to match cosmological statistics of galaxy populations; future efforts along this line will help ascertain the robustness of these models.

We present results from a continuing interferometric survey of high-redshift submillimeter galaxies with the Submillimeter Array, including high-resolution (beam size ~2 arcsec) imaging of eight additional AzTEC 1.1mm selected sources in the COSMOS Field, for which we obtain six reliable (peak S/N>5 or peak S/N>4 with multiwavelength counterparts within the beam) and two moderate significance (peak S/N>4) detections. When combined with previous detections, this yields an unbiased sample of millimeter-selected SMGs with complete interferometric followup. With this sample in hand, we (1) empirically confirm the radio-submillimeter association, (2) examine the submillimeter morphology - including the nature of submillimeter galaxies with multiple radio counterparts and constraints on the physical scale of the far infrared - of the sample, and (3) find additional evidence for a population of extremely luminous, radio-dim submillimeter galaxies that peaks at higher redshift than previous, radio-selected samples. In particular, the presence of such a population of high-redshift sources has important consequences for models of galaxy formation - which struggle to account for such objects even under liberal assumptions - and dust production models given the limited time since the Big Bang.

We present a model for the CO molecular line emission from high redshift Submillimeter Galaxies (SMGs). By combining hydrodynamic simulations of gas rich galaxy mergers with the polychromatic radiative transfer code, Sunrise, and the 3D non-LTE molecular line radiative transfer code, Turtlebeach, we show that if SMGs are typically a transient phase of major mergers, their observed compact CO spatial extents, broad line widths, and high excitation conditions (CO SED) are naturally explained. In this sense, SMGs can be understood as scaled-up analogs to local ULIRGs. We utilize these models to investigate the usage of CO as an indicator of physical conditions. We find that care must be taken when applying standard techniques. The usage of CO line widths as a dynamical mass estimator from SMGs can possibly overestimate the true enclosed mass by a factor ~1.5-2. At the same time, assumptions of line ratios of unity from CO J=3-2 (and higher lying lines) to CO (J=1-0) will oftentimes lead to underestimates of the inferred gas mass. We provide tests for these models by outlining predictions for experiments which are imminently feasible with the current generation of bolometer arrays and radio-wave spectrometers.

We analyze a sample of galaxies chosen to have F24>0.5 mJy and satisfy a certain IRAC color criterion. IRS spectra yield redshifts, spectral types, and PAH luminosities, to which we add multi-wavelength broadband photometry. Stellar population modeling and IRS spectra together demonstrate that the double criteria used to select this sample have efficiently isolated massive star-forming galaxies at z~1.9. This is the first starburst-dominated ULIRG sample at high redshift with total infrared luminosity measured directly from FIR and millimeter photometry, and as such gives us the first accurate view of broadband SEDs for starburst galaxies at extremely high luminosity and at all wavelengths. Similar broadband data are assembled for three other galaxy samples -- local starburst galaxies, local AGN/ULIRGS, and a second 24mu-luminous z~1.9 sample dominated by AGN. L(PAH)/L(IR) for the new z~1.9 starburst sample is the highest ever seen, some three times higher than in local starbursts, whereas in AGNs this ratio is depressed below the starburst trend, often severely. Several pieces of evidence imply that AGNs exist in this starburst dominated sample, except two of which even host very strong AGN, while they still have very strong PAH emission. The ACS images show most objects have very extended morphologies in the rest-frame UV band, thus extended distribution of PAH molecules. Such an extended distribution prevents further destruction PAH molecules by central AGNs. We conclude that objects in this sample are ULIRGs powered mainly by starburst; and the total infrared luminosity density contributed by this type of objects is 0.9-2.6x 10^7 Lsol/Mpc^3.

We present a merger-driven evolutionary model for the production of luminous (LIRGs) and ultraluminous infrared galaxies (ULIRGs) with warm IR colours. Our results show that simulations of gas--rich major mergers including star formation, black hole growth, and feedback can produce warm (U)LIRGs. We also find that while the warm evolutionary phase is associated with increased AGN activity, star formation alone may be sufficient to produce warm IR colours. However, the transition can be suppressed entirely - even when there is a significant AGN contribution - when we assume a single-phase ISM, which maximizes the attenuation. Finally, our evolutionary models are consistent with the 25-to-60 micron flux density ratio vs. L_HX/L_IR relation for local LIRGs and ULIRGs, and predict the observed scatter in IR color at fixed L_HX/L_IR. Therefore, our models suggest a cautionary note in the interpretation of warm IR colours: while associated with periods of active black hole growth, they are probably produced by a complex mix of star formation and AGN activity intermediate between the cold star formation dominated phase and the birth of a bright, unobscured quasar.

We present millimetre observations of a sample of 12 high redshift ultraluminous infrared galaxies (ULIRGs) in the Extended Growth Strip (EGS). These objects were initially selected on the basis of their observed mid--IR colours (0.0 < [3.6]-[4.5] < 0.4 and -0.7 < [3.6]-[8.0] < 0.5) to lie at high redshift 1.5 < z < 3, and subsequent 20-38 micron mid-IR spectroscopy confirms that they lie in a narrow redshift window centered on z=2. We detect 9/12 of the objects in our sample at high significance (>3 sigma) with a mean 1200\micron flux of <F_1200> = 1.6+/-0.1 mJy. Our millimetre photometry, combined with existing far-IR photometry from the Far-IR Deep Extragalactic Legacy (FIDEL) Survey and accurate spectroscopic redshifts, places constraints both sides of the thermal dust peak. This allows us to estimate the dust properties, including the far--IR luminosity, dust temperature, and dust mass. We find that our sample is similar to other high-z and intermediate-z ULIRGs, and local systems, but has a different dust selection function than submillimeter-selected galaxies. Finally, we use existing 20cm radio continuum imaging to test the far-IR/radio correlation at high redshift. We find that our sample is consistent with the local relation, implying little evolution. Furthermore, this suggests that our sample selection method is efficient at identifying ultraluminous, starburst--dominated systems within a very narrow redshift range centered at z~2.

We develop a physical model for how galactic disks survive and/or are destroyed in interactions. Based on dynamical arguments, we show gas primarily loses angular momentum to internal torques in a merger. Gas within some characteristic radius (a function of the orbital parameters, mass ratio, and gas fraction of the merging galaxies), will quickly lose angular momentum to the stars sharing the perturbed disk, fall to the center and be consumed in a starburst. A similar analysis predicts where violent relaxation of the stellar disks is efficient. Our model allows us to predict the stellar and gas content that will survive to re-form a disk in the remnant, versus being violently relaxed or contributing to a starburst. We test this in hydrodynamic simulations and find good agreement as a function of mass ratio, orbital parameters, and gas fraction, in simulations spanning a wide range in these properties and others, including different prescriptions for gas physics and feedback. In an immediate sense, the amount of disk that re-forms can be understood in terms of well-understood gravitational physics, independent of details of ISM gas physics or feedback. This allows us to explicitly quantify the requirements for such feedback to (indirectly) enable disk survival, by changing the pre-merger gas content and distribution. The efficiency of disk destruction is a strong function of gas content: we show how and why sufficiently gas-rich major mergers can, under general conditions, yield systems with small bulges (B/T<0.2). We provide prescriptions for inclusion of our results in semi-analytic models.

Previous models of galactic disk heating in interactions invoke restrictive assumptions not necessarily valid in modern LCDM contexts: that satellites and orbits are rigid and circular, with slow decay over many orbital times from dynamical friction. This leads to a linear scaling of disk heating with satellite mass: disk heights and velocity dispersions scale ~M_sat/M_disk. In turn, observed disk thicknesses present strong constraints on merger histories: the implication for the Milky Way is that <5% of its mass could come from mergers since z~2, in conflict with cosmological predictions. More realistically, satellites merge on nearly radial orbits, and once near the disk, resonant interactions efficiently remove angular momentum while tidal effects strip mass, leading to rapid merger/destruction in a couple of free-fall plunges. Under these conditions the proper heating efficiency is non-linear in mass ratio, ~(M_sat/M_disk)^2. We derive the scaling of disk scale heights and velocity dispersions as a function of mass ratio and disk gas content in this regime, and show this accurately describes the results of simulations with proper 'live' halos and disks. Under realistic circumstances, disk heating in minor mergers is suppressed by an order of magnitude relative to expectations of previous models. We show that the Milky Way disk could have absorbed ~5-10 1:10 mass-ratio mergers since z=2, in agreement with cosmological models. These distinctions lead to dramatic differences in which mass ratios are most important for disk heating and in the isophotal shapes of disk+bulge systems.

We present high resolution submillimeter interferometric imaging of two of the brightest high-redshift submillimeter galaxies known: GN20 and AzTEC1 at 0.8 and 0.3 arcsec resolution respectively. Our data - the highest resolution submillimeter imaging of high redshift sources accomplished to date - was collected in three different array configurations: compact, extended, and very extended. We derive angular sizes of 0.6 and 1.0 arcsec for GN20 and 0.3 and 0.4 arcsec for AzTEC1 from modeling their visibility functions as a Gaussian and elliptical disk respectively. Because both sources are B-band dropouts, they likely lie within a relatively narrow redshift window around z~4, which indicates their angular extent corresponds to physical scales of 4-8 and 1.5-3 kpc respectively for the starburst region. By way of a series of simple assumptions, we find preliminary evidence that these hyperluminous starbursts - with star formation rates >1000 $M_\odot$ yr$^{-1}$ - are radiating at or close to their Eddington limit. Should future high resolution observations indicate that these two objects are typical of a population of high redshift Eddington-limited starbursts, this could have important consequences for models of star formation and feedback in extreme environments.

The outer regions of disk galaxies show a drop-off in optical and Halpha emission, suggesting a star formation threshold radius, assumed to owe to a critical surface density below which star formation does not take place. Signs of filamentary star formation beyond this threshold radius have been observed in individual galaxies in the Halpha and recent GALEX surveys have discovered that 30% of disk galaxies show UV emission out to 2-3 times the optical radius of the galaxy. We run smooth particle hydrodynamics simulations of disk galaxies with constant density extended gas disks to test whether over-densities owing to spiral structure in the outer disk can reproduce the observed star formation. We indeed find that spiral density waves from the inner disk propagate into the outer gas disk and raise local gas regions above the star formation density threshold, yielding features similar to those observed. Because the amount of star formation is low, we expect to see little optical emission in outer disks, as observed. Our results indicate that XUV disks can be simulated simply by adding an extended gas disk with a surface density near the threshold density to an isolated galaxy and evolving it with fiducial star formation parameters.

We present a series of simulations of the self--regulated growth of supermassive black holes (SMBHs) in galaxies via three different fueling mechanisms: major mergers, minor mergers, and disk instabilities. The SMBHs in all three scenarios follow the same black hole fundamental plane (BHFP) and correlation with bulge binding energy seen in simulations of major mergers, and observed locally. Furthermore, provided that the total gas supply is significantly larger than the mass of the SMBH, its limiting mass is not influenced by the amount of gas available or the efficiency of black hole growth. This supports the assertion that SMBHs accrete until they reach a critical mass at which feedback is sufficient to unbind the gas locally, terminating the inflow and stalling further growth. At the same time, while minor and major mergers follow the same projected correlations (e.g., the $M_{BH}-\sigma$ and Magorrian relations), SMBHs grown via disk instabilities do not, owing to structural differences between the host bulges. This finding is supported by recent observations of SMBHs in pseudobulges and bulges in barred systems, as compared to those hosted by classical bulges. Taken together, this provides support for the BHFP and binding energy correlations as being more "fundamental" than other proposed correlations in that they reflect the physical mechanism driving the co-evolution of SMBHs and spheroids.

High-redshift submillimetre-bright galaxies identified by blank field surveys at millimetre and submillimetre wavelengths appear in the region of the IRAC colour-colour diagrams previously identified as the domain of luminous active galactic nuclei (AGNs). Our analysis using a set of empirical and theoretical dusty starburst spectral energy distribution (SED) models shows that power-law continuum sources associated with hot dust heated by young (<100 Myr old), extreme starbursts at z>2 also occupy the same general area as AGNs in the IRAC colour-colour plots. A detailed comparison of the IRAC colours and SEDs demonstrates that the two populations are distinct from each other, with submillimetre-bright galaxies having a systematically flatter IRAC spectrum (>1 mag bluer in the observed [4.5]-[8.0] colour). Only about 20% of the objects overlap in the colour-colour plots, and this low fraction suggests that submillimetre galaxies powered by a dust-obscured AGN are not common. The red IR colours of the submillimetre galaxies are distinct from those of the ubiquitous foreground IRAC sources, and we propose a set of IR colour selection criteria for identifying SMG counterparts that can be used even in the absence of radio or Spitzer MIPS 24 micron data.

We present high-resolution interferometric imaging of LH850.02, the brightest 850- and 1200-micron submillimetre (submm) galaxy in the Lockman Hole. Our observations were made at 890 micron with the Submillimetre Array (SMA). Our high-resolution submm imaging detects LH850.02 at >6-sigma as a single compact (size < 1 arcsec or < 8 kpc) point source and yields its absolute position to ~0.2-arcsec accuracy. LH850.02 has two alternative radio counterparts within the SCUBA beam (LH850.02N & S), both of which are statistically very unlikely to be so close to the SCUBA source position by chance. However, the precise astrometry from the SMA shows that the submm emission arises entirely from LH850.02N, and is not associated with LH850.02S (by far the brighter of the two alternative identifications at 24-micron). Fits to the optical-infrared multi-colour photometry of LH850.02N & S indicate that both lie at z~3.3, and are therefore likely to be physically associated. At these redshifts, the 24 micron--to--submm flux density ratios suggest that LH850.02N has an Arp220-type starburst-dominated far-IR SED, while LH850.02S is more similar to Mrk231, with less dust-enshrouded star-formation activity, but a significant contribution at 24-micron (rest-frame ~5-6 micron) from an active nucleus. This complex mix of star-formation and AGN activity in multi-component sources may be common in the high redshift ultraluminous galaxy population, and highlights the need for precise astrometry from high resolution interferometric imaging for a more complete understanding.

We present a scenario for the production of dynamically cold rings around the Milky Way via a high-eccentricity, flyby encounter. These initial conditions are more cosmologically motivated than those considered in previous works. We find that the encounters we examine generically produce a series of nearly dynamically cold ring-like features on low-eccentricity orbits that persist over timescales of ~2-4 Gyr via the tidal response of the primary galaxy to the close passage of the satellite. Moreover, they are both qualitatively and quantitatively similar to the distribution, kinematics, and stellar population of the Monoceros ring. Therefore, we find that a high eccentricity flyby by a satellite galaxy represents a cosmologically appealing scenario for forming kinematically distinct ring-like features around the Milky Way.

We have used the Submillimeter Array to image a flux limited sample of seven submillimeter galaxies, selected by the AzTEC camera on the JCMT at 1.1 mm, in the COSMOS field at 890um with 2" resolution. All of the sources - two radio-bright and five radio-dim - are detected as single point-sources at high significance (> 6\sigma), with positions accurate to 0.2" that enable counterpart identification at other wavelengths observed with similarly high angular resolution. All seven have IRAC counterparts, but only two have secure counterparts in deep HST/ACS imaging. As compared to the two radio-bright sources in the sample, and those in previous studies, the five radio-dim sources in the sample (1) have systematically higher submillimeter-to-radio flux ratios, (2) have lower IRAC 3.6-8.0um fluxes, and (3) are not detected at 24um. These properties, combined with size constraints at 890um (\theta < 1.2"), suggest that the radio-dim submillimeter galaxies represent a population of very dusty starbursts, with physical scales similar to local ultraluminous infrared galaxies, and an average redshift higher than radio-bright sources.

We present the rest-frame ultraviolet through near infrared spectral energy distribution for an interacting Lyman break galaxy at a redshift z=4.42, the highest redshift merging system known with clearly resolved tidal features. The two objects in this system - HDF-G4 and its previously unidentified companion - are both B_{435} band dropouts, have similar V_{606}-i_{775} and i_{775}-z_{850} colors, and are separated by 1", which at z=4.42 corresponds to 7 kpc projected nuclear separation; all indicative of an interacting system. Fits to stellar population models indicate a stellar mass of M_\star = 2.6\times 10^{10} M_\odot, age of \tau_\star = 720 My, and exponential star formation history with an e-folding time \tau_0 = 440 My. Using these derived stellar populations as constraints, we model the HDF-G4 system using hydrodynamical simulations, and find that it will likely evolve into a quasar by z\sim3.5, and a quiescent, compact spheroid by z\sim 2.5 similar to those observed at z > 2. And, the existence of such an object supports galaxy formation models in which major mergers drive the high redshift buildup of spheroids and black holes.

We use hydrodynamic simulations of minor mergers of galaxies to investigate the nature of surface brightness excesses at large radii observed in some spiral galaxies: antitruncated stellar disks. We find that this process can produce the antitruncation via two competing effects: (1) merger-driven gas inflows that concentrate mass in the center of the primary galaxy and contract its inner density profile; and (2) angular momentum transferred outwards by the interaction, causing the outer disk to expand. In our experiments, this requires both a significant supply of gas in the primary disk, and that the encounter be prograde with moderate orbital angular momentum. The stellar surface mass density profiles of our remnants both qualitatively and quantitatively resemble the broken exponentials observed in local face--on spirals that display antitruncations. Moreover, the observed trend towards more frequent antitruncation relative to classical truncation in earlier Hubble types is consistent with a merger-driven scenario.

We perform a set of non--radiative cosmological simulations of a preheated intracluster medium in which the entropy of the gas was uniformly boosted at high redshift. The results of these simulations are used first to test the current analytic techniques of preheating via entropy input in the smooth accretion limit. When the unmodified profile is taken directly from simulations, we find that this model is in excellent agreement with the results of our simulations. This suggests that preheated efficiently smoothes the accreted gas, and therefore a shift in the unmodified profile is a good approximation even with a realistic accretion history. When we examine the simulation results in detail, we do not find strong evidence for entropy amplification, at least for the high-redshift preheating model adopted here. In the second section of the paper, we compare the results of the preheating simulations to recent observations. We show -- in agreement with previous work -- that for a reasonable amount of preheating, a satisfactory match can be found to the mass-temperature and luminosity-temperature relations. However -- as noted by previous authors -- we find that the entropy profiles of the simulated groups are much too flat compared to observations. In particular, while rich clusters converge on the adiabatic self--similar scaling at large radius, no single value of the entropy input during preheating can simultaneously reproduce both the core and outer entropy levels. As a result, we confirm that the simple preheating scenario for galaxy cluster formation, in which entropy is injected universally at high redshift, is inconsistent with observations.

Forthcoming large galaxy cluster surveys will yield tight constraints on cosmological models. It has been shown that in an idealized survey, containing > 10,000 clusters, statistical errors on dark energy and other cosmological parameters will be at the percent level. It has also been shown that through "self-calibration", parameters describing the mass-observable relation and cosmology can be simultaneously determined, though at a loss in accuracy by about an order of magnitude. Here we examine the utility of an alternative approach of self-calibration, in which a parametrized ab-initio physical model is used to compute cluster structure and the resulting mass-observable relations. As an example, we use a modified-entropy ("pre-heating") model of the intracluster medium, with the history and magnitude of entropy injection as unknown input parameters. Using a Fisher matrix approach, we evaluate the expected simultaneous statistical errors on cosmological and cluster model parameters. We study two types of surveys, in which a comparable number of clusters are identified either through their X-ray emission or through their integrated Sunyaev-Zel'dovich (SZ) effect. We find that compared to a phenomenological parametrization of the mass-observable relation, using our physical model yields significantly tighter constraints in both surveys, and offers substantially improved synergy when the two surveys are combined. These results suggest that parametrized physical models of cluster structure will be useful when extracting cosmological constraints from SZ and X-ray cluster surveys. (abridged)

We have discovered 21 new Type Ia supernovae (SNe Ia) with the Hubble Space Telescope (HST) and have used them to trace the history of cosmic expansion over the last 10 billion years. These objects, which include 13 spectroscopically confirmed SNe Ia at z > 1, were discovered during 14 epochs of reimaging of the GOODS fields North and South over two years with the Advanced Camera for Surveys on HST. Together with a recalibration of our previous HST-discovered SNe Ia, the full sample of 23 SNe Ia at z > 1 provides the highest-redshift sample known. Combined with previous SN Ia datasets, we measured H(z) at discrete, uncorrelated epochs, reducing the uncertainty of H(z>1) from 50% to under 20%, strengthening the evidence for a cosmic jerk--the transition from deceleration in the past to acceleration in the present. The unique leverage of the HST high-redshift SNe Ia provides the first meaningful constraint on the dark energy equation-of-state parameter at z >1. The result remains consistent with a cosmological constant (w(z)=-1), and rules out rapidly evolving dark energy (dw/dz >>1). The defining property of dark energy, its negative pressure, appears to be present at z>1, in the epoch preceding acceleration, with ~98% confidence in our primary fit. Moreover, the z>1 sample-averaged spectral energy distribution is consistent with that of the typical SN Ia over the last 10 Gyr, indicating that any spectral evolution of the properties of SNe Ia with redshift is still below our detection threshold.

The evolution of the cluster mass function and the cluster correlation function from z = 0 to z = 3 are determined using 10^6 clusters obtained from high-resolution simulations of the current best-fit LCDM cosmology (\Omega_m = 0.27, \sigma_8 = 0.84, h = 0.7). The results provide predictions for comparisons with future observations of high redshift clusters. A comparison of the predicted mass function of low redshift clusters with observations from early Sloan Digital Sky Survey data, and the predicted abundance of massive distant clusters with observational results, favor a slightly larger amplitude of mass fluctuations (\sigma_8 = 0.9) and lower density parameter (\Omega_m = 0.2); these values are consistent within 1-\sigma with the current observational and model uncertainties. The cluster correlation function strength increases with redshift for a given mass limit; the clusters were more strongly correlated in the past, due to their increasing bias with redshift - the bias reaches b = 100 at z = 2 for M > 5 x 10^13 h^-1 M_sun. The richness-dependent cluster correlation function, represented by the correlation scale versus cluster mean separation relation, R0-d, is generally consistent with observations. This relation can be approximated as R_0 = 1.7 d^0.6 h^-1 Mpc for d = 20 - 60 h^-1 Mpc. The R0-d relation exhibits surprisingly little evolution with redshift for z < 2; this can provide a new test of the current LCDM model when compared with future observations of high redshift clusters.

We present an empirical algorithm for obtaining photometric redshifts of quasars using 5-band Sloan Digital Sky Survey (SDSS) photometry. Our algorithm generates an empirical model of the quasar color-redshift relation, compares the colors of a quasar candidate with this model, and calculates possible photometric redshifts. Using the 3814 quasars of the SDSS Early Data Release Quasar Catalog to generate a median color-redshift relation as a function of redshift we find that, for this same sample, 83% of our predicted redshifts are correct to within |Delta z|<0.3. The algorithm also determines the probability that the redshift is correct, allowing for even more robust photometric redshift determination for smaller, more restricted samples. We apply this technique to a set of 8740 quasar candidates selected by the final version of the SDSS quasar-selection algorithm. The photometric redshifts assigned to non-quasars are restricted to a few well-defined values. In addition, 90% of the objects with spectra that have photometric redshifts between 0.8 and 2.2 are quasars with accurate (|Delta z|<0.3) photometric redshifts. Many of these quasars lie in a single region of color space; judicious application of color-cuts can effectively select quasars with accurate photometric redshifts from the SDSS database. When the SDSS is complete, this technique will allow the determination of photometric redshifts for ~10^6 faint SDSS quasar candidates, enabling advances in our knowledge of the quasar luminosity function, gravitational lensing of quasars, and correlations among quasars and between galaxies.