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

We use simulations with realistic models for stellar feedback to study galaxy mergers. These high resolution (1 pc) simulations follow formation and destruction of individual GMCs and star clusters. The final starburst is dominated by in situ star formation, fueled by gas which flows inwards due to global torques. The resulting high gas density results in rapid star formation. The gas is self gravitating, and forms massive (~10^10 M_sun) GMCs and subsequent super-starclusters (masses up to 10^8 M_sun). However, in contrast to some recent simulations, the bulk of new stars which eventually form the central bulge are not born in superclusters which then sink to the center of the galaxy, because feedback efficiently disperses GMCs after they turn several percent of their mass into stars. Most of the mass that reaches the nucleus does so in the form of gas. The Kennicutt-Schmidt law emerges naturally as a consequence of feedback balancing gravitational collapse, independent of the small-scale star formation microphysics. The same mechanisms that drive this relation in isolated galaxies, in particular radiation pressure from IR photons, extend over seven decades in SFR to regulate star formation in the most extreme starbursts (densities >10^4 M_sun/pc^2). Feedback also drives super-winds with large mass loss rates; but a significant fraction of the wind material falls back onto the disks at later times, leading to higher post-starburst SFRs in the presence of stellar feedback. Strong AGN feedback is required to explain sharp cutoffs in star formation rate. We compare the predicted relic structure, mass profile, morphology, and efficiency of disk survival to simulations which do not explicitly resolve GMCs or feedback. Global galaxy properties are similar, but sub-galactic properties and star formation rates can differ significantly.

We analyze the link between active galactic nuclei (AGN) and mid-infrared flux using dust radiative transfer calculations of starbursts realized in hydrodynamical simulations. Focusing on the effects of galaxy dust, we evaluate diagnostics commonly used to disentangle AGN and star formation in ultraluminous infrared galaxies (ULIRGs). We examine these quantities as a function of time, viewing angle, dust model, AGN spectrum, and AGN strength in merger simulations representing two possible extremes of the ULIRG population: one is a typical gas-rich merger at z ~ 0, and the other is characteristic of extremely obscured starbursts at z ~ 2 to 4. This highly obscured burst begins star-formation-dominated with significant PAH emission, and ends with a ~10^9 yr period of red near-IR colors. At coalescence, when the AGN is most luminous, dust obscures the near-infrared AGN signature, reduces the relative emission from polycyclic aromatic hydrocarbons (PAHs), and enhances the 9.7 micron absorption by silicate grains. Although generally consistent with previous interpretations, our results imply none of these indicators can unambiguously estimate the AGN luminosity fraction in all cases. Motivated by the simulations, we show that a combination of the extinction feature at 9.7 micron, the PAH strength, and a near-infrared slope can simultaneously constrain the AGN fraction and dust grain distribution for a wide range of obscuration. We find that this indicator, accessible to the James Webb Space Telescope, may estimate the AGN power as tightly as the hard X-ray flux alone, thereby providing a valuable future cross-check and constraint for large samples of distant ULIRGs.

We present an analysis of the spatial distribution of star formation in a sample of 60 visually identified galaxy merger candidates at z>1. Our sample, drawn from the 3D-HST survey, is flux-limited and was selected to have high star formation rates based on fits of their broad-band, low spatial resolution spectral energy distributions. It includes plausible pre-merger (close pairs) and post-merger (single objects with tidal features) systems, with total stellar masses and star formation rates derived from multi-wavelength photometry. Here we use near-infrared slitless spectra from 3D-HST which produce Halpha or [OIII] emission line maps as proxies for star-formation maps. This provides a first comprehensive high-resolution, empirical picture of where star formation occurred in galaxy mergers at the epoch of peak cosmic star formation rate. We find that detectable star formation can occur in one or both galaxy centres, or in tidal tails. The most common case (58%) is that star formation is largely concentrated in a single, compact region, coincident with the centre of (one of) the merger components. No correlations between star formation morphology and redshift, total stellar mass, or star formation rate are found. A restricted set of hydrodynamical merger simulations between similarly massive and gas-rich objects implies that star formation should be detectable in both merger components, when the gas fractions of the individual components are the same. This suggests that z~1.5 mergers typically occur between galaxies whose gas fractions, masses, and/or star formation rates are distinctly different from one another.

(Abridged) We study relationships between the SFR and the nuclear properties of X-ray selected AGNs out to z=2.5, using far-IR data in three extragalactic deep fields as part of the PACS Evolutionary Probe (PEP) program. Guided by studies of intrinsic infra-red AGN SEDs, we show that the majority of the FIR emission in AGNs is produced by cold dust heated by star-formation. We uncover characteristic redshift-dependent trends between the mean FIR luminosity (L_fir) and accretion luminosity (L_agn) of AGNs. At low AGN luminosities, accretion and SFR are uncorrelated at all redshifts, consistent with a scenario where most low-luminosity AGNs are primarily fueled by secular processes in their host galaxies. At high AGN luminosities, a significant correlation is observed between L_fir and L_agn, but only among AGNs at low and moderate redshifts (z<1). We interpret this as a signature of the increasing importance of major-mergers in driving both the growth of super-massive black holes (SMBHs) and global star-formation in their hosts at high AGN luminosities. However, we also find that the enhancement of SFR in luminous AGNs weakens or disappears at high redshifts (z>1). This suggests that the role of mergers in SMBH-galaxy co-evolution is less important at these epochs. At all redshifts, we find essentially no relationship between L_fir and nuclear obscuration across five orders of magnitude in obscuring column density, suggesting that various different mechanisms are likely to be responsible for obscuring X-rays in active galaxies. We explain our results within a scenario in which two different modes of SMBH fueling operate among low- and high-luminosity AGNs. We postulate, guided by emerging knowledge about the properties of high redshift galaxies, that the dominant mode of accretion among high-luminosity AGNs evolves with redshift.

High-redshift submillimetre galaxies (SMGs) are some of the most rapidly star-forming galaxies in the Universe. Historically, galaxy formation models have had difficulty explaining the observed number counts of SMGs. We combine a semi-empirical model with 3-D hydrodynamical simulations and 3-D dust radiative transfer to predict the number counts of unlensed SMGs. Because the stellar mass functions, gas and dust masses, and sizes of our galaxies are constrained to match observations, we can isolate uncertainties related to the dynamical evolution of galaxy mergers and the dust radiative transfer. The number counts and redshift distributions predicted by our model agree well with observations. Isolated disc galaxies dominate the faint (S_1.1 < ~1 mJy, or S_850 < ~2 mJy) population. The brighter sources are a mix of merger-induced starbursts and galaxy-pair SMGs; the latter subpopulation accounts for ~30 - 50 per cent of all SMGs at all S_1.1 > ~0.5 mJy (S_850 > ~1 mJy). The mean redshifts are ~3.0 - 3.5, depending on the flux cut, and the brightest sources tend to be at higher redshifts. Because the galaxy-pair SMGs will be resolved into multiple fainter sources by ALMA, the bright ALMA counts should be as much as 2 times less than those observed using single-dish telescopes. The agreement between our model, which uses a Kroupa IMF, and observations suggests that the IMF in high-redshifts starbursts need not be top-heavy; if the IMF were top-heavy, our model would over-predict the number counts. We conclude that the difficulty some models have reproducing the observed SMG counts is likely indicative of more general problems -- such as an under-prediction of the abundance of massive galaxies or a star formation rate--stellar mass relation normalisation lower than that observed -- rather than a problem specific to the SMG population.

We study the future orbital evolution and merging of the MW-M31-M33 system, using a combination of collisionless N-body simulations and semi-analytic orbit integrations. Monte-Carlo simulations are used to explore the consequences of varying the initial phase-space and mass parameters within their observational uncertainties. The observed M31 transverse velocity implies that the MW and M31 will merge t = 5.86 (+1.61-0.72) Gyr from now, after a first pericenter at t = 3.87 (+0.42-0.32) Gyr. M31 may (probability p=41%) make a direct hit with the MW (defined here as a first-pericenter distance less than 25 kpc). Most likely, the MW and M31 will merge first, with M33 settling onto an orbit around them. Alternatively, M33 may make a direct hit with the MW first (p=9%), or M33 may get ejected from the Local Group (p=7%). The MW-M31 merger remnant will resemble an elliptical galaxy. The Sun will most likely (p=85%) end up at larger radius from the center of the MW-M31 merger remnant than its current distance from the MW center, possibly further than 50 kpc (p=10%). The Sun may (p=20%) at some time in the next 10 Gyr find itself moving through M33 (within 10 kpc), but while dynamically still bound to the MW-M31 merger remnant. The arrival and possible collision of M31 (and possibly M33) with the MW is the next major cosmic event affecting the environment of our Sun and solar system that can be predicted with some certainty. (Abridged)

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

In recent work (arXiv:1101.0002) we have suggested that the high-redshift (z ~ 2-4) bright submillimetre galaxy (SMG) population is heterogeneous, with major mergers contributing both at early stages, where quiescently star-forming discs are blended into one submm source (`galaxy-pair SMGs'), and late stages, where mutual tidal torques drive gas inflows and cause strong starbursts. Here we combine hydrodynamic simulations of major mergers with 3-D dust radiative transfer calculations to determine observational diagnostics that can distinguish between quiescently star-forming SMGs and starburst SMGs via integrated data alone. We fit the far-IR SEDs of the simulated galaxies with the optically thin single-temperature modified blackbody, the full form of the single-temperature modified blackbody, and a power-law temperature-distribution model. The effective dust temperature, T_dust, and power-law index of the dust emissivity in the far-IR, \beta, derived can significantly depend on the fitting form used, and the intrinsic \beta\ of the dust is not recovered. However, for all forms used here, there is a T_dust above which almost all simulated galaxies are starbursts, so a T_dust cut is very effective at selecting starbursts. Simulated merger-induced starbursts also have higher L_IR/M_gas and L_IR/L_FUV than quiescently star-forming galaxies and lie above the star formation rate-stellar mass relation. These diagnostics can be used to test our claim that the SMG population is heterogeneous and to observationally determine what star formation mode dominates a given galaxy population. We comment on applicability of these diagnostics to ULIRGs that would not be selected as SMGs. These `hot-dust ULIRGs' are typically starburst galaxies lower in mass than SMGs, but they can also simply be SMGs observed from a different viewing angle.

We present a novel pair of numerical models of the interaction history between the Large and Small Magellanic Clouds (LMC and SMC, respectively) and our Milky Way (MW) in light of recent high precision proper motions (Kallivayalil et al. 2006a,b). Given the new velocities, cosmological simulations of structure formation favor a scenario where the Magellanic Clouds (MCs) are currently on their first infall towards our Galaxy (Boylan-Kolchin et al. 2011, Busha et al. 2011). We illustrate here that the observed irregular morphology and internal kinematics of the MCs (in gas and stars) are naturally explained by interactions between the LMC and SMC, rather than gravitational interactions with the MW. This picture further supports a first infall scenario (Besla et a. 2007). In particular, we demonstrate that the Magellanic Stream, a band of HI gas trailing behind the MCs 150 degrees across the sky, can be accounted for by the action of LMC tides on the SMC before the system was accreted by the MW. We further demonstrate that the off-center, warped stellar bar of the LMC and its one-armed spiral, can be naturally explained by a recent direct collision with the SMC. Such structures are key morphological characteristics of a class of galaxies referred to as Magellanic Irregulars (de Vaucouleurs & Freeman 1972), the majority of which are not associated with massive spiral galaxies. We infer that dwarf-dwarf galaxy interactions are important drivers for the morphological evolution of Magellanic Irregulars and can dramatically affect the efficiency of baryon removal from dwarf galaxies via the formation of extended tidal bridges and tails. Such interactions are important not only for the evolution of dwarf galaxies but also have direct consequences for the buildup of baryons in our own MW, as LMC-mass systems are believed to be the dominant building blocks of MW-type halos.

Models of poststarburst (or "K+A") galaxies are constructed by combining fully three-dimensional hydrodynamic simulations of galaxy mergers with radiative transfer calculations of dust attenuation. Spectral line catalogs are generated automatically from moderate-resolution optical spectra calculated as a function of merger progress in each of a large suite of simulations. The mass, gas fraction, orbital parameters, and mass ratio of the merging galaxies are varied systematically, showing that the lifetime and properties of the K+A phase are strong functions of merger scenario. K+A durations are generally less than ~0.1-0.3 Gyr, significantly shorter than the commonly assumed 1 Gyr, which is obtained only in rare cases, owing to a wide variation in star formation histories resulting from different orbital and progenitor configurations. Combined with empirical merger rates, the model lifetimes predict rapidly-rising K+A fractions as a function of redshift that are consistent with results of large spectroscopic surveys, resolving tension between the observed K+A abundance and that predicted when one assumes the K+A duration is the lifetime of A stars (~1 Gyr). The effects of dust attenuation, viewing angle, and aperture bias on our models are analyzed. In some cases, the K+A features are longer-lived and more pronounced when AGN feedback removes dust from the center, uncovering the young stars formed during the burst. In this picture, the K+A phase begins during or shortly after the bright starburst/AGN phase in violent mergers, and thus offers a unique opportunity to study the effects of quasar and star formation feedback on the gas reservoir and evolution of the remnant. Analytic fitting formulae are provided for the estimates of K+A incidence as a function of merger scenario.

Central supermassive black holes (SMBHs) are a ubiquitous feature of locally-observed galaxies, and ample evidence suggests that the growth of SMBHs and their host galaxies is closely linked. However, in the event of a merger, gravitational-wave (GW) recoil may displace a SMBH from its galactic center, or eject it entirely. To explore the consequences of this phenomenon, we use hydrodynamic simulations of gaseous galaxy mergers that include a range of BH recoil velocities. We have generated a suite of over 200 simulations with more than 60 merger models, enabling us to identify systematic trends in the behavior of recoiling BHs -- specifically (i) their dynamics, (ii) their observable signatures, and (iii) their effects on BH/galaxy co-evolution. (i) Recoiling BH trajectories depend heavily on the gas content of the host galaxy; maximal BH displacements from the center may vary by up to an order of magnitude between gas-rich and gas-poor mergers. In some cases, recoil trajectories also depend on the timing of the BH merger relative to the formation of the galaxy merger remnant. (ii) Recoiling BHs may be observable as offset active galactic nuclei (AGN) via either kinematic offsets (v > 800 km s^-1) or spatial offsets (R > 1 kpc) for lifetimes of about 1 - 100 Myr. In addition, recoil events affect the total AGN lifetime. GW recoil generally reduces the lifetimes of bright AGN, but may extend lower-luminosity AGN lifetimes. (iii) Rapidly-recoiling BHs may be up to about 5 times less massive than their stationary counterparts. These mass deficits lower the normalization of the M - sigma relation and contribute to both intrinsic and overall scatter. Furthermore, recoil events displace AGN feedback from the galactic center, which enhances central star formation rates. This results in longer starburst phases and higher central stellar densities in merger remnants.

We perform 3-D dust radiative transfer (RT) calculations on hydrodynamic simulations of isolated and merging disk galaxies in order to quantitatively study the dependence of observed-frame submillimeter (submm) flux density on galaxy properties. We find that submm flux density and star formation rate (SFR) are related in dramatically different ways for quiescently star-forming galaxies and starbursts. Because the stars formed in the merger-induced starburst do not dominate the bolometric luminosity and the rapid drop in dust mass and more compact geometry cause a sharp increase in dust temperature during the burst, starbursts are very inefficient at boosting submm flux density (e.g., a $\ga16$x boost in SFR yields a $\la 2$x boost in submm flux density). Moreover, the ratio of submm flux density to SFR differs significantly between the two modes; thus one cannot assume that the galaxies with highest submm flux density are necessarily those with the highest bolometric luminosity or SFR. These results have important consequences for the bright submillimeter-selected galaxy (SMG) population. Among them are: 1. The SMG population is heterogeneous. In addition to merger-driven starbursts, there is a subpopulation of galaxy pairs, where two disks undergoing a major merger but not yet strongly interacting are blended into one submm source because of the large ($\ga 15$", or $\sim 130$ kpc at $z = 2$) beam of single-dish submm telescopes. 2. SMGs must be very massive ($M_{\star} \ga 6 \times 10^{10} \msun$). 3. The infall phase makes the SMG duty cycle a factor of a few greater than what is expected for a merger-driven starburst. (Abridged.)

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

We employ hydrodynamical simulations to study the effects of dissipational gas physics on the vertical heating and thickening of disc galaxies during minor mergers. For the first time we present a suite of simulations that includes a diffuse, rotating, cooling, hot gaseous halo, as predicted by cosmological hydrodynamical simulations as well as models of galaxy formation. We study the effect of this new gaseous component on the vertical structure of a Milky Way-like stellar disc during 1:10 and 1:5 mergers. For 1:10 mergers we find no increased final thin disc scale height compared to the isolated simulation, leading to the conclusion that thin discs can be present even after a 1:10 merger if a reasonable amount of hot gas is present. The reason for this is the accretion of new cold gas, leading to the formation of a massive new thin stellar disc that dominates the surface brightness profile. In a previous study, in which we included only cold gas in the disk, we showed that the presence of cold gas decreased the thickening by a minor merger relative to the no-gas case. Here, we show that the evolution of the scale height in the presences of a cooling hot halo is dominated by the formation of the new stellar disc. In this scenario, the thick disc is the old stellar disc that has been thickened in a minor merger at z>1, while the thin disc is the new stellar disc that reforms after this merger. In addition, we study the evolution of the scale height during a 1:5 merger and find that a thin disc can be present even after this merger, provided enough hot gas is available. The final scale height in our simulations depends on the mass of the hot gaseous halo, the efficiency of the winds and the merger mass ratio. We find post-merger values in the range 0.5<z0<1.0 kpc in good agreement with observational constraints by local galaxies.

Nuclear inflows of metal-poor interstellar gas triggered by galaxy interactions can account for the systematically lower central oxygen abundances observed in local interacting galaxies. Here, we investigate the metallicity evolution of a large set of simulations of colliding galaxies. Our models include cooling, star formation, feedback, and a new stochastic method for tracking the mass recycled back to the interstellar medium from stellar winds and supernovae. We study the influence of merger-induced inflows, enrichment, gas consumption, and galactic winds in determining the nuclear metallicity. The central metallicity is primarily a competition between the inflow of low-metallicity gas and enrichment from star formation. An average depression in the nuclear metallicity of ~0.07 is found for gas-poor disk-disk interactions. Gas-rich disk-disk interactions, on the other hand, typically have an enhancement in the central metallicity that is positively correlated with the gas content. The simulations fare reasonably well when compared to the observed mass-metallicity and separation-metallicity relationships, but further study is warranted.

Recent observations of excited CO emission lines from z~2 disc galaxies have shed light on the Kennicutt-Schmidt relation at high-z via observed SFR-CO (J=2-1) and (J=3-2) relations. Here, we describe a novel methodology for utilising these observations of high-excitation CO to derive the underlying Schmidt (SFR-rho^N) relationship. To do this requires an understanding of the potential effects of differential CO excitation with SFR. If the most heavily star-forming galaxies have a larger fraction of their gas in highly excited CO states than the lower SFR galaxies, then the observed molecular SFR-CO^alpha index, alpha, will be less than the underlying (volumetric) Schmidt index, N. Utilising a combination of SPH models of galaxy evolution and molecular line radiative transfer, we present the first calculations of CO excitation in z~2 disc galaxies with the aim of developing a mapping between various observed SFR-CO relationships and the underlying Schmidt relation. We find that even in relatively luminous z~2 discs, differential excitation does indeed exist, resulting in alpha < N for highly excited CO lines. This means that an observed (e.g.) SFR-CO (J=3-2) relation does not map linearly to SFR-H2 relation. We utilise our model results to provide a mapping from alpha to N for the range of Schmidt indices N=1-2. By comparing to recent observational surveys, we find that the observed SFR-CO (J=2-1) and SFR-CO (J=3-2) relations suggest that an underlying SFR rho^1.5 relation describes z~2 disc galaxies.

Major mergers of disk galaxies are thought to be a substantial driver in galaxy evolution. To trace the fraction and the rate galaxies are in mergers over cosmic times, several observational techniques, including morphological selection criteria, have been developed over the last decade. We apply this morphological selection of mergers to 21 cm radio emission line (HI) column density images of spiral galaxies in nearby surveys. In this paper, we investigate how long a 1:1 merger is visible in HI from N-body simulations. We evaluate the merger visibility times for selection criteria based on four parameters: Concentration, Asymmetry, M20, and the Gini parameter of second order moment of the flux distribution (GM). Of three selection criteria used in the literature, one based on Concentration and M20 works well for the HI perspective with a merger time scale of 0.4 Gyr. Of the three selection criteria defined in our previous paper, the GM performs well and cleanly selects mergers for 0.69 Gyr. The other two criteria (A-M20 and C-M20), select isolated disks as well, but perform best for face-on, gas-rich disks (T(merger) ~ 1 Gyr). The different visibility scales can be combined with the selected fractions of galaxies in any large HI survey to obtain merger rates in the nearby Universe. All-sky surveys such as WALLABY with ASKAP and the Medium Deep Survey with the APETIF instrument on Westerbork are set to revolutionize our perspective on neutral hydrogen and will provide an accurate measure of the merger fraction and rate of the present epoch.

Cosmological hydrodynamical simulations as well as observations indicate that spiral galaxies are comprised of five different components: dark matter halo, stellar disc, stellar bulge, gaseous disc and gaseous halo. While the first four components have been extensively considered in numerical simulations of binary galaxy mergers, the effect of a hot gaseous halo has usually been neglected even though it can contain up to 80% of the total gas within the galaxy virial radius. We present a series of hydrodynamic simulations of major mergers of disc galaxies, that for the first time include a diffuse, rotating, hot gaseous halo. Through cooling and accretion, the hot halo can dissipate and refuel the cold gas disc before and after a merger. This cold gas can subsequently form stars, thus impacting the morphology and kinematics of the remnant. Simulations of isolated systems with total mass M~10^12Msun show a nearly constant star formation rate of ~5Msun/yr if the hot gaseous halo is included, while the star formation rate declines exponentially if it is neglected. We conduct a detailed study of the star formation efficiency during mergers and find that the presence of a hot gaseous halo reduces the starburst efficiency (e=0.5) compared to simulations without a hot halo (e=0.68). Moreover we find cases where the stellar mass of the merger remnant is lower than the sum of the stellar mass of the two progenitor galaxies when evolved in isolation. This suggests a revision to semi-analytic galaxy formation models which assume that a merger always leads to enhanced star formation. We show that adding the hot gas component has a significant effect on the kinematics and internal structure of the merger remnants, like an increased abundance of fast rotators and an r^(1/4) surface brightness profile at small scales.

Gravitational-wave (GW) recoil of merging supermassive black holes (SMBHs) may influence the co-evolution of SMBHs and their host galaxies. We examine this possibility using SPH/N-body simulations of gaseous galaxy mergers in which the merged BH receives a recoil kick. With our suite of over 200 merger simulations, we identify systematic trends in the behavior of recoiling BHs. Our main results are as follows. (1) While BHs kicked at nearly the central escape speed (vesc) are essentially "lost" to the galaxy, in gas rich mergers, BHs kicked with up to about 0.7 vesc may be confined to the central few kpc of the galaxy. (2) The inflow of cold gas during a gas-rich major merger may cause a rapid increase in central escape speed; in such cases recoil trajectories will depend on the timing of the BH merger relative to the change in vesc. (3) Recoil events generally reduce the lifetimes of bright active galactic nuclei (AGN) but may actually extend AGN lifetimes at lower luminosities. (4) Recoiling AGN may be observable via kinematic offsets (v > 500 km s^-1) or spatial offsets (R > 1 kpc) for lifetimes of up to about 10 - 100 Myr. (5) Rapidly-recoiling BHs may be up to about 5 times less massive than their stationary counterparts. These mass deficits lower the normalization of the M - sigma relation and contribute to both intrinsic and overall scatter. (6) Finally, the displacement of AGN feedback by a recoil event causes higher central star formation rates in the merger remnant, thereby extending the starburst phase of the merger and creating a denser, more massive stellar cusp.

Gravitational-wave (GW) recoil of merging supermassive black holes (SMBHs) may influence the co-evolution of SMBHs and their host galaxies. We examine this possibility using SPH/N-body simulations of gaseous galaxy mergers in which the merged BH receives a recoil kick. This enables us to follow recoiling BHs in self-consistent, evolving merger remnants. In contrast to recent studies on similar topics, we conduct a large parameter study, generating a suite of over 200 simulations with more than 60 merger models and a range of recoil velocities (vk). Our main results are as follows. (1) BHs kicked at nearly the central escape speed (vesc) may oscillate on large orbits for up to a Hubble time, but in gas-rich mergers, BHs kicked with up to ~ 0.7 vesc may be confined to the central few kpc of the galaxy, owing to gas drag and steep central potentials. (2) vesc in gas-rich mergers may increase rapidly during final coalescence, in which case trajectories may depend on the timing of the BH merger relative to the formation of the potential well. (3) Recoil events generally reduce the lifetimes of bright active galactic nuclei (AGN), but may actually extend AGN lifetimes at lower luminosities. (4) Kinematically-offset AGN (v > 800 km s^-1) may be observable for up to ~ 10 Myr either immediately after the recoil or during pericentric passages through a gas-rich remnant. (5) Spatially-offset AGN (R > 1 kpc) generally have low luminosities and lifetimes of ~ 1 - 100 Myr. (6) Rapidly-recoiling BHs may be up to ~ 5 times less massive than their stationary counterparts. This lowers the normalization of the M-sigma relation and contributes to both intrinsic and overall scatter. (7) Finally, the displacement of AGN feedback after a recoil event enhances central star formation rates, thereby extending the starburst phase of the merger and creating a denser stellar cusp. [Abridged.]

We present an analysis of the large-scale galaxy distribution around two possible warm-hot intergalactic medium (WHIM) absorption systems reported along the Markarian 421 sightline. Using the Sloan Digital Sky Survey, we find a prominent galaxy filament at the redshift of the z=0.027 X-ray absorption line system. The filament exhibits a width of 3.2 Mpc and length of at least 20 Mpc, comparable to the size of WHIM filaments seen in cosmological simulations. No individual galaxies fall within 350 projected kpc so it is unlikely that the absorption is associated with gas in a galaxy halo or outflow. Another, lower-significance X-ray absorption system was reported in the same Chandra spectrum at z=0.011, but the large-scale structure in its vicinity is far weaker and may be a spurious alignment. By searching for similar galaxy structures in 140 random smoothed SDSS fields, we estimate a ~5-10% probability of the z=0.027 absorber-filament alignment occurring by chance. If these two systems are indeed physically associated, this would represent the first known coincidence between large-scale galaxy structure and a blind X-ray WHIM detection.

Theoretical models have had difficulty matching the observed number density of sub-millimeter galaxies (SMGs), causing some authors (e.g., Baugh et al. 2005) to suggest that SMGs provide evidence for a top-heavy initial mass function (IMF). To test this claim, we have, for the first time, combined high-resolution 3-D hydrodynamic simulations of isolated and merging massive, gas-rich galaxies, radiative transfer, and a semi-empirical merger rate model to predict the number density of SMGs. Our model can reproduce the observed SMG number density even when using a standard (Kroupa) IMF. The agreement is due to a combination of relatively long sub-mm duty cycles for mergers (a few times 10^8 years for our most massive models), which owe to our combination of high-resolution 3-D hydrodynamic simulations and dust radiative transfer; sufficient number densities of massive, gas-rich mergers; and the decrease in sub-mm counts observed by recent deep/wide surveys (e.g., Austermann et al. 2010) relative to previous surveys. Our results suggest that the observed SMG number counts do not provide evidence for a top-heavy IMF at high redshift.

We present a detailed analysis of the structure and resolved stellar populations of simulated merger remnants, and compare them to observations of compact quiescent galaxies at z ~ 2. We find that major merging is a viable mechanism to produce systems of ~ 10^11 Msun and ~ 1 kpc size, provided the gas fraction at the time of final coalescence is high (~ 40%), and provided that the progenitors are compact star-forming galaxies, as expected at high redshift. Their integrated spectral energy distributions and velocity dispersions are in good agreement with the observations, and their position in the (v_{maj}/sigma, ellipticity) diagram traces the upper envelope of the distribution of lower redshift early-type galaxies. The simulated merger remnants show time- and sightline-dependent M/L ratio gradients that result from a superposition of radially dependent stellar age, stellar metallicity, and extinction. The median ratio of effective radius in rest-frame V-band light to that in mass surface density is ~ 2 during the quiescent remnant phase. This is typically expressed by a negative color gradient (i.e., red core), which we expect to correlate with the integrated color of the system. Finally, the simulations differ from the observations in their surface brightness profile shape. The simulated remnants are typically best fit by high (n >> 4) Sersic indices, whereas observed quiescent galaxies at z ~ 2 tend to be less cuspy (median n ~ 2.3). Limiting early star formation in the progenitors may be required to prevent the simulated merger remnants from having extended wings.

Recent high precision proper motions from the Hubble Space Telescope (HST) suggest that the Large and Small Magellanic Clouds (LMC and SMC, respectively) are either on their first passage or on an eccentric long period (>6 Gyr) orbit about the Milky Way (MW). This differs markedly from the canonical picture in which the Clouds travel on a quasi-periodic orbit about the MW (period of ~2 Gyr). Without a short period orbit about the MW, the origin of the Magellanic Stream, a young (1-2 Gyr old) coherent stream of HI gas that trails the Clouds ~150 degrees across the sky, can no longer be attributed to stripping by MW tides and/or ram pressure stripping by MW halo gas. We propose an alternative formation mechanism in which material is removed by LMC tides acting on the SMC before the system is accreted by the MW. We demonstrate the feasibility and generality of this scenario using an N-body/SPH simulation with cosmologically motivated initial conditions constrained by the observations. Under these conditions we demonstrate that it is possible to explain the origin of the Magellanic Stream in a first infall scenario. This picture is generically applicable to any gas-rich dwarf galaxy pair infalling towards a massive host or interacting in isolation.

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.

The outer regions of galactic disks have received increased attention since ultraviolet observations with GALEX demonstrated that nearly 30% of galaxies have UV emission beyond their optical extents, indicating star formation activity. These galaxies have been termed extended UV (XUV) disks. Here, we address whether these observations contradict the gas surface density threshold for star formation inferred from Halpha radial profiles of galaxies. We run smoothed particle hydrodynamics simulations of isolated disk galaxies with fiducial star formation prescriptions and show that over-densities owing to the presence of spiral structure can induce star formation in extended gas disks. For direct comparison with observations, we use the 3-D radiative transfer code Sunrise to create simulated FUV and K_s band images. We find that galaxies classified as Type I XUV disks are a natural consequence of spiral patterns, but we are unable to reproduce Type II XUV disks. We also compare our results to studies of the Kennicutt-Schmidt relation in outer disks.

We present the first luminous, spatially resolved binary quasar that clearly inhabits an ongoing galaxy merger. SDSS J125455.09+084653.9 and SDSS J125454.87+084652.1 (SDSS J1254+0846 hereafter) are two luminous z=0.44 radio quiet quasars, with a radial velocity difference of just 215 km/s, separated on the sky by 21 kpc in a disturbed host galaxy merger showing obvious tidal tails. The pair was targeted as part of a complete sample of binary quasar candidates with small transverse separations drawn from SDSS DR6 photometry. We present follow-up optical imaging which shows broad, symmetrical tidal arm features spanning some 75 kpc at the quasars' redshift. Numerical modeling suggests that the system consists of two massive disk galaxies prograde to their mutual orbit, caught during the first passage of an active merger. This demonstrates rapid black hole growth during the early stages of a merger between galaxies with pre-existing bulges. Neither of the two luminous nuclei show significant instrinsic absorption by gas or dust in our optical or X-ray observations, illustrating that not all merging quasars will be in an obscured, ultraluminous phase. We find that the Eddington ratio for the fainter component B is rather normal, while for the A component L/LEdd is quite (>3sigma) high compared to quasars of similar luminosity and redshift, possibly evidence for strong merger-triggered accretion. More such mergers should be identifiable at higher redshifts using binary quasars as tracers.

Since the violent relaxation in hierarchical merging is incomplete, elliptical galaxies retain a wealth of information about their formation pathways in their present-day orbital structure. A variety of evidence indicates that gas-rich major mergers play an important role in the formation of elliptical galaxies. We simulate 1:1 disk mergers at seven different initial gas fractions ranging from 0 to 40%, using the TreeSPH code Gadget-2. We classify the stellar orbits in each remnant and construct radial profiles of the orbital content, intrinsic shape, and orientation. The dissipationless remnants are typically prolate-triaxial, dominated by box orbits within r_c ~ 1.5Reff, and by tube orbits in their outer parts. As the gas fraction increases, the box orbits within r_c are increasingly replaced by a population of short axis tubes (z-tubes) with near zero net rotation, and the remnants become progressively more oblate and round. The long axis tube (x-tube) orbits are highly streaming and relatively insensitive to the gas fraction, implying that their angular momentum is retained from the dynamically cold initial conditions. Outside r_c, the orbital structure is essentially unchanged by the gas. The 15-20% gas remnants often display disk-like kinematically distinct cores (KDCs). These remnants show an interesting resemblance, in both their velocity maps and intrinsic orbital structure, to the KDC galaxy NGC4365 (van den Bosch et al. 2008). At 30-40% gas, the remnants are rapidly rotating, with sharp embedded disks on ~ 1Reff scales. We predict a characteristic, physically intuitive orbital structure for 1:1 disk merger remnants, with a distinct transition between 1 and 3Reff that will be readily observable with combined data from the 2D kinematics surveys SAURON and SMEAGOL.

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

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

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

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.

Massive galaxies at high-z have smaller effective radii than those today, but similar central densities. Their size growth therefore relates primarily to the evolving abundance of low-density material. Various models have been proposed to explain this evolution, which have different implications for galaxy, star, and BH formation. We compile observations of spheroid properties as a function of redshift and use them to test proposed models. Evolution in progenitor gas-richness with redshift gives rise to initial formation of smaller spheroids at high-z. These systems can then evolve in apparent or physical size via several channels: (1) equal-density 'dry' mergers, (2) later major or minor 'dry' mergers with less-dense galaxies, (3) adiabatic expansion, (4) evolution in stellar populations & mass-to-light-ratio gradients, (5) age-dependent bias in stellar mass estimators, (6) observational fitting/selection effects. If any one of these is tuned to explain observed size evolution, they make distinct predictions for evolution in other galaxy properties. Only model (2) is consistent with observations as a dominant effect. It is the only model which allows for an increase in M_BH/M_bulge with redshift. Still, the amount of merging needed is larger than that observed or predicted. We therefore compare cosmologically motivated simulations, in which all these effects occur, & show they are consistent with all the observational constraints. Effect (2), which builds up an extended low-density envelope, does dominate the evolution, but effects 1,3,4, & 6 each contribute ~20% to the size evolution (a net factor ~2). This naturally also predicts evolution in M_BH-sigma similar to that observed.

The properties of elliptical galaxies are broadly consistent with simulated remnants of gas-rich mergers between spirals, motivating more detailed studies of the imprint of this formation mechanism on the remnant distribution function. Gas has a strong impact on the non-Gaussian shapes of the line-of-sight velocity distributions (LOSVDs) of the merger remnant, owing to the embedded disk that forms out of the gas that retains its angular momentum during the merger, and the strong central mass concentration from the gas that falls to the center. The deviations from Gaussianity are parametrized by the Gauss-Hermite moments h_3 and h_4, which are related to the skewness and kurtosis of the LOSVDs. We quantify the dependence of the (h_3,h_4)-v/sigma relations on the initial gas fraction of the progenitor disks in 1:1 mergers, using Gadget-2 simulations including star formation, radiative cooling, and feedback from supernovae and AGN. For gas fractions f_gas < ~15% the overall correlation between h_3 and v/sigma is weak, consisting of a flat negatively correlated component arising from edge-on viewing angles plus a steep positively correlated part from face-on projections. The spread in v/sigma values decreases toward high positive h_4, and there is a trend toward lower h_4 as the gas fraction increases from 0 to 15%. For f_gas > ~20% the (h_3,4)- v/sigma distributions look quite different - there is a tight negative h_3- v/sigma correlation, and a wide spread in v/sigma values at all h_4, in better agreement with observations. Re-mergers of the high-f_gas remnants (dry mergers) produce slowly rotating systems with nearly Gaussian LOSVDs. We explain all of these trends in terms of the underlying orbit structure of the remnants, as molded by their dissipative formation histories.

Transformation of disks into spheroids via mergers is a well-accepted element of galaxy formation models. However, recent simulations have shown that bulge formation is suppressed in increasingly gas-rich mergers. We investigate the global implications of these results in a cosmological framework, using independent approaches: empirical halo-occupation models (where galaxies are populated in halos according to observations) and semi-analytic models. In both, ignoring the effects of gas in mergers leads to the over-production of spheroids: low and intermediate-mass galaxies are predicted to be bulge-dominated (B/T~0.5 at <10^10 M_sun), with almost no bulgeless systems), even if they have avoided major mergers. Including the different physical behavior of gas in mergers immediately leads to a dramatic change: bulge formation is suppressed in low-mass galaxies, observed to be gas-rich (giving B/T~0.1 at <10^10 M_sun, with a number of bulgeless galaxies in good agreement with observations). Simulations and analytic models which neglect the similarity-breaking behavior of gas have difficulty reproducing the strong observed morphology-mass relation. However, the observed dependence of gas fractions on mass, combined with suppression of bulge formation in gas-rich mergers, naturally leads to the observed trends. Discrepancies between observations and models that ignore the role of gas increase with redshift; in models that treat gas properly, galaxies are predicted to be less bulge-dominated at high redshifts, in agreement with the observations. We discuss implications for the global bulge mass density and future observational tests.

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.

Dwarf spheroidal galaxies are the most dark matter dominated systems in the nearby Universe and their origin is one of the outstanding puzzles of how galaxies form. Dwarf spheroidals are poor in gas and stars, making them unusually faint, and those known as ultra-faint dwarfs have by far the lowest measured stellar content of any galaxy. Previous theories require that dwarf spheroidals orbit near giant galaxies like the Milky Way, but some dwarfs have been observed in the outskirts of the Local Group. Here we report simulations of encounters between dwarf disk galaxies and somewhat larger objects. We find that the encounters excite a process, which we term ``resonant stripping'', that can transform them into dwarf spheroidals. This effect is distinct from other mechanisms proposed to form dwarf spheroidals, including mergers, galaxy-galaxy harassment, or tidal and ram pressure stripping, because it is driven by gravitational resonances. It may account for the observed properties of dwarf spheroidals in the Local Group, including their morphologies and kinematics. Resonant stripping predicts that dwarf spheroidals should form through encounters, leaving detectable long stellar streams and tails.

An updated version of the dust radiation transfer code Sunrise, including models for star-forming regions and a self-consistent calculation of the spatially dependent dust and PAH emission, is presented. Given a hydrodynamic simulation of a galaxy, this model can calculate a realistic 2-dimensional ultraviolet--submillimeter spectral energy distribution of the galaxy, including emission lines from HII regions, from any viewpoint. To model the emission from star-forming regions, the MAPPINGSIII photoionization code is used. The high wavelength resolution (~ 1000 wavelengths) is made possible by the polychromatic Monte-Carlo algorithm employed by Sunrise. From the 2-D spectral energy distributions, images in any filter bands or integrated galaxy SEDs can be created. Using a suite of hydrodynamic simulations of disc galaxies, the output broad-band images and spectral energy distributions are compared with observed galaxies from the multiwavelength SINGS and SLUGS galaxy surveys. Overall, the output spectral energy distributions show a good match with observed galaxies in colours ranging from GALEX far-UV to SCUBA submillimeter wavelengths. The only possible exception is the 160 micron/850 micron colour, which the simulations underestimate by a factor "of order 5" compared to the SINGS sample. However, the simulations here agree with the SLUGS galaxies, which consistently have significantly larger amounts of cold dust than the SINGS galaxies. The Sunrise model can be used to generate simulated observations of arbitrary hydrodynamic galaxy simulations. In this way, predictions of galaxy formation theories can be directly tested against observations of galaxies.

We present a comparison between the observed color distribution, number and mass density of massive galaxies at 1.5 < z < 3 and a model by Hopkins et al. that relates the quasar and galaxy population on the basis of gas-rich mergers. In order to test the hypothesis that quiescent red galaxies are formed after a gas-rich merger involving quasar activity, we confront photometry of massive (M > 4x10^10 Msun) galaxies extracted from the FIRES, GOODS-South, and MUSYC surveys, together spanning an area of 496 arcmin^2, with synthetic photometry from hydrodynamical merger simulations. As in the Hopkins et al. (2006b) model, we use the observed quasar luminosity function to estimate the merger rate. We find that the synthetic U-V and V-J colors of galaxies that had a quasar phase in their past match the colors of observed galaxies that are best characterized by a quiescent stellar population. At z ~ 2.6, the observed number and mass density of quiescent red galaxies with M > 4x10^10 Msun is consistent with the model in which every quiescent massive galaxy underwent a quasar phase in the past. At z ~ 1.9, 2.8 times less quiescent galaxies are observed than predicted by the model as descendants of higher redshift quasars. The merger model also predicts a large number of galaxies undergoing merger-driven star formation. We find that the predicted number and mass density accounts for 30-50% of the observed massive star-forming galaxies. However, their colors do not match those of observed star-forming galaxies. In particular, the colors of dusty red galaxies are not reproduced by the simulations. Several possible origins of this discrepancy are discussed. The observational constraints on the validity of the model are currently limited by cosmic variance and uncertainties in stellar population synthesis and radiative transfer.

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 investigate how 'extra' central light in the surface brightness profiles of cusp ellipticals relates to the profiles of ellipticals with cores. Cusp elliptical envelopes are formed by violent relaxation in mergers acting on stars in progenitor disks, while their centers are structured by dissipational starbursts. Core ellipticals are formed by subsequent merging of (now gas-poor) cusp ellipticals, with the fossil starburst components combining to preserve a compact component in the remnant (although the 'transition' is smoothed). Comparing hydrodynamical simulations and observed profiles, we show how to observationally isolate the relic starburst components in core ellipticals. We demonstrate that these survive re-mergers and reliably trace the dissipation in the initial gas-rich merger(s). The typical degree of dissipation is a strong function of stellar mass, tracing observed disk gas fractions. We find a correlation between dissipation and effective radius: systems with more dissipation are more compact. The survival of this component and scattering of stars into the envelope naturally explain high-Sersic index profiles characteristic of massive core ellipticals. This is also closely related to the kinematics and isophotal shapes: only systems with matched starburst components from their profile fits also reproduce the observed kinematics of boxy/core ellipticals. We show that it is critical to adopt physically motivated profiles when attempting to quantify how much mass has been 'scoured' or scattered out of the inner regions by binary black holes. Estimates of scoured mass ignoring multi-component structure can be strongly biased, potentially explaining observed systems with large inferred core masses in apparent conflict with core-scouring models.

We present a detailed comparison of the CO(3-2) emitting molecular gas between a local sample of luminous infrared galaxies (U/LIRGs) and a high redshift sample that comprises submm selected galaxies (SMGs), quasars, and Lyman Break Galaxies (LBGs). The U/LIRG sample consists of our recent CO(3-2) survey using the Submillimeter Array while the CO(3-2) data for the high redshift population are obtained from the literature. We find that the L(CO(3-2)) and L(FIR) relation is correlated over five orders of magnitude, which suggests that the molecular gas traced in CO(3-2) emission is a robust tracer of dusty star formation activity. The near unity slope of 0.93 +/- 0.03 obtained from a fit to this relation suggests that the star formation efficiency is constant to within a factor of two across different types of galaxies residing in vastly different epochs. The CO(3-2) size measurements suggest that the molecular gas disks in local U/LIRGs (0.3 - 3.1 kpc) are much more compact than the SMGs (3 - 16 kpc), and that the size scales of SMGs are comparable to the nuclear separation (5 - 40 kpc) of the widely separated nuclei of U/LIRGs in our sample. We argue from these results that the SMGs studied here are predominantly intermediate stage mergers, and that the wider line-widths arise from the violent merger of two massive gas-rich galaxies taking place deep in a massive halo potential.

We present a detailed analysis of our ability to determine stellar masses, ages, reddening and extinction values, and star formation rates of high-redshift galaxies by modeling broad-band SEDs with stellar population synthesis. In order to do so, we computed synthetic optical-to-NIR SEDs for model galaxies taken from hydrodynamical merger simulations placed at redshifts 1.5 < z < 3. Viewed under different angles and during different evolutionary phases, the simulations represent a wide variety of galaxy types (disks, mergers, spheroids). We show that simulated galaxies span a wide range in SEDs and color, comparable to these of observed galaxies. In all star-forming phases, dust attenuation has a large effect on colors, SEDs, and fluxes. The broad-band SEDs were then fed to a standard SED modeling procedure and resulting stellar population parameters were compared to their true values. Disk galaxies generally show a decent median correspondence between the true and estimated mass and age, but suffer from large uncertainties. During the merger itself, we find larger offsets (e.g., log M_recovered - log M_true = -0.13^{+0.10}_{-0.14}). E(B-V) values are generally recovered well, but the estimated total visual absorption Av is consistently too low, increasingly so for larger optical depths. Since the largest optical depths occur during the phases of most intense star formation, it is for the highest SFRs that we find the largest underestimates. The masses, ages, E(B-V), Av, and SFR of merger remnants (spheroids) are very well reproduced. We discuss possible biases in SED modeling results caused by mismatch between the true and template star formation history, dust distribution, metallicity variations and AGN contribution.

We study the origin and properties of 'extra' or 'excess' central light in the surface brightness profiles of cusp or power-law ellipticals. Dissipational mergers give rise to two-component profiles: an outer profile established by violent relaxation acting on stars present in the progenitors prior to the final merger, and an inner stellar population comprising the extra light, formed in a compact starburst. Combining a large set of hydrodynamical simulations with data that span a broad range of profiles and masses, we show that this picture is borne out -- cusp ellipticals are indeed 'extra light' ellipticals -- and examine how the properties of this component scale with global galaxy properties. We show how to robustly separate the 'extra' light, and demonstrate that observed cusps are reliable tracers of the degree of dissipation in the spheroid-forming merger. We show that the typical degree of dissipation is a strong function of stellar mass, tracing observed disk gas fractions at each mass. We demonstrate a correlation between extra light content and effective radius at fixed mass: systems with more dissipation are more compact. The outer shape of the light profile does not depend on mass, with a mean outer Sersic index ~2.5. We explore how this relates to shapes, kinematics, and stellar population gradients. Simulations with the gas content needed to match observed profiles also reproduce observed age, metallicity, and color gradients, and we show how these can be used as tracers of the degree of dissipation in spheroid formation.

We develop a model for the origins and redshift evolution of spheroid scaling relations. We consider spheroid sizes, velocity dispersions, masses, profile shapes (Sersic indices), and black hole (BH) masses, and their related scalings. Our approach combines advantages of observational constraints in halo occupation models and hydrodynamic merger simulations. This allows us to separate the relative roles of dissipation, dry mergers, formation time, and progenitor evolution, and identify their effects on scalings at each redshift. Dissipation is the most important factor determining spheroid sizes and fundamental plane (FP) scalings, and can account for the FP tilt and differences between disk and spheroid scalings. Because disks at high-z have higher gas fractions, mergers are more gas-rich, yielding more compact spheroids. This predicts mass-dependent evolution in spheroid sizes, in agreement with observations. This relates to subtle evolution in the FP, important to studies that assume a fixed intrinsic FP. This also predicts mild evolution in BH-host correlations, towards larger BHs at higher z. Dry mergers are significant, but only for massive systems which form early: they form compact, but undergo dry mergers (consistent with observations) such that their sizes at later times are similar to spheroids of similar mass formed more recently. We model descendants of observed compact high-z spheroids: most will become cores of BCGs, with sizes, velocity dispersions, and BH masses consistent with observations, but we identify a fraction that might survive to z=0 intact.

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.

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

The HST proper motion (PM) measurements of the Clouds have severe implications for their interaction history with the Milky Way (MW) and with each other. The Clouds are likely on their first passage about the MW and the SMC's orbit about the LMC is better described as quasi-periodic rather than circular. Binary L/SMC orbits that satisfy observational constraints on their mutual interaction history (e.g. the formation of the Magellanic Bridge during a collision between the Clouds ~300 Myr ago) can be located within 1 sigma of the mean PMs. However, these binary orbits are not co-located with the Magellanic Stream (MS) when projected on the plane of the sky and the line-of-sight velocity gradient along the LMC's orbit is significantly steeper than that along the MS. These combined results ultimately rule out a purely tidal origin for the MS: tides are ineffective without multiple pericentric passages and can neither decrease the velocity gradient nor explain the offset stream in a polar orbit configuration. Alternatively, ram pressure stripping of an extended gaseous disk may naturally explain the deviation. The offset also suggests that observations of the little-explored region between RA 21h and 23h are crucial for characterizing the full extent of the MS.

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 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).

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 develop observational tests of the idea that dissipation in gas-rich mergers produces the fundamental plane (FP) and related correlations obeyed by ellipticals. The FP 'tilt' implies lower-mass ellipticals have a higher ratio of stellar to dark matter within their stellar effective radii. Models argue that mergers between more gas-rich (typically lower-mass) disks yield larger mass fractions formed in compact starbursts, giving a smaller stellar R_e and higher M_stellar/M_tot within that R_e. Such starbursts leave a characteristic imprint in the surface brightness profile: a central excess above an outer profile established by the dissipationless violent relaxation of disk stars. In previous work, we developed empirical methods to decompose the observed profiles of ellipticals and robustly estimate the amount of dissipation in the original spheroid-forming merger(s). Applying this to a large sample of observed ellipticals, we test whether or not their location on the FP and its tilt are driven by dissipation. At fixed mass, ellipticals formed in more dissipational events are smaller and have higher M_stellar/M_tot. At fixed degree of dissipation, there is no tilt in the FP. We show that the dynamical mass estimator R_e*sigma^2/G is a good estimator of the true mass: the observed FP tilt cannot primarily owe to other forms of non-homology. Removing the effects of dissipation, observed ellipticals obey the same FP correlations as disks: unusual progenitors are not required to make typical ellipticals. Dissipation appears to be both necessary and sufficient to explain the FP tilt.

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.

We present new data obtained with the Submillimeter Array for a sample of fourteen nearby luminous and ultraluminous infrared galaxies. The galaxies were selected to have luminosity distances D < 200 Mpc and far-infrared luminosities log(L_FIR) > 11.4. The galaxies were observed with spatial resolutions of order 1 kpc in the CO J=3-2, CO J=2-1, 13CO J=2-1, and HCO+ J=4-3 lines as well as the continuum at 880 microns and 1.3 mm. We have combined our CO and continuum data to measure an average gas-to-dust mass ratio of 120 +/- 28 (rms deviation 109) in the central regions of these galaxies, very similar to the value of 150 determined for the Milky Way. This similarity is interesting given the more intense heating from the starburst and possibly accretion activity in the luminous infrared galaxies compared to the Milky Way. We find that the peak H_2 surface density correlates with the far-infrared luminosity, which suggests that galaxies with higher gas surface densities inside the central kiloparsec have a higher star formation rate. The lack of a significant correlation between total H_2 mass and far-infrared luminosity in our sample suggests that the increased star formation rate is due to the increased availability of molecular gas as fuel for star formation in the central regions. In contrast to previous analyses by other authors, we do not find a significant correlation between central gas surface density and the star formation efficiency, as trace by the ratio of far-infrared luminosity to nuclear gas mass. Our data show that it is the star formation rate, not the star formation efficiency, that increases with increasing central gas surface density in these galaxies.

We investigate the relationship between the star formation rate (SFR) and dense molecular gas mass in the nuclei of galaxies. To do this, we utilize the observed 850 micron luminosity as a proxy for the infrared luminosity and SFR, and correlate this with the observed CO (J=3-2) luminosity. We find tentative evidence that the LIR-CO (J=3-2) index is similar to the Kennicutt-Schmidt (KS) index (N ~ 1.5) in the central ~1.7 kpc of galaxies, and flattens to a roughly linear index when including emission from the entire galaxy. This result may imply that the volumetric Schmidt relation is the underlying driver behind the observed SFR-dense gas correlations, and provides tentative confirmation for recent numerical models. While the data exclude the possibility of a constant LIR-CO (J=3-2) index for both galaxy nuclei and global measurements at the ~80% confidence level, the considerable error bars cannot preclude alternative interpretations.

We study the origin and properties of 'extra' or 'excess' central light in the surface brightness profiles of gas-rich merger remnants. Combining a large set of hydrodynamical simulations with data on observed mergers (spanning a broad range of profiles at various masses and degrees of relaxation), we show how to robustly separate the physically meaningful extra light -- stellar populations formed in a compact central starburst during a gas-rich merger -- from the outer profile established by violent relaxation acting on stars already present in the progenitors prior to the final merger. This separation is sensitive to the profile treatment, and we demonstrate that certain fitting procedures can yield physically misleading results. We show that our method reliably recovers the younger starburst population, and examine how the properties of this component scale with mass, gas content, and other aspects of the progenitors. We consider the time evolution of profiles in different bands, and estimate biases introduced by observational studies at different times and wavelengths. We show that extra light is ubiquitous in observed and simulated gas-rich merger remnants, with sufficient mass (~3-30% of the stellar mass) to explain the discrepancy in the maximum phase-space densities of ellipticals and their progenitor spirals. The nature of this central component provides powerful new constraints on the formation histories of observed systems.

We derive a physical model for the observed relations between star formation rate (SFR) and molecular line (CO and HCN) emission in galaxies, and show how these observed relations are reflective of the underlying star formation law. We do this by combining 3D non-LTE radiative transfer calculations with hydrodynamic simulations of isolated disk galaxies and galaxy mergers. We demonstrate that the observed SFR-molecular line relations are driven by the relationship between molecular line emission and gas density, and anchored by the index of the underlying Schmidt law controlling the SFR in the galaxy. Lines with low critical densities (e.g. CO J=1-0) are typically thermalized and trace the gas density faithfully. In these cases, the SFR will be related to line luminosity with an index similar to the Schmidt law index. Lines with high critical densities greater than the mean density of most of the emitting clouds in a galaxy (e.g. CO J=3-2, HCN J=1-0) will have only a small amount of thermalized gas, and consequently a superlinear relationship between molecular line luminosity and mean gas density. This results in a SFR-line luminosity index less than the Schmidt index for high critical density tracers. One observational consequence of this is a significant redistribution of light from the small pockets of dense, thermalized gas to diffuse gas along the line of sight, and prodigious emission from subthermally excited gas. At the highest star formation rates, the SFR-Lmol slope tends to the Schmidt index, regardless of the molecular transition. The fundamental relation is the Kennicutt-Schmidt law, rather than the relation between SFR and molecular line luminosity. We use these results to make imminently testable predictions for the SFR-molecular line relations of unobserved transitions.

We use a N--body/hydrodynamic simulation to forecast the future encounter between the Milky Way and the Andromeda galaxies, given current observational constraints on their relative distance, relative velocity, and masses. Allowing for a comparable amount of diffuse mass to fill the volume of the Local Group, we find that the two galaxies are likely to collide in a few billion years - within the Sun's lifetime. During the the interaction, there is a chance that the Sun will be pulled away from its present orbital radius and reside in an extended tidal tail. The likelihood for this outcome increases as the merger progresses, and there is a remote possibility that our Sun will be more tightly bound to Andromeda than to the Milky Way before the final merger. Eventually, after the merger has completed, the Sun is most likely to be scattered to the outer halo and reside at much larger radii (>30 kpc). The density profiles of the stars, gas and dark matter in the merger product resemble those of elliptical galaxies. Our Local Group model therefore provides a prototype progenitor of late--forming elliptical galaxies.

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

We assess the impact of starburst and AGN feedback-driven winds on the CO emission from galaxy mergers, and, in particular, search for signatures of these winds in the simulated CO morphologies and emission line profiles. We do so by combining a 3D non-LTE molecular line radiative transfer code with smoothed particle hydrodynamics (SPH) simulations of galaxy mergers that include prescriptions for star formation, black hole growth, a multiphase interstellar medium (ISM), and the winds associated with star formation and black hole growth. Our main results are: (1) Galactic winds can drive outflows of masses ~10^8-10^9 Msun which may be imaged via CO emission line mapping. (2) AGN feedback-driven winds are able to drive imageable CO outflows for longer periods of time than starburst-driven winds owing to the greater amount of energy imparted to the ISM by AGN feedback compared to star formation. (3) Galactic winds can control the spatial extent of the CO emission in post-merger galaxies, and may serve as a physical motivation for the sub-kiloparsec scale CO emission radii observed in local advanced mergers. (4) Secondary emission peaks at velocities greater than the circular velocity are seen in the CO emission lines in all models. In models with winds, these high velocity peaks are seen to preferentially correspond to outflowing gas entrained in winds, which is not the case in the model without winds. The high velocity peaks seen in models without winds are typically confined to velocity offsets (from the systemic) < 1.7 times the circular velocity, whereas the models with AGN feedback-driven winds can drive high velocity peaks to ~2.5 times the circular velocity.

The detection of large quantities of dust in z ~ 6 quasars by infrared and radio surveys presents puzzles for the formation and evolution of dust in these early systems. Previously (Li et al. 2007), we showed that luminous quasars at z > 6 can form through hierarchical mergers of gas-rich galaxies. Here, we calculate the dust properties of simulated quasars and their progenitors using a three-dimensional Monte Carlo radiative transfer code, ART^2 -- All-wavelength Radiative Transfer with Adaptive Refinement Tree. ART^2 incorporates a radiative equilibrium algorithm for dust emission, an adaptive grid for inhomogeneous density, a multiphase model for the ISM, and a supernova-origin dust model. We reproduce the SED and dust properties of SDSS J1148+5251, and find that the infrared emission are closely associated with the formation and evolution of the quasar host. The system evolves from a cold to a warm ULIRG owing to heating and feedback from stars and AGN. Furthermore, the AGN has significant implications for the interpretation of observation of the hosts. Our results suggest that vigorous star formation in merging progenitors is necessary to reproduce the observed dust properties of z~6 quasars, supporting a merger-driven origin for luminous quasars at high redshifts and the starburst-to-quasar evolutionary hypothesis. (Abridged)

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

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

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

We present the 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 examine the origin and evolution of correlations between properties of supermassive black holes (BHs) and their host galaxies using simulations of major galaxy mergers, including the effects of gas dissipation, cooling, star formation, and BH accretion and feedback. We demonstrate that the simulations predict the existence of a BH 'fundamental plane' (BHFP), of the form M_BH sigma^(3.0+-0.3)*R_e^(0.43+-0.19) or M_BH M_bulge^(0.54+-0.17)*sigma^(2.2+-0.5), similar to relations found observationally. The simulations indicate that the BHFP can be understood roughly as a tilted intrinsic correlation between BH mass and spheroid binding energy, or the condition for feedback coupling to power a pressure-driven outflow. While changes in halo circular velocity, merger orbital parameters, progenitor disk redshifts and gas fractions, ISM gas pressurization, and other parameters can drive changes in e.g. sigma at fixed M_bulge, and therefore changes in the M_BH-sigma or M_BH-M_bulge relations, the BHFP is robust. Given the empirical trend of decreasing R_e for a given M_bulge at high redshift, the BHFP predicts that BHs will be more massive at fixed M_bulge, in good agreement with recent observations. This evolution in the structural properties of merger remnants, to smaller R_e and larger sigma (and therefore larger M_BH, conserving the BHFP) at a given M_bulge, is driven by the fact that bulge progenitors have characteristically larger gas fractions at high redshifts. Adopting the observed evolution of disk gas fractions with redshift, our simulations predict the observed trends in both R_e(M_bulge) and M_BH(M_bulge).

We study observed correlations between supermassive black hole (BHs) and the properties of their host galaxies, and show that the observations define a BH 'fundamental plane' (BHFP), of the form M_BH sigma^(3.0+-0.3)*R_e^(0.43+-0.19), or M_BH M_bulge^(0.54+-0.17)*sigma^(2.2+-0.5), analogous to the FP of elliptical galaxies. The BHFP is preferred over a simple relation between M_BH and any of sigma, M_bulge, M_dyn, or R_e alone at >99.9% significance. The existence of this BHFP has important implications for the formation of supermassive BHs and the masses of the very largest black holes, and immediately resolves several apparent conflicts between the BH masses expected and measured for outliers in both the M_BH-sigma and M_BH-M_bulge relations.

We investigate the nature of CO emission from z~6 quasars by combining non-LTE radiative transfer calculations with merger-driven models of z~6 quasar formation that arise naturally in LCDM cosmological simulations. We consider four model quasars formed in 10^12-10^13 M_sun halos from different merging histories. Our main results follow. Owing to massive starbursts and funneling of dense gas into the nuclear regions of merging galaxies, the CO is highly excited and the flux density peaks between J=5-8. The CO morphology of z~6 quasars often exhibits multiple emission peaks which arise from H2 concentrations which have not yet fully coalesced. Quasars at z~6 display a large range of sightline dependent line widths such that the lines are narrowest when the rotating H2 gas associated with the quasar is viewed face-on (when L_B is largest), and broadest when the gas is seen edge-on (when L_B is lowest). Thus for all models selection effects exist such that quasars selected for optical luminosity are preferentially face-on which may result in detected CO line widths narrower than the median. The sightline averaged line width is reflective of the circular velocity (V_c) of the host halo, and ranges from sigma~300-650 km/s. For optically selected QSOs, 10-25% (halo-mass dependant) of sightlines have narrow line widths compatible with the sole CO detection at z>6, J1148+5251. When accounting for both the temporal evolution of CO line widths, as well as the redshift evolution of halo V_c, these models self-consistently account for the CO line widths of both z~2 sub-mm galaxies and QSO's. Finally, the dynamical mass derived from the sightline averaged line widths provides a good estimate of the total mass, and allows for a stellar bulge and SMBH consistent with the local M_BH-M_bulge relation. [abridged]

We calculate the observable properties of the most massive high-redshift galaxies in the hierarchical formation scenario where stellar spheroid and supermassive black hole growth are fueled by gas-rich mergers. Combining high-resolution hydrodynamical simulations of the hierarchical formation of a z~6 quasar, stellar population synthesis models, template AGN spectra, prescriptions for interstellar and intergalactic absorption, and the response of modern telescopes, the photometric evolution of galaxies destined to host z~6 quasars are modeled at redshifts z~4-14. These massive galaxies, with enormous stellar masses of M_star ~10^11.5-10^12 M_sun. and star formation rates of SFR~10^3-10^4 M_sun yr^-1 at z>~7, satisfy a variety of photometric selection criteria based on Lyman-break techniques including V-band dropouts at z>~5, i-band dropouts at z>~6, and z-band dropouts at z>~7. The observability of the most massive high-redshift galaxies is assessed and compared with a wide range of existing and future photometric surveys including SDSS, GOODS/HUDF, NOAO WDFS, UKIDSS, the IRAC Shallow Survey, Pan-STARRS, LSST, and SNAP. Massive stellar spheroids descended from z~6 quasars will likely be detected at z~4 by existing surveys, but owing to their low number densities the discovery of quasar progenitor galaxies at z>7 will likely require future surveys of large portions of the sky (>~0.5%) at wavelengths lambda>1 micron. The detection of rare, star-bursting, massive galaxies at redshifts z>~6 would provide support for the hierarchical formation of the earliest quasars and characterize the primitive star-formation histories of the most luminous elliptical galaxies.

The discovery of luminous quasars at redshift z ~ 6 indicates the presence of supermassive black holes (SMBHs) of mass ~10^9 Msun when the Universe was less than one billion years old. This finding presents several challenges for theoretical models. Here, we present the first multi-scale simulations that, together with a self-regulated model for the SMBH growth, produce a luminous quasar at z ~ 6.5 in the LCDM paradigm. We follow the hierarchical assembly history of the most massive halo in a ~ 3 Gpc^3 volume, and find that this halo of ~ 8x 10^{12} Msun forming at z ~ 6.5 after several major mergers is able to reproduce a number of observed properties of SDSS J1148+5251, the most distant quasar detected at z =6.42 (Fan et al. 2003). Moreover, the SMBHs grow through gas accretion below the Eddington limit in a self-regulated manner owing to feedback. We find that the progenitors experience significant star formation (up to 10^4 Msun/yr) preceding the major quasar phase such that the stellar mass of the quasar host reaches 10^{12} Msun at z ~ 6.5, consistent with observations of significant metal enrichment in SDSS J1148+5251. Our results provide a viable formation mechanism for z ~ 6 quasars in the standard LCDM cosmology, and demonstrate a common, merger-driven origin for the rarest quasars and the fundamental SMBH-host correlation in a hierarchical Universe.(Abridged)

We compare observed clustering of quasars and galaxies as a function of redshift, mass, luminosity, & color/morphology, to constrain models of quasar fueling and spheroid-BH co-evolution. High redshift quasars are shown to be drawn from progenitors of local early-type galaxies, with the characteristic quasar luminosity L* reflecting a characteristic mass of 'active' BH/host populations at each redshift. Evolving observed high-z quasar clustering to z=0 predicts a trend of clustering in 'quasar remnants' as a function of stellar mass identical to that observed for early-types. However, quasar clustering does not simply reflect observed early (or late)-type populations; at each redshift, quasars cluster as an 'intermediate' population. Comparing with the age of elliptical stellar populations reveals that this 'intermediate' population represents those ellipticals undergoing or terminating their final significant star formation at each epoch. Assuming that quasar triggering is associated with the formation/termination epoch of ellipticals predicts quasar clustering at all observed redshifts without any model dependence or assumptions about quasar light curves, lifetimes, or accretion rates. This is not true for disks or quasar halos: i.e. quasars do not generically trace star formation, disks, or halo assembly. Quasar clustering at all z is consistent with a constant halo mass ~4x10^{12} M_sun, similar to local 'group scales.' The observations support a scenario in which major mergers trigger quasar activity and dominate bright, high-z quasar populations. We show that measurements of quasar clustering versus luminosity at z~1 can be used to constrain different lower-luminosity AGN fueling mechanisms, and that high-z clustering is sensitive to whether or not systems 'shut down' growth at z>3.

(Abridged) We develop and test a model for the cosmological role of mergers in the formation and quenching of red, early-type galaxies. Making the ansatz that star formation is quenched after a gas-rich, spheroid-forming major merger, we demonstrate that this naturally predicts the turnover in the efficiency of star formation at ~L_star, as well as the observed mass functions/density of red galaxies as a function of redshift, the formation times of spheroids as a function of mass, and the fraction of quenched galaxies as a function of galaxy and halo mass, environment, and redshift. Comparing to a variety of semi-analytic models in which quenching is primarily driven by halo mass considerations or secular/disk instabilities, we demonstrate that our model and different broad classes of models make unique and robust qualitative predictions for a number of observables, including the red fraction as a function of galaxy and halo mass, the density of passive galaxies and evolution of the color-morphology-density relations at high z, and the fraction of disky/boxy spheroids as a function of mass. In each case, the observations favor a model in which galaxies quench after a major merger builds a massive spheroid, and disfavor quenching via secular or pure halo processes. We discuss a variety of physical possibilities for this quenching, and propose a mixed scenario in which traditional quenching in hot, massive halos is supplemented by the feedback associated with star formation and quasar activity in a major merger, which temporarily suppress cooling and establish the conditions of a dynamically hot halo in the central regions of the host, even in low mass halos.

(Abridged) We develop a model for the cosmological role of mergers in the evolution of starbursts, quasars, and spheroidal galaxies. Combining halo mass functions (MFs) with empirical halo occupation models, we calculate where major galaxy-galaxy mergers occur and what kinds of galaxies merge, at all redshifts. We compare with observed merger MFs, clustering, fractions, and small-scale environments, and show that this yields robust estimates in good agreement with observations. Making the simple ansatz that major, gas-rich mergers cause quasar activity, we demonstrate that this naturally reproduces the observed rise and fall of the quasar luminosity density from z=0-6, as well as quasar LFs, fractions, host galaxy colors, and clustering as a function of redshift and luminosity. The observed excess of quasar clustering on small scales is a natural prediction of the model, as mergers preferentially occur in regions with excess small-scale galaxy overdensities. We show that quasar environments at all observed redshifts correspond closely to the empirically determined small group scale, where mergers of gas-rich galaxies are most efficient. We contrast with a secular model in which quasar activity is driven by bars/disk instabilities, and show that while these modes probably dominate at Seyfert luminosities, the constraints from clustering (large and small-scale), pseudobulge populations, disk MFs, luminosity density evolution, and host galaxy colors argue that they must be a small contributor to the z>1 quasar luminosity density.

Recent proper motion measurements of the Large and Small Magellanic Clouds (LMC and SMC, respectively) by Kallivayalil et al (2006a,b) suggest that the 3D velocities of the Clouds are substantially higher (~100 km/s) than previously estimated and now approach the escape velocity of the Milky Way (MW). Previous studies have also assumed that the Milky Way can be adequately modeled as an isothermal sphere to large distances. Here we re-examine the orbital history of the Clouds using the new velocities and a LCDM-motivated MW model with virial mass Mvir = 1e12 Msun (e.g. Klypin et al 2002). We conclude that the Clouds are either currently on their first passage about the MW or, if the MW can be accurately modeled by an isothermal sphere to distances >200 kpc (i.e., Mvir > 2e12 Msun), that their orbital period and apogalacticon distance must be a factor of two larger than previously estimated, increasing to 3 Gyr and 200 kpc, respectively. A first passage scenario is consistent with the fact that the LMC and SMC appear to be outliers when compared to other satellite galaxies of the MW: they are irregular in appearance and are moving faster. We discuss the implications of this orbital analysis for our understanding of the star formation history, the nature of the warp in the MW disk and the origin of the Magellanic Stream (MS), a band of HI gas trailing the LMC and SMC that extends roughly 100 degrees across the sky. Specifically, as a consequence of the new orbital history of the Clouds, the origin of the MS may not be explainable by current tidal and ram pressure stripping models.

We calculate infrared spectral energy distributions (SEDs) from simulations of major galaxy mergers and study the effect of AGN and starburst driven feedback on the evolution of the SED as a function of time. We use a self-consistent three-dimensional radiative equilibrium code to calculate the emergent SEDs and to make images. To facilitate a simple description of our findings, we describe our results in reference to an approximate analytic solution for the far-IR SED. We focus mainly on the luminous infrared galaxy (LIRG) and ultraluminous infrared galaxy (ULIRG) phases of evolution. We contrast the SEDs of simulations performed with AGN feedback to simulations performed with starburst driven wind feedback. We find that the feedback processes critically determine the evolution of the SED. Changing the source of illumination (whether stellar or AGN) has virtually no impact on the reprocessed far-infrared SED. We find that AGN feedback is particularly effective at dispersing gas and rapidly injecting energy into the ISM. The observational signature of such powerful feedback is a warm SED. In general, simulations performed with starburst driven winds have colder spectra and reprocess more of their emission into the infrared, resulting in higher infrared to bolometric luminosities compared to (otherwise equivalent) simulations performed with AGN feedback. We depict our results in IRAS bands, as well as in Spitzer's MIPS bands, and in Herschel's PACS bands.

[abridged]We calculate multi-wavelength spectral energy distributions (SEDs) (spanning optical to millimeter wavelengths) from simulations of major galaxy mergers with black hole feedback which produce submillimeter bright galaxies (SMGs), using a self-consistent three-dimensional radiative transfer code. We reproduce correlations for local AGN observed in Spitzer Space Telescope's IRAC bands, and make definitive predictions for infrared X-ray correlations that should be testable by combining observations by Spitzer and the upcoming Herschel mission with X-ray surveys. Our dynamical approach allows us to directly correlate observed clustering in the data as seen in IRAC color-color plots with the relative amount of time the system spends in a region of color-color space. We also find that this clustering is positively correlated with the stars dominating in their contribution to the total bolometric luminosity. We present photo albums spanning the lifetime of SMGs, from their infancy in the pre-merger phase to the final stage as an elliptical galaxy, as seen in the observed $3.6 \micron$ and $450 \micron$ band to visually illustrate some of the morphological differences between mergers of differing orbital inclination and progenitor redshift. We find that SMGs are a broader class of systems than starbursts or quasars. We introduce a simple, heuristic classification scheme on the basis of the $L_{\rm IR}/L_{\rm x}$ ratios of these galaxies, which may be interpreted qualitatively as an evolutionary scheme, as these galaxies evolve in $L_{\rm IR}/L_{\rm x}$ while transiting from the pre-merger stage, through the quasar phase, to a merger remnant.

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.

(Abridged) We study the link between quasars and the red galaxy population using a model for the self-regulated growth of supermassive black holes in mergers involving gas-rich galaxies. Using a model for quasar lifetimes and evolution motivated by hydrodynamical simulations of galaxy mergers, we de-convolve the observed quasar luminosity function at various redshifts to determine the rate of formation of black holes of a given final mass. Identifying quasar activity with the formation of spheroids in the framework of the merger hypothesis, this enables us to deduce the corresponding rate of formation of spheroids with given properties as a function of redshift. This allows us to predict, for the red galaxy population, the distribution of galaxy velocity dispersions, the mass function, mass density, star formation rates, the luminosity function in many observed wavebands (NUV, U, B, V, R, I, J, H, K), the total red galaxy number density and luminosity density, the distribution of colors as a function of magnitude and velocity dispersion for several different wavebands, the distribution of mass to light ratios vs. mass, the luminosity-size relations, and the typical ages and distribution of ages (formation redshifts) as a function of both mass and luminosity. For each of these quantities, we predict the evolution from redshift z=0-6. Each of our predictions agrees well with existing observations, without the addition of tunable parameters; the essential observational inputs come from the observed quasar luminosity function. These predictions are skewed by several orders of magnitude if we adopt simpler, traditional models of quasar lifetimes in which quasars turn on/off or follow simple exponential light curves, instead of the more complicated evolution implied by our simulations.

We describe a model-independent integral constraint which defines an upper limit to the allowed degree of evolution in the ratio of black hole (BH) mass to host galaxy luminosity or mass, as a function of redshift. Essentially, if the BH/host ratio is excessive at redshift z, then it would imply that the total mass density in BHs above some M_BH(min) is larger at that redshift than at z=0, which is impossible. This argument requires no knowledge of host or BH properties, only a lower limit to the observed luminosity density in the brightest galaxies at some z. We calculate this upper limit from a variety of luminosity and mass functions in different bands from redshifts z=0-2. We show that it is consistent with passive evolution of spheroid populations (with a fixed M_BH/M_host relation) in all cases, and provides tighter constraints than have generally been obtained previously, ruling out at >6 sigma observational and theoretical estimates suggesting that M_BH/M_host was significantly larger at high redshifts than locally, although relatively weak (factor ~2 by z=2) evolution is still allowed. We discuss a variety of possible 'loopholes' or changes in the BH/host populations and correlations, and show that they typically lower the upper limits and strengthen our conclusions.

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

We use numerical simulations to study the kinematic structure of remnants formed from mergers of equal-mass disk galaxies. In particular, we show that remnants of dissipational mergers, which include the radiative cooling of gas, star formation, feedback from supernovae, and the growth of supermassive black holes, are smaller, rounder, have, on average, a larger central velocity dispersion, and show significant rotation compared to remnants of dissipationless mergers. The increased rotation speed of dissipational remnants owes its origin to star formation that occurs in the central regions during the galaxy merger. We have further quantified the anisotropy, three-dimensional shape, minor axis rotation, and isophotal shape of each merger remnant, finding that dissipational remnants are more isotropic, closer to oblate, have the majority of their rotation along their major axis, and are more disky than dissipationless remnants. Individual remnants display a wide variety of kinematic properties. A large fraction of the dissipational remnants are oblate isotropic rotators. Many dissipational, and all of the dissipationless, are slowly rotating and anisotropic. The remnants of gas-rich major mergers can well-reproduce the observed distribution of projected ellipticities, rotation parameter (V/\sigma)*, kinematic misalignments, Psi, and isophotal shapes. The dissipationless remnants are a poor match to this data. Our results support the merger hypothesis for the origin of low-luminosity elliptical galaxies provided that the progenitor disks are sufficiently gas-rich, however our remnants are a poor match to the bright ellipticals that are slowly rotating and uniformly boxy.

We examine X-ray emission produced from hot gas during collisions and mergers of disk galaxies. To study this process, we employ simulations that incorporate cosmologically motivated disk-galaxy models and include the effects of radiative cooling, star formation, supernova feedback, and accreting supermassive black holes. We find that during a merger, the colliding gas in the disks is shock-heated to X-ray-emitting temperatures. The X-ray luminosity is spatially extended, rises during the initial stages of the merger, and peaks when the galactic centers coalesce. When a physical model for accreting black holes is included, the resulting feedback can drive powerful winds that contribute significantly to the amount and metallicity of hot gas, both of which increase the X-ray luminosity. In terms of their stellar kinematics and structural properties, the merger remnants in our simulations resemble elliptical galaxies. We find that the X-ray luminosities of the remnants with B-band luminosities in the range L_B ~ 10^10 - 10^11 Lsun are consistent with observations, while remnants with smaller or larger masses are underluminous in X-rays. Moreover, because the majority of the merger remnants are broadly consistent with the observed scaling relations between temperature, B-band luminosity and X-ray luminosity we conclude that major mergers are a viable mechanism for producing the X-ray halos of large, luminous elliptical galaxies.

We study the effects of feedback from active galactic nuclei (AGN) on emission from molecular gas in galaxy mergers by combining hydrodynamic simulations which include black holes with a three-dimensional, non-local thermodynamic equilibrium (LTE) radiative transfer code. We find that molecular clouds entrained in AGN winds produce an extended CO morphology with significant off-nuclear emission, which may be detectable via contour mapping. Furthermore, kinematic signatures of these molecular outflows are visible in emission line profiles when the outflow has a large line of sight velocity. Our results can help interpret current and upcoming observations of luminous infrared galaxies, as well as provide a detailed test of subresolution prescriptions for supermassive black hole growth in galaxy-scale hydrodynamic simulations.

(Abridged) The violent hierarchical nature of the LCDM cosmology poses serious difficulties for the formation of disk galaxies. To help resolve these issues, we describe a new, merger-driven scenario for the cosmological formation of disk galaxies at high redshifts that supplements the standard model based on dissipational collapse.In this picture, large gaseous disks may be produced from high-angular momentum mergers of systems that are gas-dominated, i.e. M_gas/(M_gas +M_star > 0.5 at the height of the merger. Pressurization from the multiphase structure of the interstellar medium prevents the complete conversion of gas into stars during the merger, and if enough gas remains to form a disk, the remnant eventually resembles a disk galaxy. We perform numerical simulations of galaxy mergers to study how supernovae feedback strength, supermassive black hole growth and feedback, progenitor gas fraction, merger mass-ratio, and orbital geometry impact the formation of remnant disks. We find that disks can build angular momentum through mergers and the degree of rotational support of the baryons in the merger remnant is primarily related to feedback processes associated with star formation. Disk-dominated remnants are restricted to form in mergers that are gas-dominated at the time of final coalescence and gas-dominated mergers typically require extreme progenitor gas fractions (>80%). We also show that the formation of rotationally-supported stellar systems in mergers is not restricted to idealized orbits, or major or minor mergers. We suggest that the hierarchical nature of the LCDM cosmology and the physics of the interstellar gas may act together to form spiral galaxies by building the angular momentum of disks through early, gas-dominated mergers.

We investigate the effect of embedded active galactic nuclei (AGN) in galaxy mergers on the CO molecular line emission by combining non-local thermodynamic equilibrium (LTE) radiative transfer calculations with hydrodynamic simulations. We find that AGN feedback energy in gas rich galaxy mergers can contribute to large molecular outflows which may be detectable via velocity-integrated emission contour maps, as well as through kinematic features in the emission line profiles.

Galactic sized gravitational lenses are simulated by combining a cosmological N-body simulation and models for the baryonic component of the galaxy. The lens caustics, critical curves, image locations and magnification ratios are calculated by ray-shooting on an adaptive grid. When the source is near a cusp in a smooth lens' caustic the sum of the magnifications of the three closest images should be close to zero. It is found that in the observed cases this sum is generally too large to be consistent with the simulations implying that there is not enough substructure in the simulations. This suggests that other factors play an important role. These may include limited numerical resolution, lensing by structure outside the halo, selection bias and the possibility that a randomly selected galaxy halo may be more irregular, for example due to recent mergers, than the isolated halo used in this study. It is also shown that, with the level of substructure computed from the N-body simulations, the image magnifications of the Einstein cross type lenses are very weak functions of source size up to $\sim 1\kpc$. This is also true for the magnification ratios of widely separated images in the fold and cusp caustic lenses. This means that selected magnification ratios for different the emission regions of a lensed quasar should agree with each other, barring microlensing by stars. The source size dependence of the magnification ratio between the closest pair of images is more sensitive to substructure.

The kinematical properties of elliptical galaxies formed during the mergers of equal mass, stars+gas+dark matter spiral galaxies are compared to the observed low velocity dispersions found for planetary nebulae on the outskirts of ellipticals, which have been interpreted as pointing to a lack of dark matter in ellipticals (which poses a problem for the standard model of galaxy formation). We find that the velocity dispersion profiles of the stars in the simulated ellipticals match well the observed ones. The low outer stellar velocity dispersions are mainly caused by the radial orbits of the outermost stars, which, for a given binding energy must have low angular momentum to reach their large radial distances, usually driven out along tidal tails.

Elements of kinematical and dynamical modeling of elliptical galaxies (Es) are presented. In projection, NFW models resemble Sersic models, but with a very narrow range of shapes (m=3+/-1). The total density profile of Es cannot be NFW-like because the predicted local M/L and aperture velocity dispersion within an effective radius (Re) are much lower than observed. Stars must then dominate Es out to a few Re. Fitting an NFW model to the total density profile of Sersic+NFW (stars+dark matter [DM]) Es results in very high concentration parameters, as found by X-ray observers. Kinematical modeling of Es assuming an isotropic NFW DM model underestimates M/L at the virial radius by a factor of 1.6 to 2.4, because dissipationless LCDM halos have slightly different density profiles and slightly radial velocity anisotropy. In N-body+gas simulations of Es as merger remnants of spirals embedded in DM halos, the slope of the DM density profile is steeper when the initial spiral galaxies are gas-rich. The Hansen & Moore (2006) relation between anisotropy and the slope of the density profile breaks down for gas and DM, but the stars follow an analogous relation with slightly less radial anisotropies for a given density slope. Using kurtosis (h_4) to infer anisotropy in Es is dangerous, as h_4 is also sensitive to small levels of rotation. The stationary Jeans equation provides accurate masses out to 8 Re. The discrepancy between the modeling of Romanowsky et al. (2003), indicating a dearth of DM in Es, and the simulations analyzed by Dekel et al. (2005), which match the spectroscopic observations of Es, is partly due to radial anisotropy and to observing oblate Es face-on. However, one of the 15 solutions to the orbit modeling of Romanowsky et al. is found to have an amount and concentration of DM consistent with LCDM predictions.

(ABRIDGED) We examine the fundamental scaling relations of elliptical galaxies formed through mergers. Using hundreds of simulations to judge the impact of progenitor galaxy properties on merger remnants, we find that gas dissipation provides an important contribution to tilt in the Fundamental Plane relation. Dissipationless mergers of disks produce remnants that occupy the virial plane. As the gas content of disk galaxies is increased, the tilt of the Fundamental Plane relation increases and the slope of the Re-M_* relation steepens. For gas fractions fgas > 30%, the simulated Fundamental Plane scalings approach those observed in the K-band. In our simulations, feedback from supermassive black hole growth has only a minor influence on the stellar-mass scaling relations of spheroidal galaxies, but may play a role in maintaining the observed Fundamental Plane tilt at optical wavelengths by suppressing residual star formation in merger remnants. We estimate that \approx 40-100% of the Fundamental Plane tilt induced by structural properties owes to trends in the central total-to-stellar mass ratio M_total/M_* produced by dissipation. Lower mass systems obtain greater phase- space densities than higher mass systems, producing a galaxy mass-dependent central M_total/M_* and a corresponding tilt in the Fundamental Plane.

(Abridged) We examine the evolution of the black hole mass - stellar velocity dispersion (M-sigma) relation over cosmic time using simulations of galaxy mergers that include feedback from supermassive black hole growth. We consider mergers of galaxies varying the properties of the progenitors to match those expected at redshifts z=0-6. We find that the slope of the resulting M-sigma relation is the same at all redshifts considered. For the same feedback efficiency that reproduces the observed amplitude of the M-sigma relation at z=0, there is a weak redshift-dependence to the normalization that results from an increasing velocity dispersion for a given galactic stellar mass. We develop a formalism to connect redshift evolution in the M-sigma relation to the scatter in the local relation at z=0. We show that the scatter in the local relation places severe constraints on the redshift evolution of both the normalization and slope of the M-sigma relation. Furthermore, we demonstrate that cosmic downsizing introduces a black hole mass-dependent dispersion in the M-sigma relation and that the skewness of the distribution about the locally observed M-sigma relation is sensitive to redshift evolution in the normalization and slope. In principle, these various diagnostics provide a method for differentiating between theories for producing the M-sigma relation. In agreement with existing constraints, our simulations imply that hierarchical structure formation should produce the relation with small intrinsic scatter.

We present an evolutionary model for starbursts, quasars, and spheroidal galaxies in which mergers between gas-rich galaxies drive nuclear inflows of gas, producing intense starbursts and feeding the buried growth of supermassive black holes (BHs) until feedback expels gas and renders a briefly visible optical quasar. The quasar lifetime and obscuring column density depend on both the instantaneous and peak luminosity of the quasar, and we determine this dependence using a large set of simulations of galaxy mergers varying host galaxy properties, orbital geometry, and gas physics. We use these fits to deconvolve observed quasar luminosity functions (LFs) and obtain the evolution of the formation rate of quasars with a certain peak luminosity, n(L_peak,z). Quasars spend extended periods of time at luminosities well below peak, and so n(L_peak) has a maximum corresponding to the 'break' in the observed LF, falling off at both brighter and fainter luminosities. From n(L_peak) and our simulation results, we obtain self-consistent fits to hard and soft X-ray and optical quasar LFs and predict many observables, including: column density distributions of optical and X-ray samples, the LF of broad-line quasars in X-ray samples and the broad-line fraction as a function of luminosity, active BH mass functions, the distribution of Eddington ratios at z~0-2, the z=0 mass function of relic BHs and total mass density of BHs, and the cosmic X-ray background. In every case, our predictions agree well with observed estimates, and unlike previous modeling attempts, we are able to reproduce them without invoking any ad hoc assumptions about source properties or distributions. We provide a library of Monte Carlo realizations of our models for comparison with observations. (Abridged)

(Abridged) Based on numerical simulations of galaxy mergers that incorporate black hole (BH) growth, we predict the faint end slope of the quasar luminosity function (QLF) and its evolution with redshift. Our simulations have yielded a new model for quasar lifetimes where the lifetime depends on both the instantaneous and peak quasar luminosities. This motivates a new interpretation of the QLF in which the bright end consists of quasars radiating at nearly their peak luminosities, but the faint end is mostly made up of quasars in less luminous phases of evolution. The faint-end QLF slope is then determined by the faint-end slope of the quasar lifetime for quasars with peak luminosities near the observed break. We determine this slope from the quasar lifetime as a function of peak luminosity, based on a large set of simulations spanning a wide variety of host galaxy, merger, BH, and ISM gas properties. Brighter peak luminosity (higher BH mass) systems undergo more violent evolution, and expel and heat gas more rapidly in the final stages of quasar evolution, resulting in a flatter faint-end slope (as these objects fall below the observed break in the QLF more rapidly). Therefore, as the QLF break luminosity moves to higher luminosities with increasing redshift, implying a larger typical quasar peak luminosity, the faint-end QLF slope flattens. From the quasar lifetime as a function of peak luminosity and this interpretation of the QLF, we predict the faint-end QLF slope and its evolution with redshift in good agreement with observations. Although BHs grow anti-hierarchically (with lower-mass BHs formed primarily at lower redshifts), the observed change in slope and differential or luminosity dependent density evolution in the QLF is completely determined by the luminosity-dependent quasar lifetime and physics of quasar feedback.

We investigate the luminosity dependence of quasar clustering, inspired by numerical simulations of galaxy mergers that incorporate black hole growth. These simulations have motivated a new interpretation of the quasar luminosity function. In this picture, the bright end of the quasar luminosity function consists of quasars radiating nearly at their peak luminosities, while the faint end consists mainly of very similar sources, but at dimmer phases in their evolution. We combine this model with the statistics of dark matter halos that host quasar activity. We find that, since bright and faint quasars are mostly similar sources seen in different evolutionary stages, a broad range in quasar luminosities corresponds to only a narrow range in the masses of quasar host halos. On average, bright and faint quasars reside in similar host halos. Consequently, we argue that quasar clustering should depend only weakly on luminosity. This prediction is in qualitative agreement with recent measurements of the luminosity dependence of the quasar correlation function (Croom et al. 2005) and the galaxy-quasar cross-correlation function (Adelberger & Steidel 2005). Future precision clustering measurements from SDSS and 2dF, spanning a large range in luminosity, should provide a strong test of our model.

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

We consider implications of our new model of quasar lifetimes and light curves for the quasar luminosity function (LF) at different frequencies and redshifts. In our picture, quasars evolve rapidly and the lifetime depends on both their instantaneous and peak luminosities. The bright end of the LF traces the peak intrinsic quasar activity, but the faint end consists of quasars which are either undergoing exponential growth to much larger masses and luminosities, or are in sub-Eddington quiescent states going into or coming out of a period of peak activity. The 'break' in the observed LF corresponds directly to the maximum in the intrinsic distribution of peak luminosities, which falls off at both brighter and fainter luminosities. We study this model using simulations of galaxy mergers which successfully reproduce a wide range of observed quasar phenomena, including the observed column density distribution. By combining quasar lifetimes and the distribution of maximum quasar luminosities determined from the observed hard X-ray LF with the corresponding luminosity and host-system dependent column densities, we produce the expected soft X-ray and B-band LFs. Our predictions agree exceptionally well with the observed LFs at all observed luminosities, over the redshift range considered (z < 1), without invoking any ad hoc assumptions about an obscured population of sources. Our results also suggest that observed correlations in hard X-ray samples between the obscured fraction of quasars and luminosity can be explained in the context of our model by the expulsion of surrounding gas due to heating from accretion feedback energy as a quasar nears its peak luminosity and final black hole mass.

Based on numerical simulations of gas-rich galaxy mergers, we discuss a model in which quasar activity is tied to the self-regulated growth of supermassive black holes in galaxies. Nuclear inflow of gas attending a galaxy collision triggers a starburst and feeds black hole growth, but for most of the duration of the starburst, the black hole is heavily obscured by surrounding gas and dust which limits the visibility of the quasar, especially at optical and UV wavelengths. Eventually, feedback energy from accretion heats the gas and expels it in a powerful wind, leaving a 'dead quasar'. Between buried and dead phases there is a window during which the galaxy would be seen as a luminous quasar. Because the black hole mass, radiative output, and distribution of obscuring gas and dust all evolve strongly with time, the duration of this phase of observable quasar activity depends on both the waveband and imposed luminosity threshold. We determine the observed and intrinsic lifetimes as a function of luminosity and frequency, and calculate observable lifetimes ~10 Myr for bright quasars in the optical B-band, in good agreement with empirical estimates and much smaller than the black hole growth timescales ~100 Myr, naturally producing a substantial population of 'buried' quasars. However, observed and intrinsic energy outputs converge in the IR and hard X-ray bands as attenuation becomes weaker and chances of observation greatly increase. We obtain the distribution of column densities along sightlines in which the quasar is seen above a given luminosity, and find that our result agrees remarkably well with observed estimates of the column density distribution from the SDSS for appropriate luminosity thresholds. (Abridged)

We propose a new interpretation of the quasar luminosity function (LF), derived from physically motivated models of quasar lifetimes and light curves. In our picture, quasars evolve rapidly and their lifetime depends on both their instantaneous and peak luminosities. We study this model using simulations of galaxy mergers that successfully reproduce a wide range of observed quasar phenomena. With lifetimes inferred from the simulations, we deconvolve the observed quasar LF from the distribution of peak luminosities, and show that they differ qualitatively, unlike for the simple models of quasar lifetimes used previously. We find that the bright end of the LF traces the intrinsic peak quasar activity, but that the faint end consists of quasars which are either undergoing exponential growth to much larger masses and higher luminosities, or are in sub-Eddington quiescent states going into or coming out of a period of peak activity. The 'break' in the LF corresponds directly to the maximum in the intrinsic distribution of peak luminosities, which falls off at both brighter and fainter luminosities. Our interpretation of the quasar LF provides a physical basis for the nature and slope of the faint-end distribution, as well as the location of the break luminosity.

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

We propose a model of quasar lifetimes in which observational quasar lifetimes and an intrinsic lifetime of rapid accretion are strongly distinguished by the physics of obscuration by surrounding gas and dust. Quasars are powered by gas funneled to galaxy centers, but for a large part of the accretion lifetime are heavily obscured by the large gas densities powering accretion. In this phase, starbursts and black hole growth are fueled but the quasar is buried. Eventually, feedback from accretion energy disperses surrounding gas, creating a window in which the black hole is observable optically as a quasar, until accretion rates drop below those required to maintain a quasar luminosity. We model this process and measure the unobscured and intrinsic quasar lifetimes in a hydrodynamical simulation of a major galaxy merger. The source luminosity is determined from the black hole accretion rate, calculated from local gas properties. We calculate the column density of hydrogen to the source along multiple lines of sight and use these column densities and gas metallicities to determine B-band attenuation of the source. Defining the observable quasar lifetime as the total time with an observed B-band luminosity above some limit L_B,min, we find lifetimes ~10-20 Myr for L_B,min=10^11 L_sun (M_B=-23), in good agreement with observationally determined quasar lifetimes. This is significantly smaller than the intrinsic lifetime ~100 Myr obtained if attenuation is neglected. The ratio of observed to intrinsic lifetime is also strong function of both the limiting luminosity and the observed frequency.

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

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